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Disclinations were first observed in mesomorphic phases. They were later found relevant to a 
number of ill-ordered condensed matter media, involving continuous symmetries or frustrated 
order. Disclinations also appear in polycrystals at the edges of grain boundaries; but they are of 
limited interest in solid single crystals, where they can only move by diffusion climb and, owing 
to their large elastic stresses, mostly appear in close pairs of opposite signs. 

The relaxation mechanisms associated with a disclination in its creation, motion, change of shape, 
involve an interplay with continuous or quantized dislocations and/or continuous disclinations. 
These are attached to the disclinations or are akin to Nye's dislocation densities, which are par- 
ticularly well suited here. We introduce the notion of extended Volterra process which takes these 
relaxation processes into account and covers different situations where this interplay takes place. 
These concepts arc illustrated by some of a variety of applications in amorphous solids, mesomor- 
phic phases and frustrated media in their curved habit space. These often involve disclination 
networks with specific node conditions. 

The powerful topological theory of line defects only considers defects stable against any change 
of boundary conditions or of relaxation processes compatible with the structure considered. It 
can be seen as a simplified case of the approach considered here, particularly suited for media of 
high plasticity or/and complex structures. It cannot analyze the dynamical properties of defects 
nor the elastic constants involved in their static properties; topological stability cannot guarantee 
energetic stability and sometimes cannot distinguish finer details of structure of defects. 
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I. INTRODUCTION 

A. General considerations 

Defects in mesomorphic phases (or liquid crystals) have been the subject of numerous investigations in the last 
thirty years. This research has emphasized the importance of disclinations, line defects first defined by Volterra 
(1907) and the main specific defects of liquid crystals (Fricdcl, 1922). As a consequence, the role of defects of a similar 
nature has also been recognized in other media, most of them with non solid crystalline symmetries, but not only. 
This essay is devoted to the geometrical and topological concepts relevant to this field of research, in various media 
where disclinations are acting in interplay with other defects, mostly other types of line defects (dislocations). 

We are led to compare two different theoretical approaches to the classification of line defects, the Volterra process 
and the topological stability theory. They are not equivalent, but rather complementary. The topological theory has 
attracted attention by its property of classifying, not only line defects, but also defects of any dimensionality, in a very 
general manner, on the basis of the topological properties of the order parameter. The Volterra process only applies 
to line defects, which it classifies by the elements of the symmetry group of the ordered medium. This approach 
allows to deal with the (static and dynamical) interplay between disclinations and other line defects. Continuous 
defects (which relate to continuous symmetries) thus control the shape of quantized disclinations, to which they are 
associated, either attached to them, or accompanying them at a small distance. These specific defects of the Volterra 
process have been little studied in recent years. Yet the sustained interest in mesomorphic and frustrated phases, and 
in other media whose structural properties are remote from 3D crystalline media, justifies a reappraisal of the subject. 

As we shall sec, the Volterra process yields the same main conclusions as the topological stability theory, but at 
a finer level, by properly handling boundary conditions and all the plastic relaxations, including those related to 
line-attached defects. This approach can be particularly useful when dealing with nanostructures or with dynamical 
aspects, when the viscosity is large. 

We do not approach the subject of mesomorphic phases straightaway. The final situation depends of course on 
the symmetries of the media under consideration, a topic that is not visited in the first part (section II) of our 
essay, and on the physical nature of the relaxation processes that bring the defects towards their final, stable or 
metastable, state. We shall consider several structures characteristic of ill-condensed media. But we shall first make 
a detour through isotropic uniform solid media, with the sole purpose of understanding the generic geometrical 
relations between disclinations on one hand, dislocations on the other, and of their interplay A number of new results 
are inserted, whose weight in the general balance of the paper has forced us to limit the review of some topics (the 
foundations of the topological theory, geometrical frustration), especially when excellent reviews already exist. On 
the other hand, the basic ingredients of the theory of continuous defects, in the perspective we want to replace it, 
are rather scattered in the literature, and deserved anyhow some deeper analysis; the distinction between constitutive 
and relaxation line-attached dislocations is new and structures this topic. 

As a major application, the concept of disclination in smectic A (SmA) and other mesomorphic phases will be 
considered at some length. Other topics of some importance will also be tackled: 

- the role of disclinations in polycrystal structure and deformation, of importance for polynanocrystals, 
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FIG. 1 The Volterra process for elementary types of dislocations: (a) screw dislocation, (b) and (c) two constructions of the 
same edge dislocation; and dischnations: (d) wedge disclination, (e) twist disclination. The original cut surface E is a vertical 
half plane bounded by the axis of a cylinder of matter. View before the void is filled ((c) and (d)) and the medium relaxed 
(all). From Kleman and Lavrentovich, Soft Matter Physics, an Introduction, Springer, 2003, with permission 

- the structure of an undcrcoolcd liquid, which features a curved space crystal mapped into flat physical space. This is 
a typical example, along with liquid crystalline Blue Phases (BPs), of geometrical frustration. Our analysis of defects 
in a 3D 'amorphous' curved space S 3 is new and generalizes the results for the usual E 3 amorphous phase. 

Before listing the contents of this paper, we find it useful to recall the main feaures of the description of defects 
in ordered madia in terms of the Volterra process and its extension, the topological classification and the theory of 
continuous defects, all involved in our description of dischnations. 

We recall the main characteristics of the Volterra process for defect lines in a solid body (Volterra, 1907), cf. J. 
Friedel (1964). Cut the matter along a surface £ (the cut surface) bound by a line L (a loop, or an infinite line), 
displace the two lips £' and S" of the cut surface by a relative rigid displacement that can be analyzed as the sum 
of a translation b and a rotation 17, introduce matter in its perfect state {i.e elastically undeformed) in order to fill 
the void left by the rigid displacement, or remove matter in the regions of double covering, glue back together new 
matter and primitive matter along the lips £' and let the medium relax elastically; Fig. 1 represents the Volterra 
process for different orientations of b and *f2 with respect to a straight dislocation line. This process is certainly 
badly defined along the line L itself (an elastic singularity remains along L, the defect) and generates a singularity 
of the order parameter on the cut surface; however this latter difficulty disappears if b and f2 are translational and 
rotational symmetries of the medium {quantized, perfect defects). The strain field is small if the rigid displacement 
is restricted to a translation of small amplitude, comparable to the atomic distance say; L is then a dislocation, b its 
Burgers vector. There always are large amplitude displacements if the Volterra process applies to a rotation f2 = f2t, 
except on the axis of rotation t itself; L is then a disclination. This is the reason why research in the field of defects 
presents such a complexity and, thereby, such an accumulation of riches. 

When the Volterra process involves at the same time a translation and a rotation, one has a dispiration (Harris, 
1970, 1977) 1 . 



Dispirations have been studied, but to a limited extent, only in some smectic phases where some symmetry elements include simulta- 
neously a translation and a rotation. Tanakashi et al. (1992) have observed dispirations in a antiferroelectric smectic phase (ShiCa): 
the layer thickness do is half the repeat distance of the polarization P, which changes sign from one layer to the next. Hence a do 
translation and a n rotation of P constitute together a helical symmetry. See also Kuczinski et al. (1999) for observations in a chiral 
antiferroelectric smectic phase (SmQ), and Lcjcck (2002)for theoretical considerations, as well as references therein. 
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a. Quantized perfect disclinations. 

The situation is particularly simple if the line is straight, along the rotation vector (wedge disclination) , more 
complex in the Volterra process sense when the rotation vector is perpendicular to the line (twist disclination), in 
which case it generally involves the simultaneous presence of perfect dislocations attached to the twist line. Research 
in this domain is very lively: the approach in the Volterra process terms is here supplemented by the theory of 
topological stability, see LB below. 

Stemming from the Volterra process picture of a defect, line and surface geometry theory is specifically employed for 
columnar and lamellar media (Klcman et al. , 2004) and, more generally, for liquid media with quantized translations 
(Achard et al., 2005). Riemannian geometry is the main mathematical tool to treat perfect quantized disclinations 
in geometrically frustrated media (Klcman, 1989). 

b. Three important concepts in the development of the Volterra process. 

(i) Imperfect line singularities. If the relative displacement (f2, b) is not a symmetry of the medium, the defect line 
is bordering a surface of misfit along the cut surface of the Volterra process. One will talk of imperfect dislocation, 
disclination, or dispiration.This is for instance the case, in a crystal, of partial dislocations bordering a stacking fault 
or of disclinations bordering grain boundaries. 

(ii) Continuous distributions of defect lines of infinitesimal strength (ft, b), already considered by Volterra (1907). 
Such planar distributions were introduced by Frank (1950b) for describing pure flexion and rotation grain boundaries, 
which can be attached to straight wedge and twist disclinations respectively (Q respectively parallel and perpendicular 
to the line). Three dimensional distributions were introduced by Nye (1953) and others to describe plastic distortions 
of minimal energy. 

(iii) Plastic relaxation, leading to extended Volterra process. In a medium such as liquid crystals where some of 
the stresses can be released plastically in a viscous way, this relaxation plays a large role in reducing the energy 
and increasing the mobility and flexibility of disclinations. The classical Volterra process, which refers to a solid 
(frozen) medium, has then to be completed by a stress relaxation that can be analyzed in terms of a three dimensional 
continuous distribution of dislocations of infinitesimal strength. This is usually obtained by replacing the clastic stress 
field of a solid medium by the elastic stresses referring to the liquid crystal considered (e.g. Frank and Oseen equations 
(Frank, 1958) replacing Hooke's elasticity). We will in this case talk of an extended Volterra process. 

Somewhat similar plastic stress relaxations are possible in solids at the tip of slip lines or cracks, considered as 
dislocations or disclinations; the end of a sub-boundary produced by slip can also relax plastically by developing a 
localized crack. One can also talk in those cases of extended Volterra process, but remember that the plastic relaxation 
can have now a finite elastic limit or at least a strongly viscous friction (Fricdel, 1959b). 

The motion and bending of disclinations involve the production of two dimensional continuous distributions of 
infinitesimal dislocations that can often disperse in a liquid crystal, or be plastically relaxed in a solid, or again be 
absorbed by another defect, as we shall show on some examples. 

B. The topological classification 

a. Why a topological classification 

We outline the general principles of this theory very briefly, and illustrate them on some simple examples of liquid 
crystal phases, in order to bring out the concept of topological stability (the essential contribution of this theory to 
the physics of defects), and to comment on the effects of non-commutativity of the symmetries, which are taken into 
account in a systematic way in this theory. 

The topological theory started with the articles of Rogula (1976); Toulouse and Klcman (1976); Volovik and Mineev 
(1976); for recent reviews, with emphasis on mesomorphic phases, see Mermin (1979); Michel (1980); Trcbin (1982). 
The topological classification has cured a certain number of difficulties raised by the Volterra process, when applied to 
ill-condensed matter, mostly liquid crystals. Liquid crystal defect features have no equivalent in usual solid crystals: 
disclinations whose core singularity vanishes away, point singularities, 3D knotted non-singular configurations. 

- the Volterra process correctly describes straight wedge disclinations of strength \k\ — 1/2 (i.e., rotational symme- 
tries of angle u> = ±7r), but cannot be extended to |fc| > 1/2, i.e., angles u> = (2n ± l)ir, n ^ 0, n € Z, in its 'naive' 
version, 

- twist disclinations cannot be constructed, except locally (for an illustration, see (Harris, 1977)), 

- the 'escape in the 3rd dimension' (Meyer, 1973) obviously does not come out of the Volterra process. 
However these difficulties can also be dealt with successfully by introducing the concept of extended Volterra process. 

b. The order parameter space 
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The topological classification of defects relies on the application to ordered media of the methods and concepts of 
algebraic topology; standard references are Steenrod (1957) or Massey (1967). The manifold V of internal states, also 
called the vacuum or the order parameter space is the space of all the possible different positions and orientations of the 
perfect crystal in (flat) Euclidean space. Let H be the group of symmetry of the ordered structure, G=E 3 =R 3 DO(3) 
the group of Euclidean isometrics of E 3 . The symbol □ is for the semi-direct product of groups. R 3 is the 3D group of 
continuous translations, 0(3)=SO(3) xZ 2 the full group of rotations, with center of symmetry included. The symbol 
x is for the direct product of groups, H is a subgroup of G. V is then the quotient space G/H. Observe that V is not 
a group, generically, except if H is a normal subgroup of G. 

Examples of order parameter spaces are given in Kleman and Lavrentovich (2000); the order parameter space of 
a uniaxial nematic V(N) = P , the projective plane; the order parameter space of a 3D crystal, regarding only the 
translations, is V(Xal)=T 3 , the 3D torus. 

These concepts extend in a natural way when G is the group of a space M of constant but ^ curvature, and 
H a subgroup of G, i.e. the group of symmetry of a ordered structure of habit space M. We shall make use of this 
extension in the investigation of defects in frustrated media, see section VII. 

c. The Erst homotopy group (the Fundamental Group ) of the order parameter space 

Let us now consider how the order parameter space V enters into the issue of the topological classification of defects. 
Start from a distorted ordered medium, in which the order parameter is broken along a line L. In order to test the 
topological nature of the breaking, surround L by a closed loop 7 entirely located in the good ordered structure, good 
in the sense of the usual theory of dislocations in solids. It is possible to attach to each point r belonging to 7 a 
'tangent' perfect ordered structure, which maps in a unique way onto some point R G V. When r traverses the closed 
loop 7, R traverses a closed loop V on V. Call (f> this well-defined continuous mapping 

cf) : 7-* r. 

The function (f> can be extended continuously to the whole continuous domain D of the ordered medium in which 
the order parameter is well defined, since the order parameter is expected to vary continuously in D. Therefore 
any continuous displacement of 7 in D maps on a continuous displacement of T. There is consequently a relation 
of equivalence between the different images Ts; it is this equivalence that is described by the notion of homotopy, 
(Steenrod, 1957). All the Ts belonging to the same class of equivalence are represented by an element [r] of a 
group, the so-called fundamental group (or first homotopy group) 111(1/). More precisely, ni(V) is the group of 
classes of oriented loops belonging to V, equivalent under homotopy, and all having the same base point. This latter 
technicality has no incidence on the classification of defects. More important is the remark that, in most cases, the 
fundamental group is not commutative, a property that is related to the fact that the topological charge of a defect, 
i.e. the corresponding clement in 111(1/), is modified when the defect [T] has circumnavigated about a defect [V]; 
it is changed to [V] [T] [T']~ , an element of the same conjugacy class. It is therefore usual to consider that all the 
elements of a given class of conjugacy represent the same defect. Examples are provided later on in this article. 

C. The theory of continuous defects 

This theory flourished in the early sixties, see Nye (1953); Bilby et al. (1955); Kondo (1955-1967); Kroner (1981); 
it was at that time applied to solid crystals only. It concentrates mostly on the study of sets of line defects, whether 
these defects, of the quantized type, are considered at such a scale that the concept of defect density makes sense, 
or whether the characteristic invariants carried by the defects (in the sense of Volterra, i.e. translations, rotations) 
are continuous, with the result that the notion of infinitesimal defects — with vanishing Burgers vector dislocation or 
vanishing rotation vector disclination — is significant. And indeed most of the results of the theory of continuous defects 
relates to continuous distributions of dislocations of infinitesimal strength, to intrinsic point defects (interstitials, 
vacancies, etc), and, to a lesser extent, to disclinations of infinitesimal strength. The main geometrical ingredients of 
the continuous theory of dislocation line defects and intrinsic point defects are the concepts of torsion and curvature 
on a manifold; the points of contact with the Einstein's theory of generalized relativity are therefore numerous, 
and have been stated a number of times (Hchl et al., 1976; Kroner, 1981): cosmic strings are spacetime topological 
defects that may be described in delta- function valued torsion and curvature components, carried respectively by 
translation and rotation symmetry breaking defects in a Minkowskian manifold (Vilenkin and Shcllard, 1994); for 
recent spacetime defect calculations see e.g. Letcher (1995), Puntigam and Soleng (1997), and references therein. 
Defects such as disclinations whose characteristic invariants belong to non-commutative groups cannot so easily be 
turned into density sets. Methods borrowed from the theory of fields in high-energy physics have also been applied; 
this is the so-called gauge field theory of defects, cf. Julia and Toulouse (1979), Dzyaloshinskii and Volovik (1980), 
Dzyaloshinskii (1981). 
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It was expected that these continuous approaches would open a way to a certain number of problems difficult to 
attack by the physics of quantized defects when these defects are plenty. Thus the dynamical theory of dislocations 
in solids, the defect melting theory, the defect content of deformed frustrated phases (e.g. Frank and Kasper phases, 
quasicrystals, even amorphous media and glasses, and among mesomorphic blue phases or twist grain boundary (TGB) 
phases), also would provide a suitable definition of other defect densities: of disclinations, of point defects, etc. The 
rich and interesting courses given at the Les Houches Summer School on Defects held in 1980 (Dzyaloshinskii, 1981; 
Kroner, 1981) take stock of the various advances in continuous and gauge field theory made at the time. 

However, although the continuous theory of defects is of a rare mathematical elegance, the applications have been 
scarce, one of the most convincing being perhaps the analysis of magnetostrictive effects in ferromagnets (Kleman, 
1967). At the Discussion Meeting he organized in Stuttgart (Kroner, 1982), Kroner made the following remark: 
"Although the field theory of defects [by field theory he obviously meant the traditional theory of continuous defects as 
well as its gauge field extension] has found many applications, the early hope that it could become the basis of a general 
theory of plasticity has not been fulfilled. Among various reasons we mention first of all, that the defects namely the 
dislocations that above all are responsible for plastic flow, do not form smooth line densities that can well be described by a 
dislocation density tensor field. Direct observation of dislocations in crystals, for instance by means of electron microscopy, 
shows that dislocations rather form three dimensional networks [our sub-lining] that are interconnected in practically 
immobile nodes and other often complex local arrangements... These networks have a strong statistical component, a fact 
that shows that a real physical understanding of plasticity requires also considerations in the frame of statistical physics. 
However, a statistical theory of interacting deformable lines that can be created, annihilated and change their length, has 
never been worked out" . 

To these remarks can be added that the material science physicists, who certainly know best the problems and the 
traps of the physics of defects, are in general not fully aware of the mathematical tools (non-Riemannian manifolds, 
exterior calculus, Grassmann algebra, differentiation on manifolds and fibre bundles) that are at the very basis of 
continuous gauge theory. This language problem has little chance to be solved in the near future, inasmuch as the 
gap between material scientists and field theorists keeps widening. By all means, we here avoid as much as possible 
complex mathematical tools. 

In any case, we arc led to conclude that field theory is not a panacea. On the other hand, as a matter of fact, the 
fundamentals of plastic deformation and fracture of crystalline materials have recently undergone a revival through 
the development of new experimental methods that now explore nanometric scales, and through the improvement of 
computer power and computing methods. 

Dislocations (whose set of characteristic invariants is isomorphic to an infinite abelian group) and intrinsic point 
defects are the dominant defects in solid crystals; there thereby is no real necessity to introduce in the theory the 
group-theoretic description of other types of defects that would require finite abelian groups or non-abelian groups, 
making the theory unmanageable. Notwithstanding this simplification, Kroner's criticism stands valid; even if one 
restricts to static situations, the continuous theory neglects real dislocation networks — in the sense of Frank (1950b) 
— , and the frictional effects due to them. This is probably the true reason why the theory of continuous dislocations 
has found so little use yet. 

D. Disclinations 

a. Disclinations and continuous dislocations 

There is however a situation where dislocation densities retrieve their true importance and where continuous theory 
plays a role, it is with disclinations considered as singularities of dislocation and disclination densities; this is the point 
of view which is taken here, and which has some relevance to mesomorphic phases and perhaps also to quasicrystals, 
Frank and Kasper phases, undercooled liquids, polynanocrystals, and more generally frustrated phases. 

Disclinations can exist in solid crystals, whose building blocks are atomic, point like; but, as already stated, 
the continuous theory has been but little applied to disclinations in solid crystals, where such objects have a 
huge energy. On the other hand, disclinations are the rule in mesomorphic phases, whose building blocks are 
anisometric molecules (rod like, disk like, etc). Again, these disclinations quite often appear as the singularity set 
of a dislocation density. Hence the interest in revisiting the continuous theory of defects, although it appears at 
once that new concepts have to emerge. The case of mesomorphic phases requires an extension of the theory of 
continuous defects for solids to the case when there is locally only one physical direction (the director), i.e. no 
local trihedron of directions, as in the uniaxial nematic N, the SmA and the columnar D cases 2 ; see also a re- 



2 The continuous theory of defects makes use of lattice manifolds, whose points carry local trihedra. 
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lated remark in IV. A. 2. The role of stress relaxation is especially important and complex in such mesomorphic phases. 

b. Three dimensional networks 

The question also arises whether disclinations form 3D networks in amorphous systems, liquid crystals such as 
cholcstcric Blue Phases, undcrcoolcd liquids — as already stated, it has been hypothesized that glass disorder can be 
described in terms of disclinations in a icosahedral curved crystal (Kleman and Sadoc, 1979) — or in nanocrystals, 
clusters (Fricdcl, 1984) and quasicrystals, where they are indeed documented (Frank and Kaspcr, 1958, 1959; Nelson, 
1983a). In such 3D networks, the disclinations have to be somewhat flexible, which is possible, whether these 
disclinations are quantized or not, only if other defects, dislocations or disclinations, continuous or not, attach to 
them. Thus the question of the interplay between dislocations and disclinations goes beyond mesomorphic phases. 
This question will be presented and discussed in this essay. 

E. Outline of the article 

To summarize, the continuous theory of defects, in its primitive form, only considers dislocation densities, which 
are singularities of continuous fields. It does not consider finite defects like disclinations, neither grain boundaries, nor 
Frank networks. This essay, contrariwise, assumes the coexistence of finite and infinitesimal defects. Grain boundaries 
are introduced. The field theory instruments are not employed; but this is possible (see Kleman (1982b), in a similar 
context, but in reason of the advances presented here, the results have to be revisited). An important aspect of 
disclinations in mesomorphic phases is how their flexibility and mobility concur to the relaxation of stresses imposed 
by the boundary conditions (static or dynamic). Both dislocation and disclination densities, we shall advance, play a 
leading role in such relaxation processes. 

Section II is about the description, in geometrical terms, of the defect structure of disclinations, without taking into 
account the constraints due to the symmetries of the medium. Therefore this section applies to amorphous media 
and isotropic liquids, but at this stage it does not provide more than the geometrical tools to study the role of stress 
relaxation mediated by continuous defects in various media, and how it affects disclination line flexibility and mobility, 
and the geometrical rules to building disclination networks. It is therefore directly applicable to solid media only. 

The results of section II are employed in section III to shed a new light on the properties of nanocrystals, and are 
extended in section IV to quantized disclinations, examples being taken in mesomorphic phases. This is the place 
where we discuss the relationship between topological stability and the kind of stability that stems from the Volterra 
process. Section IV and section V (more especially devoted to focal conies in SmA's as quantized disclinations), are 
thereby about the nature of disclinations in partly ordered materials, and their interplay with continuous and/or 
quantized dislocations. 

Defects in frustrated phases are discussed in section VI. taking into account the frustrated local order, it extends 
to three dimensional spherical amorphous media the results obtained in section II for three dimensional euclidean 
amorphous media, making full use of the quaternion representation of the geometry of S 3 . Important results are (i) 
that the 'dislocations', like the translational symmetries they break, are non-commutative — we call them disvections 
— , (ii) that infinitesimal disclinations (rather than dislocations) are attached to 'twist' finite disclinations. This 
result emphasizes the role of disclination networks in frustrated media. A part of VI is devoted to the classification 
of defects in the spherical {3,3,5} polytope, which has been used as a template for amorphous media with local 
icosahedral order. Section VLB discusses some characteristics of the decurving process of the just alluded to curved 
media. 

The discussion in section VIII bears on a comparison between the extended Volterra process and the topo- 
logical theory — a question that runs as a red thread through the entire essay — and expatiates on the question 
of plastic relaxation, i.e. the role in various media of continuous — mostly — and quantized defects in stress relaxation. 



II. CONTINUOUS DEFECTS IN ISOTROPIC UNIFORM MEDIA. GEOMETRICAL INTERPLAY BETWEEN 
DISCLINATIONS AND DISLOCATIONS 

An amorphous metal, considered at a scale larger than the atomic size, is an example of an isotropic uniform solid 
medium. The Volterra process allows the consideration of continuous, non-quantized, dislocations and disclinations 
that carry stresses. On the other hand a result of the topological theory of defects is that these are not topologically 
stable. The objects to which this section is devoted are therefore, to the best, metastable. 

In subsections II. A to II. D we investigate disclinations of finite strength whose rotation vector f2 is constant in 
modulus and direction. Such objects are attended by two types of attached infinitesimal dislocations (constitutive 
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FIG. 2 Wedge disclination L and its edge dislocation content. OM = OM~ = OM+; M + M~ = 2QMsin — 



and relaxation dislocations), from which the concept of extended Volterra process. We comment on the equivalence 
between these infinitesimal dislocation sets and grain boundaries in subsection II. E. In subsections II. F and II. G we 
investigate the case when the rotation vector is varying along the line, in direction and in modulus. The ensuing 
considerations directly yield an expression for the fundamental invariant of a disclination, which we call the Frank 
vector. The Frank vector is for disclinations what the Burgers vector is for dislocations; in particular it satisfies a 
Kirchhoff relation at disclination nodes. The detailed study of quantized disclinations is postponed to sections IV and 
V. 

This case of anisotropic uniform medium is chosen for the simplicity of the relaxations in the extended Volterra 
process. As stressed above, it is somewhat artificial to distinguish the relaxations of stresses involved in the motions of 
a given (non quantized) disclination from that for the dispersion of the disclination itself: the same plastic properties 
are involved in both and should occur in similar lengths of time in non viscous liquids. It is only for very viscous 
liquids or better for amorphous solids with slow atomic diffusion that one can reasonably assume the short range 
relaxation of stresses due to the motion of the dislocations to be more rapid than their dispersion life time, which we 
shall consider here as infinite. The dynamical competition of the two processes in somewhat viscous liquids has so far 
not been much studied. 

The models developed and applied in sections II and III will be extended in sections IV and V, with no such 
difficulty, to quantized and topologically stable disclinations. 



A. Dislocation content of a straight wedge disclination 

Let us consider an infinitely long wedge disclination line L, of rotation angle S~2 — fit, the rotation axis t being 
along L, Fig. 2. 

The Volterra process consists in opening a dihedral void of matter (for the sake of simplicity we consider a disclination 
of negative strength). The relative displacement of the lips of the cut surface E, which we take to be a half-plane, at 



a point M of E, is written as equal to sin — t X OM, O being any origin on L; this displacement can also be the 
result of a set of edge dislocations — (Friedel, 1964), chapter 1 — , located uniformly in E, and whose total Burgers 

vector bM is precisely 2 sin — t x OM, for those dislocations lying between the edge of the dihedron and M, i.e. with 
a density 

db m = 2 sin — txdM (1) 

thus producing a tilt boundary of rotation along E. 
Let us give a few examples: 

(i) In an amorphous solid or a glass the angle f2 can take any value; a continuous distribution of dislocations thereby 
yields a continuous wedge disclination. This tilt boundary introduces a mismatch of short range atomic order, which 
can be suppressed if some local atomic diffusion is allowed at short range. 

(ii) In a ferromagnetic solid, the meeting line of several magnetic walls is a continuous wedge disclination whose 
angle il relates directly to the magnetoelastic constants (Klcman, 1974). Again, this disclination can be analyzed in 
terms of continuous dislocations. 
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(iii) N on- quantized wedge disclinations in a crystalline solid are the limits of tilt boundaries, which are split in the 
usual way into finite dislocations parallel to that limit; this is the type considered up to now. f2 can take any value; 
it is tuned by the density of edge dislocations. For small values of f2, the continuous distribution of infinitesimal 
dislocations can also regroup into parallel dislocations of finite strength allowed by the crystal structure (cf III.C 
for a more detailed discussion). In the general case, this is an imperfect disclination with a stacking fault that is 
a tilt boundary. Such disclinations have very large energies, as soon as f2 is finite, in the absence of any plastic 
relaxation. They can nevertheless be produced, for instance when a slip line crosses a low angle grain boundary, 
during plastic deformation of polygonized crystals; they can also be associated in parallel pairs of equal rotation 
strength and of opposite signs. The stress concentrations produced at the cores of the disclinations are often relaxed 
by the development of cracks (Fricdcl, 1964). 

Remark: The Volterra process is properly defined for |fi| < n. Observe indeed that Eq. 1 does not distinguish 
between (i)- = 2tt and f2 = 0, (ii)- between the angles f2 — 2tt — a and 17 = a. This equation puts limits to the 
application of the Volterra process and to the use of Eq. 1. Notice also that it is inconsistent to consider a unique 
Volterra process with an angle |fi| > 2tt, since this requires removing matter (for f2 > 0) or adding matter (for il < 0) 
in at least a full space. Hence, for any angle = (2n + 1)tt + a, \a\ < tt, one has to consider (2n + 1) successive 
applications of the Volterra process, followed by a Volterra process of angle a. 

We now deepen the relationship of a wedge disclination with its accompanying dislocations by first considering the 
displacement of the entire line parallel to itself ( II. B ), secondly displacing a part of the line only (II. C). 

B. Emitted/absorbed dislocations. Constitutive and relaxation dislocations 

It is clear that the energy of the disclination becomes prohibitive if the rotation axis t does stay in place when the 
disclination is displaced; if some plastic relaxation is allowed, t moves to a new position by the emission or absorption 
of a certain number of dislocations. We have to specify the direction of this motion, as the stress about L depends on 
the position of the cut surface S in space, even if the medium is isotropic. 

1. Motion in the cut surface. 

If L is moved in the plane of X by a displacement vector S, the total Burgers vector hg of the emitted/ absorbed 
dislocations is equal to ±2 sin | txif: these are edge dislocations; is the sum total = ^]b e of elementary 
dislocations b e allowed by the symmetries of the medium, Fig. 3. 

In an amorphous solid or a glass the b e s may have any modulus. In an ordered solid, the b e 's have to be equal to 
translation symmetries of the medium. 

Notice that in both cases the emitted/absorbed dislocations contribute to the relaxation of the sample that has 
suffered the displacement of the disclination. In order to illustrate this point in the amorphous, solid, case, observe 
that it would make no sense if the emitted dislocations stay in place in the continuation of the cut surface, because 
this would not modify the strains and stresses previously carried by the medium, before the line has moved. In other 
words, one is led to recognize the existence of two types of dislocation densities: those belonging to the actual cut 
surface of the disclination, which we call constitutive dislocations, and those left away in the wake of the moving 
disclination, which we call relaxation dislocations. In principle, since the stress field attached to a wedge line L, 
measured in a frame of reference attached to L, is independent of the position of L, the relaxation dislocations should 
carry no stress at all at complete relaxation; they are dispersed in the entire space with vanishing Burgers vectors, if 
they are continuous, or vanish at the boundaries of the sample, if they are quantized. 

Of course, a part of the stresses carried by the disclination can also be relaxed by defects that are not attached to the 
line, e.g. infinitesimal dislocations nucleated in the bulk; these are Nye's dislocations, which we study in some detail 
later when we come to layered media (smectics) (IV. B. 2), and in Appendix A. But we shall not consider non-attached 
defect densities in the present section. The final result depends on the material properties of the medium, whether it 
has solid elasticity (amorphous) or viscous behaviour (liquid). In this latter case the relaxation can be complete. 

2. Motion off the cut surface 

If L is translated off the plane of £ by 8, the total Burgers vector of the absorbed dislocations is also equal to 
±2 sin ^ v X S, where u is the new axis of rotation. A new piece of cut surface, parallel to the {v, 8} plane, is created. 
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FIG. 3 Displacement of a wedge disclination from L (a) to L' (b) by emission of dislocations that disperse away 




conslitutivtf edge rittliwutior^ 



b, b, 



FIG. 4 Kink AB linking the two wedge half-lines L~ and L + . The cut surface E is supposed to be on the right of the disclination 
L _ , AB, L + . The upper segments of the dislocations that traverse AB tend to disperse away (plastic relaxation), while keeping 
attached to the constitutive segments on the kink. 



C. Twist component of a disclination 

We now turn our attention to a line L made of three segments, namely two parallel semi-infinite wedge segments 
L~ and L + , joined by a third perpendicular segment AB of small length, called a kink, Fig. 4. We assume that the 
cut surface £ is a plane that contains the three segments. 

The L + segment is the result of a displacement of a part of the entire line parallel to itself by a translation S = AB, 
by emitting or absorbing dislocations. In this process, f2 stays parallel to itself. According to the results above, we 
have, on the S side of AB, constitutive dislocations of total Burgers vector 

b e (AB) = ±2 sin - t x AB (2) 
on the other side of AB, relaxation dislocations of total Burgers vector 



b d (AB) = ±2 sin - txS 



(3) 
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These two quantities being equal, we see that the dislocations cross the segment AB, but behave quite differently 
on either side. 

A classical Volterra process, acting once for all on the cut surface of the disclination line, is not relevant to the 
present geometry, because such a Volterra process can be performed only if the rotation vector is fixed in space. 



D. £1 constant. Generic disclination line 

Again, we restrict our attention to the isotropic case. Consider a curved disclination line: Fig. 5 shows a disclination 
with constant in length and in direction from point to point along L. The angle between L and varies from point 

dP 

to point. Let P and Q = P H — —ds be two points infinitesimally close on L. The Burgers vector of the infinitesimal 

dislocation introduced by the variation of position of fl from P to Q is, by reasoning on the cut surface £ as above, 
equal to dbpQ = c1q(M) — dp(M), where 

d P (M) = 2sin^t x PM, 
d Q (M) = 2sin — t X QM 

are the displacements of the cut surface at M, any point on S, seen respectively from P and Q. Hence: 

dq(M) - d P (M) = -2 sin - t X ~-ds (4) 

This dislocation, which we denote d£pQ, can be thought of as attached to the line at the infinitesimal arc PQ. Of 
course dbpQ has to be a translation allowed by the symmetry of the phase. The shape taken by dl pq resulting from 
plastic relaxation optimizes the energy carried by the disclination. 



1. Two types of dislocation continuous distributions 

The question arises whether any infinitesimal relaxation dislocation d£pQ, attached to L at the arc PQ, crosses the 
line L, and transforms 

on the other side of L into a constitutive infinitesimal dislocation, with the same Burgers vector, as it does in the 
simple case investigated in II. C. The answer is positive and the demonstration is as follows. M being on the cut 
surface £ of L, Eq. 4 is as a result valid on the full area of E. In other words, the infinitesimal dislocation line d£pQ 
with Burgers vector dbpQ has the same cut surface as L, and consequently meets L, is fixed between P and Q, and is 
closed in the manner of L, Fig. 5(a). In fact, with the local rotation vector defined as above, the disclination L is the 
result of creating a density of infinitesimal d£pQ by infinitesimal Volterra processes on the same cut surface. 

Now let us deform these d£pQ dislocation lines (this is an allowed operation), opening them into infinite lines (or 
line segments ending on the boundaries of the sample) in such a way that now they cross the imaginary line L, which 
divides each of them into two semi infinite arcs. This process traces out the bounds of the cut surface £ along L if one 
imposes different types of distribution for the line arcs on both sides of L, yielding different elastic distributions. In 
the spirit of Fig. 4, one can imagine on one side a 2D surface tiled with constitutive dislocations, a kind of generalized 
misorientation boundary, on the other side, relaxation dislocation segments dispersed in space. See Fig. 5(b). 



2. Line tension of twist vs. wedge segments 

Notice that we have realized a dislocation geometry that display two metastable configurations of a set of infinitesi- 
mal dislocations, namely those of the misorientation boundary and those fully dispersed. Therefore the core region of 
a twist disclination line, where these two configurations merge, has a contribution to the total energy that scales as its 
length. The other contribution is the energy of the dislocation lines, essentially that of the constitutive dislocations, 
which scales as the area of the cut surface, i.e., the square of the length of the line. One expects that this contribution 
is larger than the first one. Therefore the line energy per unit length of line is proportional to the length of the 
line. The latter contribution is worth comparing with the energy per unit length of a wedge line, which scales as the 
square of the transverse size: the line tension of a disclination is thus a linear function of their length which reduces 
to a proportionality for wedge ones, whereas, for twist ones, one must add a small constant term due to their relaxed 
dislocations. 
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FIG. 5 (a)- The infinitesimal dislocations dlpQ are constructed assuming first that the cut surface bound by L is common 
to all the alls attached along L. (b)- Then each MpQ is deformed at fixed P, Q; the two parts of dlpQ on both sides of the 
infinitesimal arc PQ do not have the same elastic distribution. L is the limit between the two distribution types 



E. Disclinations and grain boundaries 

As in crystalline solids (see II. A), wedge lines in solid amorphous materials carry a large energy, except in the 
same special circumstances as indicated above. The existence of twist lines or partially twist lines is even less 
probable, and their mobility and change of curvature is certainly negligible, since it would require the climb of the 
attached dislocations, which requires plastic deformation. Hence a caveat: except in the later case of polycrystals, 
the discussion that follows assumes implicitly that there is no restriction to reaching low energy states by plastic 
relaxation; it therefore applies to some sort of amorphous material endowed with a finite viscosity, which operates 
through the existence and mobility of infinitesimal relaxation dislocation densities. One does not expect that such 
processes are possible in a solid crystal. But the comparison between crystals and amorphous media is worth carrying 
out, especially through the parallel concepts of grain boundary and cut surface. 



1. Frank's grain boundary and Friedel's disclination compared 

Equation 1, integrated along a segment MN of the disclination line L, gives the total Burgers vector of the 
dislocations that are attached to any segment MN of L (and lie along its cut surface £_l). 

Ab M N = 2 sin ^ t x MN (5) 

This expression is similar to Frank's formula (Frank, 1950b), for a crystal grain boundary J^gb of angle of misori- 
entation uj = ujs, |s| = 1; Frank's formula yields the total Burgers vector of the (quantized) dislocations that cross 
any segment PQ belonging to Egb- 

A/3 PQ = 2 sin | s x PQ (6) 

This similarity does not come as a surprise, after the foregoing discussion. In particular, one can deduce from 
Frank's approach that any closed line in a grain boundary can be chosen as a disclination line, provided the 'exterior' 
dislocation segments are let to relax. Identifying Eq. 5 and Eq. 6 and assuming M = P, N = Q, we have 

sin ^ t x MN = sin % s x MN (7) 



which yields either 



sin — t — sin — s = (8) 
2 2 w 
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FIG. 6 Burgers vector. A crystal structure has been superimposed to show more clearly the Burgers circuit 0'M'NM"0"00' 
— in the final state (a) and in the initial state (b) — and the rotation fl of the Volterra process 



sin — t sin — s oc MN (9) 
2 2 w 

Assume that a; is a constant vector, i.e. that the cut surface of the disclination line belongs to a unique grain 
boundary. Equation 8, which yields f2 = to, t — s, applies when 17 is a constant vector, which is what we have 
assumed up to now. Equation 9 on the other hand expresses that other possibilities exist, with f2 ^ to, t ^ s, with 
the tangent t to the disclination (MN = rds) belonging to the plane {t, s}. Although a; is a constant vector, Eq. 9 
applies to a variable S7, a situation more fully discussed in II. F and II. G. We give later on an important example 
obeying Eq. 9 (see Appendix C). 

Remark: The 2 sin — (or 2 sin — ) factor in Eq. 1, 5, or 6 deserves some comments. 

The integral Burgers vector b of the cut S is counted as a lack of closure of a Burgers circuit in the final state 
(more properly in the final state of the Volterra process before elastic relaxation of the sector M'OM" here introduced, 
Fig. 6a). 

It is usual for dislocations to count the Burgers vector bo as a lack of closure in the initial state; one would then 
have (Fig. 6) 

n 

|b| = M'M" = 2 0M'sin- 

and 

|b | =00 ! = 20M' tan-. 

Such a difference originates from the non-commutativity of the Burgers circuit and the rotation of the disclination. 
Obviously, the formula in 2 sin — is to be preferred as it correctly describes the final state, the only one of interest 

here. The difference is only noticeable for large fis, where the Burgers vector b of the constitutive dislocation is 
smaller by a factor cos ^ than that bo of a crystal dislocation. In the sequel we use the Burgers vector formula b; it 
appears as follows in the Frank vector introduced in II. F. 2: 



{(ft) = 2 sin^ t. (10) 
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FIG. 7 An elementary disclination of mixed character. The constitutive dislocations (inside) are drawn, but not the relaxation 
dislocations (outside) 

2. Isolated twist segments 

Equation 6 has been derived by Frank for a finite segment PQ, independently of any loop to which PQ could 
pertain; the dislocations which cross this segment (the crossing set) in the grain boundary have two parts (relaxation 
on one side of PQ, constitutive on the other), with no violation of any conservation law on Burgers vectors. Such 
a segment and its set of attached dislocations constitute a sort of stripe in the grain boundary, geometrically and 
topologically independent of it (certainly not energetically, but we shall pay no attention to this question) . 

Therefore, proceeding with our analysis, we conclude that a continuous disclination line can be made of a sequence 
of independent segments M,Mj +1 , possibly infinitesimal, each of them carrying a different rotation vector fij 
Each stripe MjM^i defines a grain boundary of finite width. Therefore, any isolated, finite, segment MN can divide 
each of the dislocations (continuous or quantized) belonging to the crossing set into constitutive (on one side) and 
relaxation (on the other side) dislocation segments. Each stripe has two edges parallel to f2 and two edges of mixed 
character, see Fig. 7. 

A stripe is an elementary disclination, with a twist (or mixed, but not pure wedge) segment transversal to the 
stripe. The longitudinal boundaries of the stripe, along the dislocations that construct the stripe, can be considered 
as wedge segments. Thus, a stripe, considered as a disclination loop, is necessarily of mixed character. 

F. Polygonal disclination lines. Attached disclinations 

Disclinations with a rotation vector varying in length and direction are possible; they require attached disclinations 
(or attached disclination densities). For the sake of clarity we do not introduce disclination densities straightaway 
and develop the theory for attached disclinations of finite strength. A disclination can be thought of as the sum 
of infinitesimal stripes that tile its cut surface and partition it into long stripes, elongated along the constitutive 
dislocations; the situation is reminiscent of the tiling of the cut surface of a dislocation into elementary dislocation 
loops. 

1. Wedge polygonal loops and bisecting disclination lines 

Consider for instance Fig. 8(a), which represents a disclination line made of two semi-infinite wedge segments L and 
L'. The stripes divide into two sets parallel to the disclination segments; the continuity of the constitutive dislocations 
is insured, when f2' = i~J, if the twist edges of the stripes are along the line that bisects L and L'; it is indeed easy to 
show that the Burgers vector of the constitutive dislocation segments parallel to L and L' are continuous across this 
line. 

A wedge loop is nothing else than a continuous generalization of Fig. 8(a), applied to a closed polygonal disclination; 
is constant in length and everywhere tangent to the loop, the constitutive dislocations close into loops entirely 
located in the cut surface. 

The bisecting line just introduced has a special stability because it has no relaxation dislocations attached to it. It 
is indeed a wedge disclination, as established by the analysis that follows. 



1G 




constitutive dislocations (Burgers 
vectors perpendicular to Q r Q us ) 




FIG. 8 A polygonal disclination; (a) two semi-infinite wedge lines L and L' meeting at O; (b) polygonal disclination made of 
segments Li of mixed character (twist-wedge); the Burgers vectors are continuous across the line parallel to fi 



2. Disclinations meeting at a node. Kirchhoff relation. Frank vector 



The foregoing considerations generalize to a polygonal disclination made of segments Li of mixed character (twist- 
wedge), with r?i varying in direction but also in modulus (S~2 j =/= Qj), Fig- 8(b). The considerations that follow apply 
to a situation where there are no restrictions on tti, ttj, and the directions of the segments Li, Lf, they are not 
necessarily in the same plane. Let Li(S~2,;), Li + i(S~2i + i) be two consecutive disclination segments. The constitutive 
dislocation segments meet without discontinuity of the Burgers vector on the half line parallel to the direction 

fi = 2 sin y U - 2 sin -yi t i+ i , (11) 

as, following Eq. 10, the density of dislocations must be counted along the common edge Lj of the two boundaries. 
In effect, if 0,P is parallel to fi , we have 

2 sin y U x O t P = 2 sin y±i t i+1 x O t P (12) 

Equation 12 means, according to Frank's formula, that the Burgers vector is continuous across any segment parallel 
to fj . The half line Li is obviously a generalization of the bisecting line introduced just above. Notice that in general 
the planes {Li, fi} and {L^+i, fi+i} are not tilt planes for the constitutive dislocations; this happens only if f2.;, iT;+i, 
and the directions of the segments Li, L.; + i are all four coplanar directions. 

a. Kirchhoff relation. Let us introduce, instead of the rotation vector S~2^ = f2jtj, the expression 

fi = 2 sin y ti (13) 



Equation 11 then takes the form 

, . 7 i o ■ ^i+ 

Z Sill * 

which can be written 



2 sin y ti = 2 sm y tj + 2 sin — — t i+i (14) 



Qi = n t + n l+1 (is) 
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FIG. 9 Polygonal disclination as in Fig. 8: splitting of Li 

if the rotation angles are small. In these expressions, the signs are such that the line Lj is oriented inward (toward 
Oj), whereas Lj+i and Lj are oriented outward. If the orientations are so chosen that they all are outward or all 
inward, one gets 

^2 sin^ t P = 0. (16) 
p 

This equation is valid for any number of disclination segments meeting at the same point O (see below). 

Equation 16 does not contain any reference to the directions of the disclinations segments meeting in O; it is akin 
to a Kirchhoff relation for the vectors 2 sin tp. Notice in particular that Eq. 15, which is valid for small angles, 

obtains straightforwardly by considering Frank circuits about the lines Lj, Lj+i and Lj. We shall christen the vector 
f = 2 sin ^ t that plays for a disclination line the same role as the Burgers vector b plays for a dislocation line, the 
Frank vector. It is oriented in the same direction as the rotation vector fit. 

By insuring that line Lj has no attached relaxed dislocations, condition 12 expresses the fact that L has a wedge 
character. Such a choice of L assumes that the twist lines have a larger line tension than the wedge ones because of 
their core energy, as discussed in II. E. 



b. Lines meeting at a node. 

We now deepen the disclination nature of Lj in the case considered above when three disclinations meet in Oj. As 

n» - 

defined above, it is obviously a wedge line — the Frank vector f \ = 2 sin — tj is along the line — that is split into 

~ t f2j ~ ~ t i^j_j_^ ~ 

two mixed lines, with fj = 2 sin — tj Frank vectors and fj + i = 2 sin — - — tj+i, Fig. 9. Again, there are no relaxation 

dislocations along Lj. It is therefore a very special wedge line. 

Now, line fluctuations <5Lj would break the continuity of the Burgers vector on Lj , and generate relaxation disloca- 
tions with Burgers vector 

Sb = 2 sin tj X <5Lj — 2 sin tj+i x <5Lj , 

according to Eq. 1 (8b — 8bi — <5bj + i), by adding the effects of the two lines involved in the splitting), i.e., from 
Eq. 11, 

^ ~ ■ — ^ ^ _i_ 1 ~* ' — ' — 

8b = 2 sin — tj X Shi — 2 sin — — tj+i x <5Lj, = Tj x <5Lj (17) 

from which equation it emerges that Tj is the Frank vector of the Lj disclination. These considerations also confirm 
that Eq. 16 is a Kirchhoff relation at a node where three disclinations meet. 

Notice that our reasoning has given a specific role to one of the three disclinations, L, but this restriction can be 
easily removed. Consider three Frank vectors 0i, 02, 03, and construct three disclination segments meeting at a 
common point O, such that: 



fl = 02 + 03, f2 = -03 + 01: f3 = —01-02- 



(18) 
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FIG. 10 Three broken disclinations of Frank vectors </>i, <j>2, 4>3, meeting at a point and composing three disclinations obeying 
Kirchhoff relation fi + f2 + f3 = (insert) 

This can easily be done by giving a sharp corner to the three disclinations (pi at O and joining their straight segments 
two by two, as in Fig. 10. If all the segments Lj are along f^, the L^s are all of wedge character, and there are 
no relaxation dislocations attached to them. The appearance of relaxation dislocations attached to the L,s would 
make them change direction. The final geometry depends on the energy balance between the grain boundaries (the 
constitutive dislocations), the relaxation dislocations, and the core energy of the disclinations. Notice that the three 
LjS are coplanar if they are of wedge character. 

The extension to any number of disclinations is obvious, which justifies our claim concerning Eq. 16. Notice however 
that, whereas Fig. 10 represents three grain boundaries merging two by two along three disclinations, the case of four 
disclinations (say) meeting at a node requires, in the most general case, four grain boundaries merging three by three 
along the four disclinations; each of them is then split into three subdisclinations. Such a geometry occurs by nature 
in ideal polynanocrystals, see III.C.2. 

Remark: In accordance with the remark at the end of II. A, Eq. 16 does not apply properly if one of the angles 

\Q p \ > 7T. 

3. Disclinations merging along a line 

The situation where three lines Li, L2, L3 merge along a unique line L is also worth considering. One expects that 
the Kirchhoff relation 

fi+f 2 + f 3 = (19) 

is satisfied. This case is physically represented in a ferromagnet by three Bloch walls merging along a Bloch line, 
the dislocations being the sources of the magnetoelastic stresses (Klcman, 1974). More generally, one expects that, 
in an amorphous medium, n disclinations • • • f j • • • merging along a line yield a unique disclination of Frank vector 

G. Generic disclination lines. Disclination densities 

We now come to the generic case when a disclination line L is smoothly curved and its Frank vector varies smoothly. 
Consider two infinitesimally close points P and Q on L, with Frank vectors fp and fq; we write, with obvious notations, 

Q = P + s 8s, fq — fp = 2 sin — — tq — 2 sin — - tp, Sq — Sp = ~Ss = —8s. (20) 

The variation between P and Q of the displacement on the cut surface of L in M, can be written: 

<5b = f Q x QM - fp x PM = -f P x s Ss + (f Q - f P ) x PM (21) 
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The first term (— fp x s Ss) measures the relaxation dislocation densities attached to the line between P and Q. We 
specialize to the second one, which measures the relaxation disclination densities. According to Kirchhoff's relation, 

we have —5f = fcj — fp, the sign being chosen so that the attached disclination densities — are oriented outward (as 

fcj), and fp inward. 

If we assume that L is a wedge disclination, i.e. tp = sp, tQ = sq, we have 

St = 5fl P t P + - sin — n Ss (22) 
R 2 

The first term of the right member (S£lp tp) measures the effect of the variation in modulus of the rotation vector. In 
its absence, the rotation vector of the attached disclinations is along the principal normal. This is obviously reminiscent 
of the bisecting disclination line. If the attached disclinations have a wedge character, i.e. are along the principal 
normal, again as above, there are no supplementary dislocations accompanying them, apart those constituting L, and 
we might expect that the energy is minimized. Another result, not visible in the previous analysis (II. F. 2), is that 
the curvature of the disclination line L is directly related to the presence of the attached disclinations. 

In the generic case (tp ^ sp, tQ ^ sq, S£lp ^ 0), the Frank vectors of the attached disclinations are no longer 
along n, but, as shown in the previous analysis, there is still a choice for the direction of the attached disclinations 
for which the supplementary dislocations are cancelled. And the conclusion on the relation between curvature and 
attached disclinations is still valid. 



III. COARSE-GRAINED CRYSTALLINE SOLIDS, GRAIN BOUNDARIES, POLYNANOCRYSTALS 
A. Coarse-grained crystalline solids 

Non-quantized wedge disclinations in a crystalline solid have been mentioned previously; they are akin to limited 
tilt boundaries discussed in II. 

Grain boundaries of small misorientation angle (subgrain boundaries) are well documented; they are limited by 
wedge, twist or mixed disclinations, according to the geometrical interactions between Burgers vectors and Frank 
vectors discussed previously. These interactions are restricted to those that imply Burgers vectors equal to translation 
symmetries of the medium. 

Large misorientation angle grain boundaries are especially important in polynanocrystals. 

A quantized disclination in a crystalline solid carries an angle f2 of rotational symmetry of the crystal; in the wedge 
case, the line is itself the axis of this symmetry, see Fig. 11. As already emphasized, the related energy is extremely 
large and thereby their existence improbable. A related imperfect disclination occurs when the grain boundary is 
a plane of geometry of large atomic density for the two grains. If this imperfect disclination is rejected outside the 
sample, one has a low energy twin. 

Quantized disclinations in mesomorphic media are discussed in IV and V. 
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FIG. 12 Classification of continuous disclinations in a solid: (a) wedge disclination; (b) normal tilt disclination; (c) pure twist 
disclination 



B. Grain Boundaries 

1. Classification of grain boundaries and of continuous disclinations 

Grain boundaries in solids are classified according to the orientation of the rotation vector ui with respect to the 
plane of the boundary Sgb : tilt grain boundary when uj is in EgBj twist grain boundary when uj is perpendicular to 
Y,qb- This classification makes sense: a tilt grain boundary can be split into a set of parallel identical edge dislocations 
whose Burgers vectors are perpendicular to the boundary; a twist grain boundary can be split into two sets of parallel 
identical screw dislocations whose Burgers vectors belong to the boundary. Such splittings are currently observed in 
small misorientation grain boundaries (also called sub-boundaries). Each set of screw lines of a twist grain boundary 
carries non- vanishing stresses, but the two sets cancel their stresses at long distance 3 . 

We have classified disclination lines according to the orientation of the rotation vector uj with respect to the line 
direction L: wedge line when uj is along L, twist line when uj is perpendicular to L. This classification is perfectly 
adequate when no account is taken of the presence of a grain boundary attached to the line, e.g. quantized disclinations 
(no grain boundaries), but is not consistent with the grain boundary classification. For instance a tilt grain boundary 
can be limited either by a wedge disclination or a twist disclination. 

A finer classification of continuous disclinations (i.e. carrying a grain boundary) seems therefore appropriate: 

(i) wedge disclination line, Fig. 12(a): L is parallel to the constitutive dislocations of a tilt boundary; uj is parallel 
to L. 

(ii) normal tilt disclination line, Fig. 12(b): L is perpendicular to the constitutive dislocations of a tilt boundary; 
uj is perpendicular to L. 

(iii) pure twist disclination line, Fig. 12(c): L belongs to a twist boundary; uj is perpendicular to the boundary, 



3 It is believed that there is only one set of screw lines in the twist grain boundaries of the TGBA liquid crystalline phase, see Renn 
and Lubcnsky (1988), Kamien and Lubensky (1999); in that case the long distance cancellation of the stresses carried by each grain 
boundary is due to the presence of a family of parallel grain boundaries at periodic distances. 
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FIG. 13 The three parallel disclinations D^, with rotations Wj, along an edge E between three grain boundaries £,i 
thus to L. We have taken |b>i| = |b>a | , hence b = bi + b2 perpendicular to L, as required. 

2. Polycrystals as compact assembly of polyhedral crystals 

Nearly perfect polycrystals, as possibly created by annealing, can be viewed as compact assemblies of polyhedral 
grains. Their common facets are commonly triangular and fairly flat grain boundaries E, each surrounded by a 
disclination line L of strength equal to the rotation uj of the grain boundary. 

Three grains meet along a fairly straight edge E bordering such a facet, where the three parallel disclinations 
combine along the edge, Fig. 13. In such stress free annealed polycrystal, each such triplet of parallel disclinations 
must compensate their long range stresses. 

a. Kirchhoff relations. To analyze the stresses due to the three disclinations of an edge E, we have again to 
distinguish the contributions of the wedge components (parallel to E, cf Fig. 12(a)) to those of the normal tilt 
and pure twist components u^x, with attached dislocations, Fig. 12(b) and (c). 

For this second part u^x, the compensation of the three families of relaxation dislocations lead to the same condition 
as above for a node: 



It is clear from Fig. 14 that the Volterra process which produces u^n and ui2,\\ by moving M' to M" and M" to M'" 
sum up to an effect opposite to that which u>3 n produces by moving M'" to M'. Thus, using Eq. 1 and the fact that 
M', M" and M'" are on a circle centred at E, one gets 




(23) 



with f (oj) = 2 sin — t. 

v ' 2 

To have a completely stress free edge, the wedge components Wj ii must also compose to zero: 




(24) 
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b. Subboundaries. For small misorientation boundaries, and also for boundaries whose orientation does not differ 
much from a small energy twin, all or part of the continuous dislocations cluster into a periodic distribution of 
quantized dislocations (Burgers, 1939; Read and Shocklcy, 1950). This polygonization was first observed by X rays 



22 




*3 

FIG. 14 Composition of the three edge disclinations of parallel rotations Wj ii 

and described in these physical terms by Crussard (1944), after annealing of FCC single crystals strained in multislips 
(stages II and III). Later observations after etching of low angle boundaries (Lacombe and Bcaujard, 1948) provided 
the first experimental proof of the decomposition of these boundaries into rows of dislocations, and various techniques 
such as electron microscopy for metals and semiconductors and pinning of dislocations by precipitates in transparent 
ionic solids analyzed the details of the dislocation networks on the subboundaries and the way these dislocations 
connect at the edges of the grains, as assumed above, cf. J. Friedel (1964, 1987). These dislocations can slip under 
stress, especially after annealing of crystals strained in single slip (stage I), which produces especially simple networks 
such as pictured Fig. 14 (Washburn and Parker, 1952); in the more general case, the bowing under stress of the 
dislocations of the various subboundaries decreases by a large fraction the effective clastic moduli (Friedel et ai, 1955). 

c. Large misorientation boundaries. Following G. Friedel (1926) who calls them 'macles par meriedrie'. Bollmann 
(1970) has established similar dislocation arrangements for large misorientation boundaries, in terms of a crystal- 
lographic network common to both grains in contact along the boundary. A common crystallographic network has 
also been put forward by G. Friedel (1926) for what he calls 'macles gcncralcs'. But generally it is believed that 
the boundary is an amorphous contact on an atomic thickness, with possible ledges along which one of the grains 
can overlap into the other. These configurations, as a whole, obey the same conditions of stability as small angle 
boundaries, but allow more stressed states than the former ones, like roughness, glide, lateral motions of the grains, • • • 

d. Specific complications that arise from the crystal structure. The coalescence of continuous distributions of in- 
finitesimal dislocations, as mostly considered in this paper, into quantized crystal dislocations can introduce some 
complications that should be stressed, as they have no equivalent in liquid crystals or magnetic structures. Some of 
the following are presented in Friedel (1964). 

(i) In the simplest cases, the infinitesimal dislocations coalesce into quantized dislocations all with the same direc- 
tions of line and Burgers vector. This condition optimizes the energy of contact between the grains at the expense 
of an elastic distortion of the crystal, over a distance from the grain boundary of the order of the distance I = ON 
between dislocations, Fig. 6. This is the case for three grain boundaries meeting along wedge disclinations, as pictured 
Fig. 12a and 14, when one set of edge dislocations rearrange after straining in stage I of a single slip system; another 
example is given by two (or three) systems of screw dislocations building a network of increasing density into a twist 
boundary, such as in Fig. 12c: this can be obtained by torsion of an hexagonal lattice along the hexagonal axis of 
symmetry, e.g. in graphite (de Gennes and Friedel, 2007) and in HCP metals (M. Fivel, private communication). 

(ii) In most cases, however, where dislocations of a number of slip systems have been developed by straining, the 
subboundaries af the polygonal structures obtained by recovery are composed each of a distribution of two or more 
systems of dislocations, so as to produce subboundaries of mixed natures and more or less random orientations. 
Such subboundaries are somewhat less stable, as their elastic distortions average out at a larger distance from the 
subboundaries for a given rotation, owing to the mixing of dislocations of several slip systems. 

(iii) Even in the simpler case first considered, Fig. 6 shows that the possible positions of the crystal dislocations can 
only occur at specific positions such as O and N along the subboundary. The periodic distribution of such dislocations 
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FIG. 15 Ledge in a merihedric twin, cf Friedel (1926), Bollmann (1970). The superlattice is underlined by black dots 

must be coherent with the crystal structure along the boundary; it must correspond to a discrete series of angles Wi, 
with intervals increasing with the angle at least for small angles u>i. A subgrain boundary with an angle between w; 
and u>i + i can only be built of successive segments corresponding to and Wj+i, and the distance of elastic relaxation 
away from the boundary is of the order of two such successive domains. 

(iv) Fig. 6 shows schematically the structure of such a symmetric tilt boundary, in its initial state with a continuous 
distribution of infinitesimal dislocations. Their regrouping into quantized crystal dislocations can only occur at 
distances Zj, where the Zj are integer multiples of the lattice period along the boundary (n = 2 in Fig. 6); Zj is related 
to oji by 

fj = — = sinwj. (25) 

H 

For increasing u>j, this special 'Frank vector' is increasingly smaller than Eq. 10 for the general case. 

(v) These WjS correspond to especially coherent structures, of lower energy. In fact, in crystal structures with a 
center of symmetry, they are twins; and Eq. 25 leads then for uji — n to a perfect coherence of the two crystals, 
with vanishing boundary tension. In the absence of a center of symmetry, un = 7r leads to a twin with especially low 
boundary tension and again no crystal dislocations. In both cases, the increase of U when uii decreases from 7r can be 
described in terms of an increasing density of crystal dislocations. More generally, perfect matching without crystal 
dislocations and vanishing grain boundary tension occurs when uji is an allowed rotational symmetry of the crystal, 
such as u>i = 7r/2 for the case of Fig. 6 with a cubic crystal symmetry. The same occurs of course in all cases for 
u>i = ±2n7r. 

(vi) 

In cubic structures, such twins are of the kind considered by G. Friedel (1926) and Bollmann (1970), in that they 
can be built on a common superlattice of the two crystals. Such twins, with not too large ks so uji is not too far from 
7r/2, can present ledges as sketched Fig. 15, which introduce supplementary dislocations such as D'2. 

C. Polynanocrystals 

Polynanocrystalline materials are assemblies of small polyhedral grains. The grain size ranges from a few nm 
to 1/xm. Their plastic properties, especially investigated in pure metals, pure semiconductors, and ceramics, show 
several features. One example is the presence of huge internal stresses, that distinguish them from the classic picture of 
dislocation driven phenomena in usual coarse-grained crystalline materials — see e.g., Weertman et al. (1999); Kumar 
et al. (2003); Ovid'ko (2005); Wolf et al. (2005); Van Swygenhoven and Weertman (2006). Polynanocrystals provide a 
remarkable physical example to discuss some of the notions just introduced. We are interested in the relation between 
the polynanocrystal plastic properties and the disclination and grain boundary structure. 
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1. Data on the plastic deformation of polynanocrystalline materials 

A dominant mechanism of coarse-grained crystal plastic deformation (work-hardening) at low temperature is the 
slip of dislocation pilc-ups; this mechanism obeys the well-known Hall— Pctch relation ay = oy a + kl^ 1 / 2 — where 
ay is the single crystal yield stress (the friction contribution), k a material dependent constant and I the grain size, 
see e.g. J. Friedel (1959a). 

But the yield stress does not increase without limit when I approaches atomic sizes. After reaching considerable 
values in the nanoscopic range, it decreases somewhat below some cross-over size l c —often referred as the strongest 
size (Yip, 1998); the material becomes ductile and even 'superplastic'. This latter property is not met in all samples, 
and depends crucially on the absence of porosity (high density) and of nanocracks. It has been observed e.g. in an 
iron alloy, (Branagan et ai, 2003), in nickel, (MacFadden et at, 1999; Schuh et ai, 2002), in copper, (Champion 
et ai, 2003; Koch, 2003; Lu et ai, 2000; Wang et ai, 2002; Zhu and Liao, 2004). Nanocrystallinc nickel, for instance, 
exhibits a Hall-Petch strengthening as the grain size decreases down to l c « 14 nm, thus reaching internal stresses of 
order at least ten times those observed in usual coarse-grained samples. 

High ductility requires suppression of plastic flow localization, i.e., strain hardening stabilizing the tensile defor- 
mation, see Ovid'ko (2005). Several mechanisms at the origin of this strain hardening are currently discussed in the 
literature. 

(i) Partial dislocations emitted by grain boundaries. The size range just below l c is characterized in FCC metals 
by the appearance of Shockley partials with 1/6 < 112 > Burgers vectors bordering {111} stacking faults and the 
(correlated) formation of twin lamellae, during which process the yield stress is still increasing; see e.g., for AI, 
Chen et al. (2003). In effect, one can imagine that the formation of partials on some glide system impedes the 
motion of partials on another one. A large number of atomistic calculations, critically reviewed and developed in Van 
Swygenhoven et al. (2006), supports a few experimental results pointing in this direction. 

We shall keep in mind that this type of dislocation-biased plastic deformation must put into play glides and 
probably lateral displacements or growth/shrinking processes of grain boundaries; the published data do not clearly 
assert whether these glides are analogous to low temperature glide or employ diffusion mechanisms, at the emis- 
sion/absorption of dislocations. 

(ii) Hardening by annealing and softening by deformation. Huang et al. (2006) have recently reported on the 
necessity of extremely high stresses in order to nucleate partials in well annealed, equilibrated, ultrafine nc grains of 
Cu with no intragranular isolated dislocations nor twins; on the other hand, the appearance of these grain boundary 
or/and disclination (probably) nucleated dislocations softens the material. This is reminiscent of the plastic behavior 
of whiskers (Ma et ai, 2006), and suggests that grain boundary perfection is an important factor. Indeed, the computer 
simulations alluded to above employ grain boundaries that are not totally relaxed and show up ledges that act as 
sources of partials. 

(iii) grain boundary sliding and grain rotation. It seems that this is the dominant mechanism at high temperatures 
and/or grain size below l c , with Coble (1963) diffusion inside the grain boundaries. See Schi0tz et al. (1998); Ovid'ko 
(2002); Van Swygenhoven (2002); Shan et al. (2004); Ma (2004). 

We distinguish in the sequel an ideal polynanocrystal configuration, with continuous disclinations and their consti- 
tutive dislocations, from actual configurations, with their defects and unusual large internal stresses. 

2. Structure of the ideal polynanocrystal 

The picture developed above for annealed coarse-grained poly crystals should apply to polynanocrystals. Thus the 
computer simulations of Van Swygenhoven et al. (2000) indicate that there is no difference between boundaries in 
polynanocrystals and in coarse-grained materials, and that the degree of organization is often rather high. 

Two main differences arise from the conditions of preparation: 

- sintering can produce large internal stresses or intergranular cavities. These large stresses should displace possible 
low angle boundaries, which are fairly common in coarse-grained polycrystals (Friedel et ai, 1953): such subboundaries 
should be mobile enough to glide towards large angle boundaries with which they should be integrated. 

- the usual size of the mosaic structure of crystals, whether in the form of a simple Frank network or of a finely 
polygonized structure, should not be present in polynanocrystals near their maximum elastic limit, both because of 
their destruction by the high internal stresses just referred to, but more generally because annealing conditions should 
destroy mosaic structures of sizes less than typically 10 microns (J. Friedel (1964) p. 240). With these provisoes, the 
conditions of stability of a polynanocrystal should follow the same Kirchhoff relations developed above for coarse- 
grained materials. 

A*b 

The surface tension E of a subboundary depends on the misoricntation uj; E Eo^, where Eo = — -, r; on 

47r(l — v) 
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FIG. 16 A modification process frequent in foams and possible in polynanocrystalline media 

the other hand, for large misorientations, above some value to > uj max ss 20 degrees, E can be taken independent of 
the direction of the grain boundary with respect to the two grain orientations (E w Eo ui max ), — see J. Fricdcl et al. 
(1953), J. Fricdel (1964), chapt. 10 — as long as the grain boundary is not along or close to a lattice direction common 
to the two grains. 

Most of the grain boundaries in polynanocrystalline materials arc large angle ones; as a consequence they have 
an energy fairly independent of the angle of misorientation. If submitted to this sole constant surface tension, the 
grain boundaries should form angles of ~ 120 degrees at triple junctions. Other forces are originating from the wedge 
disclination segments. The equilibrium can be reached, in principle, by processes of the form described Fig. 16, 
proposed by J. Fricdcl (1987) for solids and well known in foams (Weaire and Hutzlcr, 2000), which consists, first, 
in the vanishing of the boundary segment AB, and subsequently, after a passage through an unstable quadruple 
junction, the appearance of the boundary segment A'B'. In this figure A, B, A', B' are triple junctions seen end on, 

and u> a , uj a 2 , u>g\ u)jj\ ^A' ' w > w ' are ^ ne we dge component misorientations of the GBs whose sections with the 
plane are drawn; these rotation vectors are all directed along the normal to the figure. 

The triple junction Kirchhoff relations can be written: 
along A: to + ur]/ + co^ — Suj a, 

along B: — u> + Wg' + uig = 5u>b, 
along A': w' + u)^ + w s = Suja>, 

(2) (2) 

along B': — w + uj\' + uj B ' = 5ub>- 
They yield 

5u>a + Slub — 5^>a' + Su>b' (= 2Soj). 
Hence we can write 

Slua = Suj + 5m, Sujb — Suj — 8vj, 

5u>a' = Suj + dm', 6lob> = Slu — 8vj' , 
where 

25w — Slu a + Su>b, 25w' = Sua' + Slob 1 - 
Suj measures the repulsive terms between A and B, A' and B'; they are equal. The attractive terms Szu, 8m' are 
different; thus the modification is favored if \Sm'\ < \8m\. In the case when Suj = 8m = 8m' = 0, the modification is 
ruled by the surface tension. 

There is probably no chance, although it is in principle possible, that a polynanocrystal reaches the same structural 
configuration as foreseen for a foam, and even less that this structure coarsens, as 2D foams do (von Neumann, 1952), 
a result extended recently to 3D foams (Hilgcnfcldt et al., 2001). 
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FIG. 17 Slip of a pile up, left half-plane, against a normal tilt disclination 
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FIG. 18 Double kink on a disclination. The dislocation Burgers vector is in the plane of the grain boundary or close to it; see 
text. 

3. Plastic deformation of a polynanocrystal 

The lack of a mosaic structure in fine polynanocrystal prevents the presence of internal Frank Read sources in the 
grains. This is well recognized by most authors. We shall consider two possible processes: production of partials and 
grain boundary sliding, grain rotation and more generally grain boundaries as sources and sinks of dislocations. 

a. Production of partials. There probably are several mechanisms possible for the production of partials, some relating 
to the grain boundaries, other ones to the disclination lines (triple junctions). 

(i) Constitutive edge dislocations: the role of the applied stress in a possible bowing of the dislocations of the 
boundary has already been put forward. This is probably a low temperature mechanism, by glide, but a number of 
computer simulations (see Van Swygcnhovcn et al. (2006)) as well as experimental works indicate that the process is 
thermally activated (Van Pctegcm et al., 2006), and displays a small activation volume . The same dislocations are 
free, given suitable stresses, to climb in the plane of the boundary and thus to displace the bordering disclination, 
thereby displacing the two other boundaries that merge along the common edge along which the disclination is located. 

(ii) Slip of dislocation in the boundary, whose Burgers vectors are in the plane, Fig. 17. 

Such dislocations can pile up along the edges and be at the origin of stress concentrations large enough to nucleate 
new dislocations and/or move the other grains, as proposed by Ovid'ko (2005). Notice that the dislocations considered 
here are not constitutive or relaxation dislocations of the boundary to which they belong (they do not obey Eq. 1). 
In fact they form a classical piling up. 

(iii) grain boundaries and triple junctions as sources and sinks for dislocations. The emitted/absorbed dislocations 
are not constitutive dislocations of the grain boundaries, neither pairs of dislocations of opposite signs (dipoles). 
Assume the presence of a ledge on a grain boundary bordered by a triple junction. If this ledge affects only one 
(possibly two) of the three grain boundaries merging along the triple junction, there is necessarily some kind of 
splitting of the triple junction into its constitutive disclinations, in the region of the ledge. Consider one of these 
'stripped' disclination segments; the ledge determines on it a double kink AB, A'B', Fig. 18. According to Eq. 2, the 
Burgers vector of the dislocation which necessarily joins the two kinks is perpendicular to the kinks and to the rotation 
vector ui of the kinked disclination; uj is certainly close in direction, if not parallel, to the disclination line, according 
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to the analysis of III.C.2. Thus a large component of the Burgers vector of the dislocation is in the plane of the grain 
boundary; this might favor slip in this plane. Also, the smaller the kink lengths AB, A'B', Fig. 18, the smaller the 
Burgers vector of the dislocation attached to them, the smaller the energy to nucleate the double kink. One can 
therefore speculate that this is the origin of partials (rather than perfect dislocations) and of the related small ac- 
tivation energy and volume (by two orders of magnitude as compared to the values observed in coarse-grained metals) . 

b. Grain Boundary Sliding and Grain Rotation. These mechanisms of deformation have been mentioned, III.C.l. 
Notice that a relative rotation Slu of two grains with a common grain boundary requires that the constitutive dislocation 
densities be modified accordingly, i.e. that an exchange of grain boundary dislocations with the intergranular medium 
takes place. Huang et al. (2006) have observed that the number of subboundaries decreases under high T annealing, 
so that the grains necessarily 'roll' (Shan et al., 2004), the trend being seemingly toward a polynanocrystal with large 
misorientation grain boundaries only. 

IV. QUANTIZED DISCLINATIONS IN MESOMORPHIC PHASES 

There is no conceptual difficulty to construct in the manner of the Volterra process a quantized wedge disclination 
in a medium endowed with finite rotational symmetries; e.g., a nematic phase £1 = ±rr, a smectic phase (f2 = ±7r), 
liquid crystals in general, a solid crystal (ft = i 2 ^-, n= 1, 2, 3, 4, or 6), a quasicrystal (Bohsung and Trebin, 1987) 
(n — 5,8,10,12, etc) or a 3D spherical or hyperbolic curved crystal (Kleman, 1989). In 3D solid crystals, the line 
energies are so large that disclinations are observed in special conditions only. See II. A (3D solid crystal wedge lines) 
and II. E (continuous twist lines) for comments. 

The disclinations that are observed in various liquid crystals usually differ widely from what is expected to result 
from a pure (not extended) Volterra process. These differences originate in various liquid crystal symmetries, thereby 
various types of relaxation defects. For, as in the previous section, one can consider the interplay of the quantized 
wedge disclination, either with attached dislocations, which transform it into a twist or mixed disclination with the 
same rotation vector fi, or with attached disclinations, which yield a f2 variable in direction that illustrates the large 
rotation deformations that a liquid crystal can suffer, or again with unattached dislocations that result from their 
motion. Quantized disclinations constructed by such extended Volterra processes can be described in terms of twist, 
wedge, or mixed segments, in addition to their physical property of topological stability. 

Two questions therefore arise; a)- how do the Volterra characteristics of wedge, twist, or mixed character (i.e. 
different types of extended Volterra processes) reflect in the topological classification, b)-whether any quantized 
disclination, empirically given, can be constructed in a systematic way by an extended Volterra process. To point a) 
we have a partial answer, namely that disclinations, when differing only by constitutive dislocations, belong to the 
same conjugacy class of IIi(V). Point b) is reported to the final discussion (VIII). 

A. Quantized wedge disclinations and their transformations 

This subsection is devoted to the molecular configurations carried by quantized disclinations, when the limitations 
that are imposed in the Volterra process by the specific symmetries of the medium are taken into account. Examples 
given in this subsection relate to nematics (N) and cholesterics (N*), and in the next subsection to SmAs. It appears 
that important disclination properties (shape, flexibility, interplay between them and with other defects, etc.) escape 
an analysis based on the sole topological classification. 

1. N phase 

The order parameter space of the nematic phase is the projective plane P 2 , whose first homotopy group IIi(P 2 ) is 
Z-2, the group with two elements {e, a}, a 2 = e, e being the identity. All topologically classified defects belong to a 
unique class of homotopy, namely a; all Volterra defects of strength \k\ = n + |, |f2| = (2n + 1)tt, n G Z + 1J{0} can 
be mapped on a. Most experimental observations yield n = 0, i.e. two different Volterra disclination types, k = +| 
and k — — \, which the topological theory classifies under the same heading; indeed, the topological theory predicts 
that it is possible to transform smoothly a k — +| into a k — —\ (Bouligand, 1981). In fact, since the Volterra 
process predicts also k = n (these lines, again, are not topologically stable), the cases k = +\ and k = — ^ differ by 
the not topologically stable but Voterra line k = 1. The Volterra classification thereby proves useful when analyzing 
experimental results. 
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Consider a wedge line. Any axis orthogonal to the director is a twofold symmetry axis. The wedge line has to be 
along such an axis. As a consequence, there is a director that is orthogonal to the line in its close vicinity. 

The deformation of a wedge straight line implies a translation and/or a rotation of along the line. Let us apply 
the considerations of II. G, Eq. 21 and 22. The translation and rotation of f2 along the line bring 

(i) a non-vanishing density of attached dislocations 

Shxr — — 2tp x sp Ss (26) 

that vanishes if the Frank vector is along the tangent sp in P, i.e. if the deformed disclination is still of wedge 
character, 

(ii) a non-vanishing density of attached disclinations 

Sf P = |n P Ss, (27) 

(this expression being valid in the pure wedge case only) whose infinitesimal Frank vectors <5fp are along the principal 
normal np to the line in P. Such a defect configuration is licit if the direction of np is along an actual director, because 
any director is an axis of continuous rotation symmetry for the N phase. This is in agreement with our remark above, 
according to which there is a director that is orthogonal to a wedge line in its close vicinity. The density of attached 
disclinations is vanishing if suffers a pure translation along the disclination line. 

Continuous defects belong to the identity element of the first homotopy group; this is why their distribution, which 
is a function of the shape of the line and of the field, does not influence the topological invariant, always a for any 
|fc| = ^ disclination line. 

A | A: | — 1 disclination line does not require special configuration rules for the director in the near vicinity of the 
line, because any axis is a 2tt rotational symmetry axis. It is then possible to align the director along the line, at the 
expense of special densities of defects (Klcman, 1973). 

Remark: The above description of the curvature of a disclination line in a N phase in terms of defect densities 
might look overdone; but, at least in what concerns dislocation densities, it is no more so than the description of 
the strains and stresses in an amorphous medium in terms of defects. There is no difference of nature between the 
non-vanishing density of dislocations in a nematic and that one of an amorphous medium, except that in the first 
case they originate on disclinations and that the rules of elastic relaxation differ. In both cases the 3D continuous 
translational symmetries render the dislocation densities trivial (this is not so in most liquid crystalline phases, the 
biaxial nematicNs being a exception, see below) . The disclination densities (continuous rotational symmetries about 
the directors) are superimposed on the dislocation densities but are independent; it would be interesting, knowing 
the field of distortions of a N phase, to separate what is due to dislocation densities from what is due to disclination 
densities. 

2. N* phase 

The interplay between continuous dislocations and quantized disclinations in cholcsterics has been discussed by the 
present authors (Friedel and Kleman, 1969) in the very paper where this concept was first introduced. There are 
3 types of quantized rotational symmetries ft in a N* phase, all multiple of the angle it: 1)- along the molecular 
axis, which we denote A; 2)- along the helicity axis, x; 3)- along the transverse axis, r. See Fig. 19. The related 
Volterra processes yield similar results (similar limitations) to those above for the N phase, since A, % and r are 
directors; being a twofold axis, any direction orthogonal to A, e.g. x an d t, is necessarily singular on the core of 
a Volterra-constructcd related disclination, except if the rotation vector is an integer multiple of 4-7T, as shown by 
Anderson and Toulouse (1977). Hence Volterra processes are perfectly defined for k\ = 2n + y, n, m G Z (and by 
an obvious extension for k T = 2n + y , and k x = 2n + y ) . Topologically stable disclinations are classified by the 
elements of the non-Abelian quaternion group Q%, whose elements are usually denoted {±1, ±z, ±j, ±fc} in quaternion 
notation, or {±e, ±i<7i, ±«02, ±icr 3 } in a 2x2 matrix representation where the cr^'s are the Pauli matrices and e the 
unit 2x2 matrix. We employ the quaternion notation, which is more appropriate to crystals in curved spaces of 
constant positive curvature (Coxeter, 1991), see Section VI. {+1} corresponds to k\, k x , k T = 2n, {—1} to k\, k x , 
k T = n odd, {±i} , {±j} , {±k} to k\, k x , k T half integers, respectively (Klcman et at, 2004; Mcrmin, 1979). 

Employing the Volterra method, we now look for the possible attached continuous defects and their role in the 
flexibility of quantized disclination lines. 

a. Attached defects: continuous dislocations; k\, k x , and k T — ±| lines. 

The only possible continuous dislocation Burgers vectors are parallel to the cholesteric planes (orthogonal to the 
helicity axis x), which are invariant under any in-plane translation. Consider then a disclination line L, with tangent 
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FIG. 19 Cholesteric phase N*. The A and t directors indicated on the left side belong to the extreme left molecules. The x 
director is constant all over the figure. From Kleman and Lavrentovich, Soft Matter Physics, an Introduction, Springer, 2003, 
with permission 

vector t at some point P of L. The associated attached Burgers vectors are along the direction i~2 x t, according to 
Eq. 1. Therefore there is no topological obstruction to the flexibility of a line L in a plane perpendicular to t if L 
is a A disclination, and in a plane perpendicular to A if L is a r disclination, but no other types of flexibility are 
allowed for these lines. On the other hand, a \ line could curve in any plane; another way to state this latter result 
is to notice that a \ disclination line of strength k is also a dislocation of Burgers vector b = —kpx, because of the 
equivalence of a tt rotation along the x axis with a \p translation along the same axis (p is the pitch) , cf. Fricdel and 
Kleman (1969) and Bouligand and Kleman (1970). 

b. Attached defects: continuous dispirations; k x and k T — ±| twist lines. 

The other continuous symmetries in a N* phase are helical rotations {(5f2, — p^p-}, 5tt = 5&X along the x axis, 
i.e. the combination of a translation and of a rotation. The corresponding Voltcrra defect is a continuous dispiration 

that combines a dislocation and a disclination. Applying Eq. 1 to the translational part — p — , the tangent t„ to 
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the line at P has to be in the {A,r} plane (perpendicular to x)- The rotation vector direction ftp, which we write 
rip = il rop, |tnp| = 1, is along A or r. The vector — -c&q — np (it appears in Eq. 21, where it is denoted St) has 
to be along X- We have, after Eq. 21, and assuming that the disclination line strength is k = |: 

sn 



2 w X t Ss = — p (dislocation component) 

2 8zu — Sfl (disclination component, 



(28) 



and by elimination of 5£l, 
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This equation means that the rotation rate of Qp is — tp. Therefore the rotation rate of the Frenet trihedron 
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attached to the disclination line at P is 
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Because, according to Frenet's formulae, we have uj — r tp + p bp (with obvious notations), we infer that the 

P 

disclination line has a constant torsion (r = — ), the same for all lines of this type, and that z&p = bp is along the 

27T 

binomial of the disclination line. The rotation vector ftp is therefore along the binormal, which is either a A or a r 



30 



direction. Because of Eq. 28, and because 6£lp oc x , it follows that the local x ax is is along the principal normal np, 
that the rotation vector (A or t) is along the binormal, and that the tangent to the line is a t or A direction. The 
disclination is of pure twist character. The search for constant torsion curves starts with Darboux (1894) and has been 
the subject of recent investigations, relating them to the celebrated Backhand transformation and the classification of 
surfaces of constant negative curvature; see e.g. Calini and Ivey (1998). Among the solutions, the simplest ones are 

4-7T 2 / 1 47T 2 \ - 1 

those for which p(s) is a constant; the curve is then a circular helix with pitch P = ( — H — 1 and radius 



R = — I — =• H ) . The p = oo limit case is simple and interesting. The disclination line is straight and orthogonal 



to the x axis. Such an object is highly energetic, but can be stabilized by the presence of another straight disclination 
of opposite sign, parallel to the first one and at a short distance. Continuous dispirations link the two lines in the 
ribbon in between. Now, because the set of two disclination lines of opposite signs is equivalent to a dislocation, the 
ribbon can take any shape, which is comparable to the case discussed by J. Friedel and Kleman (1969). 

As reviewed by Calini and Ivey (1998), there is a considerable variety of closed constant torsion curves with different 
knot classes. The search for disclinations affecting such shapes in N* phases remains to be done. 

The foregoing analysis of the flexibility of disclinations does not leave place for curved lines in the (A, t) plane 
orthogonal to the x axis. But such lines are known to exist, e.g. ellipses belonging to focal conic domains (FCD), 
much akin to focal conic domains in SmA's (Bouligand, 1972, 1973, 1974). In the limit where the pitch is small 
compared to the size of the sample, and the size of the FCD so large compared to the pitch that the inner helicity 
of the N* layers is a negligible phenomenon, the analysis of these ellipses (and of their conjugate hyperbolae) can be 
conducted similarly to FCDs in SmAs (see below). But in most experimental cases the situation is more delicate and 
a thorough investigation is all wanting. 

Remark 1. The above discussion relates to disclinations of strength |fc| = |, of homotopy classes {±i},{±j}, or 
{±k}. The case of \k\ — 1 defects (homotopy class {—1}) requires a different approach, since Eq. 1 is no longer valid, 
and we have to treat the |fc| = 1 defect as a sum of two |fc| = \ defects. There is no objection in principle to uniting 
two disclinations of the same strength ^ along the same line, making then a disclination of unit strength. For instance, 
in the constant torsion case, such a process of adding two \k\ — g is a way of canceling the singularity on the core, 
if the rotation vector is along the r director (the A director is then along the line); however the singularity of the 
order parameter itself, which consists of the three director trihedron, is not canceled. One can infer that, in the set 
of defects just investigated (IV. A. 2), the twist = 1 line is favored. 

A situation where a non-singular defect of seemingly |fc| = 1 strength cannot be split into two |fc| = | defects is 
discussed by Bouligand et al. (1978), but in fact it relates to non-singular topological configurations, in the sense of 
Michel (1980), classified by the Hopf index for the director field, not line defects. 

Remark 2. The biaxial nematic phase (Nb ) has the same topological classification of disclinations as the N* phase 
(Toulouse, 1977a; Volovik and Mineev, 1977). But the symmetry group is different. Nb is invariant under any 
translation (i.e. belonging to E 3 ) but there are no continuous rotation symmetries. Hence, the only possibility to 
curve for a Nb line is by attaching continuous dislocations, in strong contrast with N. This is an example where the 
Volterra process appears to give a more detailed view of the defect conformation than the topological theory. 

B. SmA phase 

1. Wedge disclinations 

Figure 20 shows a k = — | wedge line in a Sm phase 4 , displaced from L to L' by the absorption of an edge dislocation 
|b| = do whose Burgers vector is twice as large as this displacement \d$. 

The nature of the core has changed. The absorption of a second dislocation equal to the first one, along the same 
route, would displace the line by the same amount, the total effect of the two displacements being equal to the repeat 
distance d Q of the layers, L being moved to a position L" (not drawn Fig. 20), where the original core is retrieved. 
The analogy with the displacement of a wedge continuous line described in II. B is striking; the configuration of L", 
compared to L, displays a full new layer equivalent to a dislocation of Burgers vector |b| = 2d$ for a displacement 
do of the disclination line, as obtained by applying Eq. 1. The inverse displacement requires the appearance of a 
dislocation line on the core, which relaxes and eventually disappears far away from the line 5 . 



4 The same picture is valid for a 3D crystal, see the 2D cut along a lattice plane Fig. 11. 

5 It is often taken for granted that dislocations always nucleate by pairs of opposite signs. But we have here a situation where a dislocation 
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FIG. 20 The disclination located at L (a, b)) is displaced to L' (c) by the absorption of a dislocation of Burgers vector unity 
(6=1, LL' = |). The dislocation visible in (b) has disappeared in the process (c) 
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FIG. 21 The k = | wedge SmA disclination; (a)- Nye's infinitesimal dislocations as agents of layer curvature; (b)- Calculating 
the sum total of the infinitesimal dislocations belonging to one layer 

F. C. Frank was the first to point out, in an oral communication at the 1969 Montpellicr Symposium on Liquid 
Crystals, that the displacement of a disclination in a liquid crystal (he used a cholesteric phase N* to illustrate his 
comment) involves the emission/absorption of dislocations, which are quantized in the case he considered. 

2. Nye's relaxation dislocations 

There is a feature not apparent in the analysis of constitutive and relaxation dislocations carried out in II, namely 
the possibility of continuous relaxation dislocations that are directly related to the curvature of the layers, and not 
attached to the line. Figure 21(a) is illustrative of the case k = |. Let us turn our attention to one of the layers 
inside the disclination wedge. As a liquid layer, its 'inner' group of symmetry is E 2 ; it thereby admits continuous 

dislocations whose Burgers vectors are parallel to the layer. Those dislocations determine the curvature — of the 

P 



line is nucleated with no partner of the opposite sign. We believe that this possibility is relevant in some important cases. For example, 
one might this way nucleate screw dislocations all of the same sign at the SmA — + TGBA transition. 



32 




\ 

11 



-±-- ± -J- | 



FIG. 22 The k = — | wedge SmA disclination; A, nA and B, iib are at infinity along the asymptotic directions of the disclination 
configuration. The edge dislocation density vanishes at infinity 



layers, the relation between the dislocation density and the curvature being — = — , where do is the layer thickness. 

ds p 

This relation was first established by Nye (1953) for solid crystals (the curvature of the lattice planes is a function of 
the dislocation densities), see Appendix A. 

Let t be a unit vector along the tangent to the layer in a section perpendicular to the wedge line, and let us traverse 
a path AB in this section, everywhere tangent to the layer along t, in the part of the layer which is curved, see 
Fig. 21(b). The total Burgers vector measured along the path from A to B ( A, B any points on the upper and lower 
horizontal parts of the path) is: 



i.e. bAB = 2e?o no, where no is the normal to the layers far from the disclination (we have employed one of Frenet's 

formulae to transform the second integral into the third, namely — — = , where r(s) is the torsion and b(s) the 

ds t p 

binormal at a point s on the path). 

Observe that bAB is of a sign opposite to the Burgers vector of the constitutive dislocations of the disclination 
line; observe further that the result does not depend on the precise shape of the layer, whose possible strain at finite 
distance does not invalidate the result of Eq. 31, as long as the layers are parallel and planar far away from the 
disclination. 

Equation 31 establishes that the set of infinitesimal dislocations attending a k = \ line relaxes the elastic stresses 
due to the constitutive dislocations, in the geometry we are considering. This is also true for a k = — | line; the 
application of the method of Eq. 31 to Fig. 22 shows that the equivalent Burgers vector is bAB = 2dono (= 2donA). 

Any layer curvature can be analyzed in terms of Nye's dislocations, see section V for the case of Focal Conic 
Domains in SmAs. 

The example that we have just developed evidences the main features of the relaxation processes relating to Nye's 
dislocations. 

(i) the relaxation of the elastic stresses carried by a disclination which results from a pure (non extended) Volterra 
process is obtained by the glide of the Sm layers past each other, which accumulates infinitesimal dislocations; 

(ii) the accumulation of these dislocations along the layers is equivalent to a set of infinitesimal disclinations attached 
to the master disclination, analogous to those discussed in II. G (see a similar discussion relating to focal conies, which 
are special types of disclinations, in V). It is also equivalent to a piling of subboundaries; 

(iii) the stresses resulting from the fact that the molecules inside the layers are compressed at one end and stretched 
at the other introduce an elastic constant at the origin of the Frank— Oseen splay constant for smectics; 

(iv) finally, the Nye's dislocations geometry screens the long distance stresses carried by the master disclination, 
which considerably reduces its line tension. 




(31) 
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3. Topological stability and Volterra process compared in SmAs. Twist disclinations 

How do the topological stability approach and the Volterra process approach compare in SmAs? The elements of the 
first homotopy group IIi(V) classify line defects (dislocations and disclinations). For a SmA, IIi(V SmA) ~ Z\Z\Z2- Let 
(n, a) denote an element of ZOZ2, with n € Z and a € Z 2 (= {e, a}). With these notations, the identity (null) defect 
is denoted (0,e), a dislocation is denoted (n, e), where n stands for the Burgers vector, and a |f2| = n disclination is 
denoted (0, a), irrespective of the sign of fi; a 2 = e. 

Clearly enough, (n, a) is a |f2| = tt disclination that has absorbed a dislocation (n, e), i.e. that has suffered a 
translation nd,Q. Therefore in Fig. 20, if L is represented by the homotopy class (0, a), (which we write L 1— > (0, a)), 
then L' 1— * (1, e)(0, a) = (1, a). Notice that the nature of the core has changed. The absorption of another dislocation 
of the same sign yields a disclination L" 1— > (2,e)(0,a) = (2, a) (not represented Fig. 20), with the same core as L. 
The product of two disclinations yield the identity (0,e). All these operations are summarized in the multiplication 
rules: 

(n, a)(m, 0) = (n + a(m), a/3), e(m) = m, a(m) = — m. (32) 

Observe that any disclination, whatever the nature of its core, can be chosen arbitrarily as the origin disclination 
(0,a). 

With the above analysis, L, L' and L" are given three different topological invariants. But L and L" can also be 
gathered under the same heading in the frame of the topological theory; they belong indeed to the same conjugacy 
class of IIi(V SmA) —see Appendix B— and this is enough to consider them as the same topologically stable defect, 
according to the general topological theory (Klcman, 1989; Michel, 1980). Consider the two equalities: 

(2, a) = (0,a)(2,e),(2,o) = (-1, e)(0, a)(l, e). (33) 

The first equality means that L" is obtained by 'adding' the dislocation 1— > (2, e) to L. In the Volterra sense, it 
is an addition; in the topological stability theory sense, it is the product of two homotopy classes (0,a) and (2,e). 
The second equality, which expresses that (0,a) and (2, a) are conjugate in LTi(V SmA), has a simple physical image: 
it expresses the effect of a complete circumnavigation of the disclination L 1— > (0, a) about a dislocation 1— * (1, e). 
Such an operation cannot change the nature of the circumnavigating defect, although its homotopy class is modified. 
We refer the reader to Mermin (1979) for a nice pedagogical discussion and an illustration of this property. 

In terms of Volterra invariants, L, L" and all disclinations of the same conjugacy class carry rotations about twofold 
axes located in planes between layers; L' and all disclinations of the same conjugacy class carry rotations about twofold 
axes located in the middle planes of layers. 

C. Nature of the defects attached to a quantized disclination 

Taking stock of the specific examples discussed above, we now derive some general properties relating to attached 
defects. Various cases arise: 



1. Continuous attached defects (dislocations, disclinations, dispirations) and kinks 

a. Topological stability. Continuous defect densities belong to the identity homotopy class, therefore they do not 
modify the homotopy class of the master disclination L all along it. Or, said otherwise, they are not visible when 
mapping a closed loop of the deformed medium into the order parameter space V. 

b. Volterra process. Because the existence of continuous defects has to comply, in an ordered medium, with the 
existence of broken symmetries, not all continuous defects are realizable, and thereby limits are put on the possible 
realizations of master disclinations, in particular their shapes (in dynamic terms, their flexibility). Equation 21, which 
derives from Eq. 1 and has been established for an isotropic uniform medium, is still valid, these limits being taken 
into account. 



2. Quantized attached defects of the first type; full kinks 

a. Topological stability. Inspired by the SmA example, we first refer to a case where L 1— > {a} is transformed into 
a disclination L" 1— > {a"} belonging to the same conjugacy class {a"} = {u}{a}{u -1 } of the first homotopy group 
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FIG. 23 L is a quantized disclination, AB a kink. The relaxation dislocations no longer cross the master line 

Ili(y) as L, by the absorption or the emission of a dislocation {v}. This dislocation hits the master disclination at 
some 'node', where we have a relation of the Kirchhoff's type, which is written in the language of the topological 
theory {a"} — {a}{v}; this relation also reads: 

{v} = {a" 1 }^"} = {a-^MiaUu- 1 }, (34) 

therefore the attached defect is a commutator of LTi(V). The result is not restricted to attached dislocations; it is 
valid for a defect {v} of any type. We now establish a reciprocity theorem, which extends the as yet not understood 
concept of node. 

Let us recall that the commutators of a group G generate an invariant subgroup D[G], also called the derived group. 
Not all elements of the derived group are commutators, but all are products of commutators. From the point of view 
of the physics of defects, it is important to observe that D contains entire conjugacy classes of II^V), and that the 
cosets of D in LTi(V) are also composed of entire conjugacy classes (Klcman, 1977; Trebin, 1984). We take an example. 
In a SmA, the dislocations having the same Burgers vector parity all belong to the same coset of Hi(V smA) let 
us, for the sake of clarity, consider only even dislocations: the two homotopy classes (2r, e) and (— 2r, e), r^O, form 
an entire class of conjugacy and they belong to the same coset as the homotopy class (0, e) - the null defect, which 
constitutes by itself a full conjugacy class. This very fact means that the even dislocations are equivalent to the null 
dislocation (0,e), in the following sense. It is possible to split the even dislocation (2r, e) into two equal dislocations 
(r, e); let one of them circumnavigate about a (m, a) disclination, which circumnavigation brings it to the conjugate 
state (— m, a)(r, e)(ra, a) = (— r, e), while the other one stays in place. It is then possible for (2r, e) to self-annihilate 
by letting the fixed and the returned mutually annihilate, i.e. generate a defect of homotopy class (0, e). Similarly, 
any odd dislocation (2r + 1, e) is equivalent to (1, e); but odd dislocations are not commutators. 

One can therefore establish the following reciprocity theorem: 
if the attached defect {v} is a commutator, then {a} and {a"} are either in the same conjugacy class, or belong to two 
conjugacy classes belonging to the same coset of IIi(V)/D. 

This is also true, by an easy extension, for any clement of D, not only commutators. We can thereby state in all 
generality that 

any defect whose homotopy class belongs to the derived group D[II 1 (V)] is eligible as an attached defect, and separates 
the master disclination into segments whose homotopy classes belong to the same coset of IIi(V)/D. 
We name such a defect an attached defect of the first type. 

b. Volterra process. Insofar as defects of the first type are commutators, they can terminate on a singular point, 
since they are equivalent (in the sense just defined) to the identity homotopy class, and this singular point can be the 
node where they meet the master disclination. This is the main result one can reach from the analysis of topological 
properties, but not more, because the topological theory does not separate properly dislocations and disclinations. In 
other words, the topological analysis does not say anything about the shape (the flexibility) of the master line, even 
though there is no doubt that the attachments are the tools for its changes of shape and the relaxation of the stresses. 

Let us consider the geometry of a kink on a wedge disclination (Fig. 23); this geometry differs somewhat from 
that of Fig. 4 for a continuous disclination. According to the above discussion, the same dislocations that have been 
absorbed by the wedge segment L~, say, are still outside the wedge segment L + ; we still have relaxation dislocations, 
terminating on the kink, an allowed process insofar as these dislocations are of the first type. 
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The Friedel's relation (Eq. 1) is established for an a priori Volterra description of line defects. Let us assume that, 
in Fig. 23, L + and L~ are k = ±| wedge disclination segments in a SmA, as in Fig. 20. According to Friedel's 
relation, which can be written, with f2 = ±7rt 

hi = ± 2 1 x AB (35) 

i 

the relaxation dislocation Burgers vectors are perpendicular to the figure plane. If we assume that we are in the 
conditions of application of this relation in a SmA, the lines are drawn in a medium that is not yet deformed by the 
Volterra process, the layers are parallel to the L + and L~ lines. We do not lose generality assuming further that the 
layers are perpendicular to the plane of the figure; AB is along the layer normal. Therefore the Burgers vectors b; 
are parallel to the layers; the related dislocations are continuous. The layers in the transition region between L + and 
L~ suffer extra curvatures, which represent these dislocations. This is certainly not a small energy geometry. But the 
geometry of Fig. 23 can be understood differently. If the relaxation dislocations are quantized, there should be layers 
parallel to the plane of the figure in the AB region (perpendicular to the Burgers vector) , whereas AB is along the 
normal to the layers. This is possible only if we are considering the medium already deformed by the disclinations. 
The Friedel's relation still works for AB joining a segment along L to a segment along L , for example, in fig. 20. 
It therefore works when applied locally to the 'tangent' undistorted medium, on either sides of the disclination. We 
give an example in a SmA phase section V. We name full kinks the kinks that separate master disclination segments 
belonging to the same coset of ITi(V)/D. 

3. Quantized attached defects of the second type; partial kinks 

The discussion of the SmA case has shown that it is perfectly licit to consider a master line made of two wedge 
segments L and L' that do not belong to the same coset of Ili(V)/D. The dislocation {v} that hits the (now partial) 
kink is not a commutator, and the relation {a'} — {a}{v} is no longer a trivial relation, in a way. In the above case it 
was possible, at least in a thought experiment, to abolish the node by smoothly turning the two segments to the same 
homotopy class - abolishing simultaneously the full kink. This is now forbidden. On the other hand, the attached 
defect exists only on one side of the master line, because {a 1 } = {a}{v} is now the TS expression of a true Kirchoff 's 
relation for three defects meeting at a node. In the Volterra process language, {a'} = {a}{^} is nothing else than 
Eq. 35. 

V. FOCAL CONICS IN SMECTIC A PHASES AS QUANTIZED DISCLINATIONS 

The most original defects in SmA phases are focal conic domains (FCD), whose geometrical properties were first 
investigated by Fricdel and Grandjcan (1910) and Fricdcl (1922). A FCD consists of a pair of confocal conies (an 
ellipse E and an hyperbola H) , which are the focal lines of the set of normals to a family of parallel smectic layers folded 
into Dupin cyclides (Hilbert and Cohn-Vosscn, 1964). For a recent account of the physics behind this geometry, see 
Klcman and Lavrentovich (2003); Kleman et al. (2004); an essential property of E and H is that they are disclination 
lines. 

Consider the simple case where E is degenerate into a circle C and thereby its conjugate conic C is a straight line 
orthogonal to the plane of C, passing through its centre. The Dupin cyclides are then nested tori (Fig. 24). They are 
restricted in this figure to the cyclides with a Gaussian curvature of negative sign, with planar continuations outside. 
This is the most frequent, if not the only one, empirical occurrence of toric domains. 

It is apparent that C and C are both wedge disclination lines, C of strength k = h, C of strength k = 1. According 
to the analysis in II. G, there are no attached dislocations, only attached disclinations, whose density can be written 

Sf = |nr5s = 2nr)tf (36) 

The Nye's edge dislocations that follow the latitude lines of the tori (call them of the first type) are of the type that 
attends such disclinations; their (infinitesimal) Burgers vectors Sh rot are along the meridian lines; they obviously 
introduce extra matter that curves the C disclination. On the other hand the Nye's edge dislocations (of the second 
type) that follow the meridian lines of the tori, whose (infinitesimal) Burgers vectors Sh tr are along the latitude lines, 
relax the quantized dislocations that attend the C wedge disclination, after the manner already discussed for a straight 
k — \ line, Fig. 21. Note that the Burgers vector 8h tr is variable along a dislocation of the second type, since the 
local frame of reference is continuously rotated by the dislocations of the first type. 
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FIG. 24 Toric FCD 




FIG. 25 FCD with positive and negative Gaussian curvature surface elements From Kleman and Lavrentovich, Soft Matter Physics, 

an Introduction, Springer, 2003, with permission 

All these results extend to more general FCDs. Figure 25 represents such a FCD with positive and negative Gaussian 
curvature surface elements. 

In a local frame of reference where the 1- and 2-axes are along the lines of curvature of the Dupin cyclide, the 3-axis 
along the normal, the dislocation density tensor (defined in Appendix A) has components 
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where and are the geodesic curvatures of the curvature lines. The geodesic curvatures vanish in the toric case, 

Pgl Pg2 

since then the lines of curvature are also geodesic lines. We recall that otij measures a dislocation content; the integral 
J J ctijdSi over an area bound by a loop 7 is the j-component of the Burgers vector of the dislocations going through 

this area along the i-direction. The dislocation tensor being a tensor is invariant under a change of coordinates; the 
choice of the curvature lines as coordinate lines has a simple physical interpretation in terms of <5b tr and Sh rot , for 
the components a\% and a 2 i- The extra dislocation densities CH31 and a.32 have also Burgers vectors parallel to the 
smcctic layers, as they must, for the symmetry reasons already alluded to. They correspond to edge lines along the 
normals to the layers, and contribute to the stability of the shape of the conies and their relaxation. 

As most generally observed (Kleman et ai, 2004), FCDs are restricted to the negative Gaussian curvature parts 
o\ 02 < of the cyclides, Fig. 26. The FCD is therefore confined inside a double cylinder lying on the ellipse, whose 
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FIG. 26 Generic FCD with negative Gaussian curvature layers, see text 

generatrices are parallel to the asymptotes of the hyperbola. Any layer inside this double cylinder is bordered by a 
curvature line of the cyclides, i.e. a circle, Fig. 26. The plane II containing the circle is tangent to the cyclide, so 
that the outside extension of the layer can be along II, without any cusp. II is perpendicular to one or the other 
of the asymptotic directions of the hyperbola, thus the set of all planes II form two families of parallel planes that 
eventually cross on the plane of the ellipse. By limiting the corresponding material layers to the half space above or 
below the ellipse plane IIe, IIe appears as a tilt grain boundary of misorientation u — 2 sin -1 e, where e is the ellipse 
eccentricity. E itself is a k = | disclination line to which are attached the constitutive dislocations of the tilt grain 
boundary. See Appendix C for a detailed investigation of its characteristics. 

VI. GEOMETRICAL FRUSTRATION. ROLE OF DISCLINATIONS 

A. Geometrical frustration, a short overview 

The concept of frustration covers a number of inhomogeneous structures which are all describable in terms of 
defects, in fact disclinations; see Klcman (1989) for a review. The word itself was introduced by Toulouse (1977b) in 
the frame of the theory of antiferromagnetic systems, where there exist closed paths of atoms with nearest neighbour 
exchange interactions that cannot be satisfied simultaneously. 

1. Unfrustrated domains separated by defects 

By geometrical frustration (Klcman, 1985a, 1987), we mean an extension of the concept of frustration such that: a) 
it connotes systems where the short range interactions are so much dominating that they completely determine the 
local configuration, and b) it is possible to describe these interactions in geometrical terms — in a sense geometrical 
frustration is an extension of the notion of steric hindrance. The concept is therefore of interest when the local 
configuration is incompatible with long range Euclidean ordering, i.e. is non-cristallographic in the usual sense of this 
term. The frustrated medium is thereby divided into small, unfrustrated, domains, of size £ say, separated by defect 
regions where the short range order is broken. 

In three dimensions, it is fruitful to introduce a crystal template where the unfrustrated domains spread without 
obstruction, if such a description is feasible. Such templates, where this local order extends homogeneously without 
distortion, are necessarily curved, non-Euclidean, with Riemannian habit spaces of constant curvature. The decurving 
of the template into an Euclidean medium employs disclinations of these curved crystals; they delimit the unfrustrated 
domains of the actual medium. Geometrical frustration connotes the existence of a particular type of incompatibility 
resulting from the different interactions in competition. Disorder at scale £ does not prevent that at larger scales 
there exist correlations between the unfrustrated domains; the resulting frustrated medium can be cither a crystal 
with broken translations — if the weak long-range interactions take over at some scale — or a medium truly disordered 
at all scales > £ — if these long-range interactions are very small. Very similar but simpler approaches in one or 
two dimensions were developed earlier to describe approximate epitaxy and surface reconstruction of crystals, helical 
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FIG. 27 Double-twist configuration. From Kleman, Lavrentovich and Nastishin, Dislocations in Solids, vol. 12, Elsevier, 2005, with 
permission 

magnetic order, charge or spin density waves (cf Fricdcl (1977b) for an introduction). Chemists have also talked in 
the same spirit since before World War II of non alternate electronic structures of molecules and solids. 

We first present a short overview of frustrated media with emphasis on the presence of disclinations. In most of 
the examples below, the 3-dimensional sphere S 3 is employed as a template habit space. 



2. Covalent glasses, disclinations 

For covalent glasses, frustration originates in the constancy of the coordination number z, which can be incompatible 
with certain ring configurations. For example 5- and 7- membered rings do not generate Euclidean order. A large 
literature exists on the subject, starting with the celebrated model of Zachariasen (1932) for continuous random 
networks; cf. Mosseri and Sadoc (1990) for covalent frustration models. 

As pointed out by Rivicr (1987), the underlying geometrical structure of a covalently bound material is a graph 
with a constant coordination at each node, except for dangling bonds and possible double bonds. There is no reason 
why covalent bonds should form polyhedra. A stacking of polyhedra is a particular graph. In that latter case, Rivier's 
theorem, which states that 

odd-mem bered rings are not found in isolation, but are threaded through by uninterrupted lines which form closed loops 
or terminate on the boundaries of the specimen (Rivier, 1979), 

is quite useful. But it applies also to non-polyhedral structures, as soon as rings are recognizable. Rivier's lines are 
obviously disclinations. See also Toulouse (1977b). 



3. Double-twisted configurations of liquid crystal directors and polymers, disclinations 

The most studied phenomenon of frustration in liquid crystals is the double-twist molecular arrangement of short 
molecules met in Blue Phases, and whose presence has been advocated in cholesterics of biological origin (DNA, etc.), 
see Livolant and Bouligand (1986); Livolant (1987); Giraud-Guille (1988). The chromosome of dinoflagellate (Livolant 
and Bouligand, 1980) displays also a double twisted geometry of DNA molecules of a type somewhat different from 
the type met in Blue Phases (Friedel, 1984; Kleman, 1985b). It has also been advanced that double-twist is present 
in amorphous molten polymers (Kleman, 1985a). 

Consider a liquid crystal of chiral molecules. If the director is an axis of cylindrical symmetry, all the directions 
orthogonal to any director can act effectively as axes of helicity — in a classical N* phase, this symmetry is broken 
and there is only one axis of helicity. The local unfrustrated arrangement can be described as in Fig. 27, where the 
integral lines of the director are helices of chirality opposite to the chirality of the rotation of the director about the 
radii. We have n r = 0, = sintp(r), n z — cosip(r), with ip(0) = 0. The molecules rotate with the inverse pitch 

2ir dip sin2ijj 

q(r) = — = -n-V><n= — + — ^ 

p ar 2r 

which is the same order of magnitude as the inverse pitch of the cholesteric phase. Double- twist is entirely satisfied 
only for the director along the z-axis (r = 0); at a distance r = p/A from the z-axis, the director has rotated by tt/2 
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along the r-axis and is now perpendicular to the z-axis; its vicinity is no longer double-twisted; frustration sets in. 

Blue phases are made of elements of double- twisted cylinders of matter assembled in space. Three cylinders can 
stack along 3 orthogonal directions; the region of highest frustration, in-between, may show up a singularity of 
the director field, a k — — | disclination. This local arrangement is found in several cubic symmetries showing a 
3D disclination segment network, (Barbct-Massin et al, 1984; Meiboom et al, 1983). The blue fog is amorphous, 
and the disclination lines arc seemingly at random; it has also been suggested that the blue fog is icosahedral, as 
in quasi-crystals (Hornrcich and Shtrikhman, 1986; Rokhsar and Scthna, 1986), achieving thereby another type of 
frustration. 

Double-twist can be thought of as resulting from a competition between a tendency to dense packing, insuring 
parallel alignment of the integral lines, and a tendency to chirality. But a geometry with equidistant helices, which 
would result from such a competition (Klcman, 1985a), is not homogeneous in Euclidean space. 

The template proposed by Sethna (1985) in a S 3 curved space insures homogeneous unfrustrated double-twist; we 
denote it {dtw/S 3 }. The director is along a family of great circles of the habit 3- sphere, all those great circles being 
equidistant and twisted, and parallel in the sense of spherical geometry. This is described in Appendix D (the Hopf 
fibration) and studied for its defects in VII. C. See also Dubois- Violette and Pansu (1988). 

4. Tetrahedral and icosahedral local orders, disclinations 

For amorphous metals and Frank and Kasper phases, the origin of frustration is the tendency towards dense 
packing of equal or quasi equal spheres, representing atoms. 

a. Frank and Kasper phases. In complex metallic alloy structures, particularly those of transition metals, it occurs 
frequently that the structure is entirely determined by the requirements for sphere packing, i.e. that the atoms form 
tetrahedral clusters, and that the coordination polyhedra are triangulated. 

Frank and Kasper (1958, 1959) thoroughful topological and geometrical study of the crystalline structures submitted 
to such constraints shows — if one admits that the number of neighbors of an atom on the coordination polyhedra 
to which it belongs is cither Z5 = 5orZs = 6 — , that there are only 4 types of coordination polyhedra, with Z = 12, 
14, 15 and 16. Frank and Kasper distinguish the sites Z — 12 as minor sites, and the other ones as major sites. The 
edges which join neighbouring major sites form a skeleton; sites of Z = 14, 15 or 16 are meeting points of two, three 
or four bones. This skeleton is much simpler to study than the structure as a whole. The description of the Frank and 
Kasper phases in terms of a skeleton of bones (i.e. a network of line defects) is contemporary to the development of 
the concept of disclination in liquid crystals (due to Frank (1958) himself), but it was only later that the topological 
nature of these defects as true disclinations was recognized (Nelson, 1983a). 

The Frank and Kasper networks constitute a remarkable example where the main characteristics of geometrical 
frustration show up; the frustrated atoms are along lines which structure a sea of unfrustrated atoms Z = 12, and 
there is a typical distance between lines which here scales with the lattice parameter. The existence of the skeleton 
of major sites is not dependent on the existence of a periodic lattice, and the only necessary hypothesis is that the 
medium be polytetrahedral. 

b. Amorphous metals. Icosahedral order (Z — 12) is met in Frank and Kasper phases, but also amorphous metals, 
undercooled atomic systems (Bernal, 1959, 1964; Frank, 1950a), and quasi-crystals (Schcchtman et al., 1984). It is 
also valid for small clusters less than a few hundreds atoms; for metals and rare gases, cf J. Friedel (1977a, 1984). 

Bernal has used a polyhedral approach in order to analyze hand-made systems of equal spheres, and has shown 
that a large majority of the polyhedra (ca. 86 %) are tetrahedra. Hence the predominance of local icosahedral order. 
Furthermore, these tetrahedra arrange frequently into pseudo-nuclei that are aggregates of face-sharing tetrahedra, 
two by two, and tend to build a connected lattice in the whole structure, wrapping themselves so-to-speak around the 
larger holes (i.e. rather low density polyhedra with V = 8, 9 or 10 vertices). Of course the tetrahedra cannot fill in 
the whole 3-dimcnsional space; however, since they can extend freely along one direction, one notices a large number 
of three-stranded spirals, right- or left-handed, formed by a one dimensional array of regular tetrahedra. Their local 
density is large. They are reminiscent of the twisted great circles of S 3 , and therefore fit locally quite well into the 
curved space crystal representative of the local order. The pseudo-nuclei are the regions of less frustrated order. 

Numerous analyses of atomic packings have followed Bernal's pioneering work; we refer the interested reader to 
Zallcn (1979) 's review on dense random packings. 

c. The {3,3,5} template. This is an example of a crystalline structure in a curved space. There are no icosahedral 
Euclidean crystals, but icosahedral order is compatible with a space of constant curvature, namely the 3-sphere S 3 , 
which can be tiled with regular tetrahedra, 20 of them meeting at a vertex, i.e. generating an icosahedron. This 
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structure is known after Coxeter (1973), as the {3,3,5} polytope: 
elementary facets have 3 edges, hence are equilateral triangles, 
3 facets meet at a vertex, hence elementary cells are regular tetrahedra, 
5 cells share a common edge, hence an edge is a five- fold axis. 

This polytope has No — 120 vertices, Ni = 720 edges, N2 — 1200 faces, and N3 — 600 cells. The passage from this 
template to a disordered system, i.e. the decurving of {3,3,5} and its infinite extension to a Euclidean space, occurs 
through the introduction in the perfect {3,3,5} crystal (Klcman and Sadoc, 1979), (Klcman, 1989) of disclinations 
of negative strength — i.e. which introduce extra matter — forming a 3D network in physical space, like the F&K 
network. For a description of quasicrystals in terms of frustration see Klcman and Ripamonti (1988) and Klcman 
(1989, 1990). 

B. The decurving process 

1. Rolling without glide and disclinations 

Let us now focus our attention to the transformation of a curved template into an actual flat medium, under the 
constraint that the local order of the template is conserved. An isometric mapping continuous all over the medium is 
not possible, because it would require that the Gaussian curvatures be equal at corresponding points (Darboux, 1894; 
Hilbert and Cohn-Vossen, 1964; Singer and Thorpe, 1996). 

However, isometric mapping can be achieved locally, by parallel transport along a line L, according to E. Cartan 
(1963). Letting M (S 3 , say) roll without glide upon E 3 , along any path L c S 3 , lengths and angles along the path 
are conserved at corresponding points in E 3 . If L is a geodesic of S 3 , it maps along a straight line L in E 3 : this is the 
so-called Levi-Civita connection. In general, such a mapping transforms a closed line LcM into an open line in E 3 . 
The closure failure can be described as a disclination. 

Consider a material cone: the curvature is concentrated at the apex, a useful feature here. For it maps isometrically 
on the plane as a whole, except at the apex, which becomes the vertex of an empty wedge bordered by two generatrices. 
This is clearly the picture of the Volterra process for a disclination, whose angle, called here the deficit angle, is a 
measure of the concentrated curvature. In order to complete the mapping, it suffices to fill the void with perfect 
matter. The final object does not carry stresses, if the material cone is amorphous; but it cannot be so if the medium 
is ordered, since then the disclinations are quantized. 

We show in the next paragraph that the decurving process yields generically two sets of disclinations: 

(i) those resulting from the mapping of M onto the Euclidean space E 3 ; these disclinations are of negative (resp. 
positive) strength if M has positive (resp. negative) Gaussian curvature; 

(ii) those resulting from an clastic relaxation (disclinations of a sign opposite to the former ones) . 

The relationship between disclination lines and curvature has been emphasized long ago by Kondo (1955-1967) and 
Bilby (1960). Their approach is opposite to the present one; Kondo and Bilby start from E 3 , which they consider as 
the habit space of the physical crystal, and map it on a space that is curved due to the presence of disclinations in 
E 3 . In contrast, the physics of the disordered system, which lives in E 3 , is contained in its curved representation. 

2. The Volterra process in a curved crystal 

It is useful to approximate a Ricmannian manifold by a piecewise flat manifold. Such a process of 'triangulation' 
has been proposed by Rcggc (1961) to calculate properties of curved manifolds in general relativity without using 
coordinates. For instance, in {3,3,5}, the edges of the lattice are replaced by straight lines in the embedding euclidean 
space, and the faces and cells by euclidean faces and cells. The edges, on which all the curvature is now concentrated, 
form the skeleton of the triangulated manifold (the 1-bones), which is articulated at the vertices (the 0-bones). 
Consider a D = 2 example, {5,3} (three pentagons at each vertex); its Regge image is a dodecahedron with flat 
pentagonal faces. By mapping the pentagons around a vertex (where the curvature is concentrated) onto the plane, 
one lets appear a deficit angle, : 

cj+ = 27t-3x^ = ^ (37) 
5 5 

Fig. 28 which is not equal to the angle of a quantized disclination in {5,3}, namely ft = 3n/5. The stress field produced 
by the Volterra process at such a vertex is that one of a negative disclination of strength w_ = £1 — lo + = — 2ir/5. 
This angle also measures the deficit angle of the local negative Gaussian curvature generated by the introduction of 
extra matter; this negative Gaussian curvature would manifest itself as a locally hyperbolic surface element, if the 
disclinated {5,3} is allowed to relax elastically in 3 dimensions. 
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FIG. 28 Deficit angle at a merging vertex of three regular pentagons, lu+ — j 

It is energetically unfavorable that the only disclinations present in the actual medium be of negative strength. 
Therefore all the vertices are not the seats of mapping disclinations, and those which are not are flattened by force, 
which yields stresses characteristic of positive disclinations of strength uj + . 

In 3 dimensions, the curvature is concentrated along the 1-bones; if one moves a vector by parallel transport on a 
closed circuit L about an edge i, it rotates by an angle u>i (the deficit along i) 

LOi = 2ir - 2J a k (38) 

which does not depend on the precise location of L as long as it is traced in cuclidean space and does not cross another 
edge; a k is the dihedral angle of the local flat polyhedron k with edge i. One can perform a Volterra process, now 
truly reminiscent of the usual Volterra process in flat space, either by gluing the 2 lips across the angle u>i, or by 
inserting (removing) a lattice unit cell f2 in the space separating the lips. According to the case, one introduces a 
topological disclination of strength ±|wj| or =p|wj — f2|. These disclinations are wedge disclinations, since they lie along 
a rotation vector n. The Volterra process for twist disclinations does not add or remove matter; therefore it does not 
change the curvature. 

The deficit vectors a>i that are merging at a vertex and are parallel to the corresponding edges obey the Kirchhoff 
relation 

X>i = 0, ( 39 ) 

here written for small WjS. This expression is in fact the Bianchi identity for curvatures when the deficit angles are 
small (Regge, 1961). 

To summarize, the actual medium contains two sets of disclinations, call them D s -lines when they carry spherical 
curvature and are positive strength disclinations, D/j-lines otherwise. Notice the duality between D s - and D/j- lines; 
it means that one can start equally from a spherical crystal or from a hyperbolic crystal to construct a disordered 
system (Klcman, 1982a, 1983). The precise location and density g s and Qh of the D s - and D/j- lines is of course subject 
to a great arbitrariness, and it is the best elastic balance which decides of the final choice. From that point of view 
an amorphous solid with local icosahedral order is certainly closer to a spherical crystal than an hyperbolic crystal, 
and one expects g s > Qh- 

C. The concept of a non-Euclidean amorphous medium 

The foregoing model of an amorphous medium, based on the existence of a frustrated order, does not forbid the 
conceptual possibility of a homogeneous, isotropic, structureless, medium in S 3 , whose defects we can investigate after 
the manner of the amorphous medium discussed in II. 

The model we have in mind has no relation whatsoever with any kind of local order. We denote it {am/S 3 }. 
Because of the curvature imposed by the habit space S 3 , the singularities that break the continuous rotation and (non- 
commutative) translation symmetries are somehow at variance with disclinations and dislocations in flat Euclidean 
space. This is also true for the {3,3,5} case, but the continuous case is expected to be easier to grasp. Furthermore, 
the concept itself of a curved amorphous medium is new and worth investigating in its own right. As we have just 
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observed, a pending question in the {3,3,5} case is how to decurve such a template in order to get an atomic liquid 
or amorphous medium with icosahedral local order. The difficulty lies in the fact that the {3,3,5} disclinations are 
quantized, and the answer is not unique; stresses remain. The same question for {am/S 3 } yields a unique answer with 
no stresses remaining, since defects are continuous in strength and distribution. 

Finally there is the question of disclination networks, analogous in spirit to Frank dislocation networks. Disclination 
networks are apparent in polynanocrystals (III.C) and in Frank and Kasper phases (VI. A. 4). They might also be 
important in undercooled liquids, but the true local geometry is that one of a spherical crystal with icosahedral 
symmetry (VI. A. 4). This situation obviously requires reconsidering the Kirchhoff relations in a curved habit space, 
either amorphous (as a generalization of II. F. 2, see VILE) or icosahedral (i.e. with quantized disclinations, see VII. F). 



VII. DEFECTS IN 3-SPHERE TEMPLATES 

A. Geometry and topology of a 3-sphere; a reminder 

A point M on the 3-sphere S 3 of unit radius will be defined by its Cartesian coordinates in E , {xo,xi,X2,X3} or 
more concisely by the unit quaternion x 

x = x + Xi« + x 2 j + x 3 fc|x| 2 = xx — Xq + x 2 + x\ + x 2 = 1 (40) 

where x stands for the (complex) conjugate of x. 



1. The rotation group SO(3) in quaternion notation 

The set of unit quaternions forms a group that is related to the group of rotations SO(3) in E 3 , as follows. The 
unit quaternion x 

| ■ 4 XlZ+X 2 j+X 3 fc 

x = cosi? + gsini?, q = —===== (41) 

y/x{ + x l 2 + x\ 

is representative of a rotation <& along the direction {xi,X2,X3}; two antipodal points x, —x on S 3 embedded in 
E 4 are representative of the same rotation ?? and d + tt along the same direction xi,X2,x^; q is called a pure unit 
quaternion; its real part vanishes, q 2 = — 1, qq = 1. All SO (3) rotations are therefore represented by a sphere S 3 with 
antipodal points identified, namely P 3 = S 3 /^, the projective plane in three dimensions. Reciprocally, Q = S 3 , the 
multiplicative group of unit quaternions, is the double covering of P 3 and, as a topological group, is isomorphic to 
SU(2) and homomorphic 2:1 to the group SO(3) of all rotations that leave the origin fixed, 

SO(3) = SU(2)/Z 2 , (42) 

the kernel of the homomorphism being generated by the rotation 2ir. 

Because of the validity of Moivre's formula for unit quaternions, x can also be denoted 

x = exp|g. (43) 



2. The rotation group SO(4) in quaternion notation 

The quaternion notation provides an easy analysis of the basic isometric transformations of E 4 (Coxeter, 1991) that 
conserve a fixed point (the centre of S 3 ). These are: 

(i) the single rotation, with one pointwise fixed 2D plane, the so-called axial plane containing the origin O; this 
plane is the E 4 -generalization of the rotation axis in E 3 . For that reason, one shall often call the single rotation in E 4 
as the rotation. 

(ii) the double rotation (with one fixed point only, the centre of S 3 ), which is the commutative product of two 
rotations about two completely orthogonal axial planes. 

a. The single rotation. The basic formula for a rotation of angle a about the axial plane H\^ w — (0, l,w), defined 
by three points in E 4 , namely a) the origin {0,0,0,0} —denoted {0}, b) {1,0,0,0} —denoted {1}, c) {0,wi,W2,W3} 
—denoted {u>}, is: 



a a 

x — e 2 xe% 



(44) 
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The transformation expressed by Eq. 44 leaves invariant any point x of IIi !tu , i.e. x — A + /iw (A, /i real), and no 
other points. 

If a; is a pure quaternion (not necessarily a unit one), Eq. 44 gives the expression of the rotation of its representative 
point x by an angle a about the axis w, on a two-sphere S 2 of radius |x|. Therefore a pair of two conjugate unit 

quaternions (e~ 2 W ,e2 w ) represents one element of SO(3); the pair (— 2 w , — e 2 *") represents the same element; the 
representation is 2:1, as already noted. 

Hi, w intersects the habit sphere S 3 of {3,3,5} or of {am/S 3 }, radius R, along a great circle Ci jt0 , radius R, which 
is thereby pointwise invariant in the rotation given by Eq. 44; {R} and {Rw} belong to this great circle. 

The vector at tangent to Ci tW at M is the local rotation vector in S 3 induced by the single rotation in E . 

With Eq. 44, we have considered a special axial plane. A rotation of angle a about a generic axial plane passing 
through the origin in E 4 is of the form 



p =/= ±q, with axial plane: 



' x e 2 9 (45) 



n M = (0,1 -pq,p+q). (46) 



It is easy to show that 1 — pq and p + q are invariant in the transformation given by Eq. 45. Notice that p and q 
are pure unit quaternions. 

This is the place to introduce the plane 

II^ = (0,1 +pq,p-q) (47) 

which is the plane completely orthogonal to n Pi(3 ; the directions denoted by 1 + pq and p — q arc both orthogonal to 
1 — pq and to p + q. A rotation of angle (3 about lip g taken as axial plane reads: 

x' = e2Pxe2 q (48) 

Since the four directions 1 — pq, p + q, 1 +pq and p— q are mutually orthogonal, they can be used as the directions 
of a Cartesian frame of reference in E 4 . 



b. The double rotation. Right and left helix turns. The rotations about two completely orthogonal planes are 
commutative. Consider the product of two rotations of the same angle a about H p >q and np~ ? . According to Eq. 45 
and Eq. 48, we have 

a a a a 

x' = e-2 v {e2 p xe2 q )e2\ 



i.e. 



i' = ie"« (49) 

This double rotation conserves S 3 (and any three-sphere centered at the origin) globally, and leaves only one point 
invariant, the intersection of the two planes n Pjg and lip q , i.e. the centre {0} of S 3 . It leaves no point invariant in S 3 ; 
thereby it is akin to a translation in Euclidean space; E. Cartan introduced the term transvection to connote such an 
operation in a Riemannian space (Cartan, 1963). 

p does not appear any longer in Eq. 49, hence this transformation can be given a geometric interpretation in any 
pair of a large set of completely orthogonal plane pairs. For instance, it is a transformation that, in the 2-plane Tli. q 
containing the origin {0} and the directions {1,0,0,0} and q = {0, q 1; q 2 , q 3 }, is a rotation of angle a, and in the 
completely orthogonal 2-plane 11^ is a rotation of the same angle a (Montesinos, 1987). 

It can also be called a right helix turn, for the following reason. Consider the transvection x' — xe m . In order to 
visualize it, let us employ the stereographic projection of S 3 from its pole {—1,0,0,0} onto the hyperplane spanned 
by the directions i, j, and k (the 'imaginary' part of the quaternion set); x 1 — xe m turns the equator plane spanned 
by the directions j and k (i.e. the plane 14^ completely orthogonal to III j) to the right by an angle a about the 
axis i, and pushes it forward by the same angle, Fig. 29(a). The same is true, also in stereographic projection, for 
the transformation x' — xe aq , where q is any pure unit quaternion chosen as the axis of the rotation, since any axis 
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FIG. 29 (a) Right helix turn x' — x e aq and (b) left helix turn x' = e ap x in stereographic projection 



q results from the axis i by a rotation in the 3D (i,j,k) hyperplane; all pure unit quaternions are visualizable in 
stereographic projection, as well as their corresponding equator planes. 
The left helix turn 

x' = e~ ap x (50) 

a a a a 

results from the product of two rotations of the same angle and of opposite signs, x' — e~ 2 p (e~ 2 p x e 2 q ) e~ 2 9 , 
performed on two completely orthogonal planes; it turns the equator of p also to the right by an angle a, but pushes 
its equator plane backward by the same angle, Fig. 29(b). A left helix turn and a right helix turn are commutative, 
but two left (or right) helix turns are not. Therefore these transformations are akin to non-commutative translations. 



3. Group of direct isometries in the habit 3-sphere S 3 

Any product of a right translation by a left translation is a direct isometry in S 3 : 

x' = e' f3p xe aq ; (51) 

it is also the product of two commutative rotations of angles a + (3 and a — f3 about two completely orthogonal 
axial planes n p 9 and 11^ ; a right (or left) helix turn is of the form exhibited by Eq. 51, with (3 (or a) = 0,2ir. 
{—e ap , —e~P q } generates the same isometry as {e ap , e~@ q } . 

Any isometry in S 3 being generated by a combination of two commutative helix turns {e ap ,e~P q }, one right, one 
left — each belonging to the group SU(2) — the maximal group that leaves S 3 invariant is 

SO(4) = G(S 3 ) - SU(2) x SU(2)/Z 2 = S 3 x SO(3) (52) 

It is a direct product. Its universal cover is G(S 3 ) = S 3 x S 3 , i.e. the group of unit quaternions, squared. 
A reminder of some geometric characteristics of S 3 can be found in Appendix D. 



B. Disclinations and disvections in S 

The notion of line defect easily generalizes to any ordered medium in S 3 . In terms of the Volterra process, we have 
two types of defects: 



1. Disclinations in S 3 

on 

Disclinations break a single rotation in E , defined by the SO(4) group element {e 2 p , e 2 9 } about some axial plane 
n Pjg ; it conserves the habit sphere S 3 globally. The plane H p , q intersects S 3 along a great circle C. For an observer 
confined to the habit sphere, this rotation appears as a set of rotation vectors fit tangent to C all along, as already 
indicated. The points x of the cut surface are displaced by the single rotation 

_n n 
x' = e 2 p xe2 q . 
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Remember that any great circle is a geodesic of S 3 , so that this operation is clearly a generalization of a disclination 
in E 3 , where the rotation vectors are along straight lines, i.e. Euclidean geodesies. 

2. Disvections in S 3 

Disvections break transvection symmetries and are thereby the generalizations of dislocations, which break trans- 
lation symmetries. A right (say) transvection 

X^ — Xi 6 ^, \Xi\ R. 

brings Mj^i) to M'j(a^), at a distance |6| = \x\ — Xi\ = 2 R sin ^, along the great circle C, through Mj and M'j. 
M'i in turn is brought by the same distance |6| along the same great circle. Any point Mj outside likewise follows 
another great circle Cj which is equidistant to C^, with the same |b|, i.e. which turns about C, in a double-twisted 
manner. The points x of the cut surface are displaced by the double rotation 



This operation is clearly a generalization of a dislocation in E 3 , with b the analog of the Burgers vector. 

The term of disvection was introduced to connote line defects that break non-commutative translation symmetries 
in quasicrystals (Klcman, 1992). 

C. Defects of the double-twist S 3 template 

The topological classification of line defects in the double- twist template {dtw/S 3 } is the same as that one of the 
uniaxial nematic (Scthna, 1985), namely 

nx(V dtw ) = Z 2 , (53) 

since Vdtw — „ TT f3 'L — = P2, the projective plane. In this expression, 
SU(2) x D 2 h 

SO(4)=SO(3)xSU(2) is the group of isometries of S 3 (Eq. 52), 

SU(2), in the denominator, is isomorphic to the group of left or right transvections — the template has a definite 
chirality, 

D2/1 is the rotation point group, as in a N phase. These results can be obtained using the geometric picture of the 
Hopf fiber bundle (see Fig. 36 in Appendix D), which is the geometrical representation of {dtw/S 3 }. 

Whereas the topological stability analysis provides only one class of topological defects, namely \k\ — ^, the Volterra 
process provides many more, exactly as for the N phase. Therefore one has to differentiate k = — \ disclinations, 
which add matter, and k = | disclinations, which remove matter. In principle, only the first category of disclinations 
are liable to decurve {dtw/S 3 } to an Euclidean medium; this is precisely the result which is claimed in the current 
BP structural models (Meiboom et at, 1983). However a caveat is in order: because viscous relaxation operates at 
constant density, a negative sign for the decurving disclinations should not be a prerequisite; Scthna (1985) has indeed 
proposed a model where there are no disclinations at all, with however the same cubic structure. 

D. Continuous defects in a 3D-spherical isotropic uniform medium 

1. The wedge disclination 

First, a few remarks about the great circles of S 3 , which are geodesies of S 3 , i.e. the equivalent of straight lines in 
flat space. 

a. Wedge disclinations are along great circles. Let H pq = {0, 1 — pq,p + q} be the plane that intersects the habit 
sphere along the great circle C. We have |1 — pq\ 2 = \p + q\ 2 — 2 — pq — qp, where — pq — qp = 2p ■ q = 2 cos 4> P ,q, § v , q 
being the angle between p and q; (1 — pq) and (p + q) are orthogonal. We then have, for any point u on C, 

u = a(l — pq) + b(p + q) \u\ = R 
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We choose the origin of the angles in the H p , q plane along the direction (1 — pq); one then gets: 

w(i9) = i?N Pi(3 el p (l -pq)e^ q = i?N M [(1 - pq) cosd +(p + q) sintf] . (54) 

where 

N p , 9 = (2 - pq - qp)- 1 ' 2 = (2 cos %*)-*. 

The first equality expresses the fact that the points of C are obtained by a rotation {e2 p ,e2 9 } about 

- - 

= {0,1 + pq,p — q} in E . Observe that the rotation {e2 p ,e2 9 } keeps pointwise invariant the intersec- 
tion of Iip q with S 3 , which is a great circle C 1 " orthogonal to C. The angular distance between any pair of points 

M e C, M 1 - e C , is equal to tt/2 (the two vectors OM and OM 1 are orthogonal). C and C , which have the 
same center and the same diameter, have no point in common. 

_ n n 

b. Volterra elements of a wedge disclination. Consider now the rotation {e 2 p 7 e2 9 } about n p 9 ; it conserves S 
globally. The plane H p , q intersects S 3 along C. In S 3 , this rotation is a set of rotation vectors fit tangent to C all 
along. It is therefore a rotation that builds C Volterra-wise as a line of wedge character, strength f2. Its constitutive 
dislocations are as follows. Let x = x" + x 1 - be the quaternion representation of a point P = x" + x 1 - on the cut 
surface Ec of C (it is not useful at this stage to select a special cut surface), split into its two components, one x" 

belonging to H p , q — invariant in the rigid Volterra rotation {e 2 p j e2 ? } that moves apart the two lips of the cut 
surface — , the other one x^, belonging to Hp q , not invariant in the rotation under consideration. Hence 

n n _n n 

.t'H = x 11 = e 2^x11 e 2 x ,J - = e 2P x ± e^ q . 

Thus, la;- 1 - 1 is the distance of P to the axial plane il p ^ of the Volterra process we are considering. 

The sum total of the Burgers vectors of the constitutive dislocations between the wedge line and P is, in quaternion 
notation, 

_ n n 

b ~ x' — x = e 2 p xe2 q — x, (55) 

which can also be written 

& = sin^e 2 p (xq-px) = sin ^ (xq-px)e^ q . (56) 
But, because x" belongs to H p , q and x 1 - to LTp- , we have 

ft ■& i9 

/ = e2Ve2', X W = e -2P x \\ e 2 q , 

for any Hence 

x ± q + px ± = x^q -par 11 = 0. (57) 



and eventually 



b= -2sm^px ± e2 q = 2sin f e" 2 P x x q. (58) 



To make contact with the discussion about Eq. 1 , we make use of two remarks 

(i) the modulus of b, calculated from Eq. 58, is 

|6| = 2|ar L |sin§ (59) 

(ii) the Burgers vector b is perpendicular to H p , q - Notice that it is also orthogonal to x 1 - and to x ,J -, which both 
belong to U^ q ; hence b OP = b OP' = 0. 

These two properties belong also to the constitutive dislocations of a wedge line in {am/E 3 }. Note that Eq. 59 
reproduces the integral of Eq. 1 performed on the interval between the line and a point P on the cut surface for a 
wedge disclination, since {x^l is in both cases the distance to the invariant manifold (in Eq. 1 it is the distance to the 
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rotation axis of the disclination.) The fact that b is orthogonal to the axial plane compares to the fact that in 
Eq. 1 the Burgers vector is orthogonal to the axis of rotation. 

The Burgers vector db of the constitutive dislocations situated between x and x + dx on the cut surface is 

_ n n 
db=e 2 p dxe2 q -dx, 

which can also be written 

db= -2sinf pdx ± e% q = 2sm i je'^ p dx ± q. (60) 

To summarize: Eq. 60, valid in {am/S 3 }, is the equivalent of Eq. 1, valid in {am/E 3 }. It describes the constitutive 
dislocations of a wedge disclination located at the intersection of the habit 3-sphere and the 2-plane H p ,q = {0,1 — 

pq,p + q}- 

It remains to find the manifolds along which dh is not only constant in magnitude, but also in direction, i.e. the 
geometry of the constitutive dislocation lines. The relation x 1 - = constant defines a two-plane parallel to n Pi<? , which 
intersects the habit sphere along a circle C^j., of radius which is not a great circle, except if x 1 - = 0, in which 
case C^-l is the disclination C itself. C and C x ± are parallel, perpendicular to some direction d x± belonging to 
6, which depends only on x 1 - , is constant on C x ± , which is indeed a constitutive dislocation of the wedge disclination. 
The set of constitutive dislocations C x ± for < \x^\ < R, so chosen that these circles are in parallel planes, all 
perpendicular to the same direction d x ± in 11^ q , describes half a great sphere, which is a geodesic manifold of S 3 , and 
consequently the equivalent of a half-plane in fiat space. There is a one-parameter family of such half great spheres, 
depending on the direction d x ± in Hp g , each of them playing the role of a possible cut surface £c for the wedge 
disclination C. 

Equation 58 can be given an interpretation in terms of disvections. We introduce the quaternion z± — —px = x q; 
it results geometrically from a rotation by an angle of ^ of the vector x 1 - in the ILj plane, so that we also have 

_ ZL i — 

z = e 4 p x e 4 q ■ z is constant all along C x ± . Hence 

db = 2sin^(dz ± e^ q )=2sm^(e-^ p dz ± ) (61) 
appears either as a left screw or a right screw. 



2. Defects attached to a disclination in {am/S 3 } 

a. Useful identities and relations. 

We introduce two great circles at each point M of a generic disclination line L; one of them C m is tangent to L at M 
along the direction m, the other one C M tangent to the local rotation vector i~2 = fi/x; mm = fxfi = 1. Their analogs 
in {am/E 3 } are the straight lines along m — the tangent to the disclination line — and fj, — along the rotation vector. 
The 2-planes to which these circles belong are denoted n p g = {0, 1 — pq,p+ q} for C m , H e . a = {0, 1 — gcr, g + a} for 
C^ t ; p, q, Q, o~, are defined as in Eq. 62 and 63. 

1 - 1 - 1 - 1 - 

p = um — — mu, q = — um = mu, (62 

R R R R y ' 

1.1- 1- 1. 

8 = ~R U ^ = R^ Ul ° = R Ufl = ~~R^ U ' ^ ' 

The demonstration is given in Appendix D. Notice that the relations mu + um — and um + mu = 0, easily 
deduced from Eq. 62, express the orthogonality of m and u, m ■ u = 0. Likewise, p,u + ujl = 0, etc. 

Observe that m is invariant in any rotation about n Pjg , and that \x is likewise invariant in any rotation about n ejCT ; 
this yields the relations 

_■& §_ _■& 
e 2 p me2 q = m, e % e [ie2 a = u, 

for any i9; hence: 

pm — mq = 0, gu~fia — 0. (64) 

These relations also stem from Eq. 62 and 63, but it is worth retrieving them this way, in order to emphasize the 
role of invariance by rotation. Since M belongs to both axial planes n Pj9 and Hg^ , we have 
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# ± _<& 

e 2 p ue2 9 =u, e ^ e ue^ a = u, 

for any hence: 

pu— uq = 0, gu — ua = 0. (65) 

also steming from Eq. 62 and 63. 

More generally, the quaternion components a" and a 1 - , resp. in H p , q and 11^ „, of any vector a = a" + a 1 - obey the 
relations 

pj - J q = 0, pa 1 + a 1 5 = 0, (66) 
and similar relations for its components ay and a± in H e ,a and 11^ CT , 

gan — an a = 0, ga±+a±a^0. (67) 

By differentiating, one also gets: 

± ± ± 

e2 p udqe2 q — udq e 2 p dp it e 2 9 = dpu 



e2 e uda e2 a = uda e2 e dgue2 a = dgu 



e 2 p mdqe2 q = mdq e 2 p dpme2 q — dpm 



e2 e ^ da e2 CT = uda e.i e dgue2 a = dgu (68) 
All these relations are easy to establish directly, by employing identities of the type 



e2 f dg = dge^z da e2 CT = e~ '2° do : (69) 

which express rotational invariance about n^ g and n^ CT : 
the quaternions 

udq, dpu, mdq, dpm 

belong to the plane 11^ q ; they are invariant in any rotation about the axial plane Hp q . 
the quaternions 

uda, dgu, uda, dgu 

belong to the plane II^; they are invariant in any rotation about the axial plane n^ CT . 

More generally, the quaternions b" and b 1 - in H pq and 11^ , and the quaternions by and b± in n e o - and n^ ff , 
b = foil + b 1 - = 6|| + bj_, are such that 

b" dq, dpb\ b X q, p b ± , belong to the plane 11^ , 

b^ dq, dpb^, b" q, pb", belong to the plane H p , q . (70) 

by da, dgb\\, bj_a, g bj_, belong to the plane n^~ CT , 

b±da, dgb±, by <r, g by , belong to the plane H e ,cr- (71) 
Other useful relations obtain by considering small rotations d§ about an axial plane. For instance we have 

d0 §0 

u + du = e2 p ue2 q , 
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which yields: 

du = \(pu + uq)dd. (72) 

Since pu — uq — (Eq. 65), Eq. 72 yields 

du = pudd = uqdt) — Rm d$, (73) 
the last equality also originating in the expressions for p and q, Eq. 62. 



b. A general expression for the attached defect density. 

We reproduce the extended Volterra process approach we have used for Euclidean crystals. Because {am/S 3 } is 
amorphous, any disclination L, whatever its strength, can be constructed as a sum of infinitesimal constitutive defects, 
attached to L if L has twist character. The wedge disclination case has already been discussed; the constitutive 
defects are disvections (one set of disvections, whose 'chirality' is ambiguous, left or right, Eq. 61). Disvections are 
the analogs in {S 3 } of dislocations in {E 3 }, so that one can consider this result as the {S 3 } analog of the result in 
{E 3 }. Contrariwise, as we show, in the twist or mixed line case, the attached defects are gcnerically disclinations, 
which however can be defined as the sum of two sets of disvections of opposite 'chiralities'. 

Now, assume that L is any loop in the habit 3-sphere, and consider two neighboring points M and M + <iM on 
this line, u and u + du in quaternion notations, with uu = (u + du)(u + du) = R 2 , i.e. duu + udu — 0. In order to 
compare the Volterra processes at two neighboring points M and M + dM, we proceed as follows. 

Let f2 = fi/i. be the rotation vector of the disclination at M. Being a rotation vector, fi belongs to the axial plane 
n e ,cr = {0,1 — ga, g + er} that intersects the habit sphere along the great circle running through M, and to which 
fi is tangent. This is akin to the situation investigated in the case of the wedge disclination, where C M is in fact the 
wedge line itself and where (j, is tangent to the wedge line (observe that fi does not belong to the habit 3-sphere but 
to the 3D flat tangent space to the sphere at M). We therefore have fi ■ OM = 0. 

Consider now two close points M and M + dM on the disclination L, with rotation vectors along u and u + du, 
and a point P (x in quaternion notation) on the cut surface El- The variation in displacement observed from M to 
M + dM at the same point P is: 

db M = e~% {e+de) x ef (<T+d<T) - e~% e xe% a = sin § (e~% Q xda - dgxe^°). (74) 

If the line L is a great circle and the local rotation vector is along the tangent of this circle, which means that 
one can choose g, a constant (independent of the point on C), then dg = da = ; the line is of wedge character, as 
expected, and db~M = . There are no attached defects. We retrieve the results of VII.D.l. 

Equation 74 is the fundamental equation related to attached defect densities in {am/S 3 }. 

We show now that it cannot be interpreted the same way as Eq. 4 in {am/E 3 }, although the extended VP from 
which it results is similar. 



c. Attached disclination densities. 

As a matter of fact, Eq. 74 appears as the difference between two disclinations of axial planes H e . a and Ii g +dg,a+da] 
this difference can be expressed as a unique disclination carrying an infinitesimally small angle of rotation: 

The quaternions g+ dg and a + da are pure unit quaternions, if second order terms in da 2 , dg 2 , etc. are neglected; 
dg and da are pure quaternions, of equal moduli \dg\ — \da\ — dX, real, positive. We introduce the two pure unit 
quaternions t g , t a such that dXt g = dg, d\t a = da. 

Let us define 

d M = e~ 2 f e xe 2 ta -x. (75) 

Equation 75 is the expression of the displacement d^/i{x) at P(x) due to an infinitesimal disclination of angle df2p 
about the axial plane H tg .t a - Notice that dyi{u) — 0, which confirms that the just defined infinitesimal disclination 
is attached at M to L. We now show that Eq. 74 can be given the same form as Eq. 75, with the following choice of 
dflp: 

dfl F = 2 sin § cos § dX, (76) 

justified below. Equation 76 can also be written cK7p = 2ci(sin-2-), dflp thus appearing as the differential of the 
modulus of the Frank vector introduced in II. F. 2, if one assumes that the variation d£l of the argument equals 
2 dX sin ^; we see later on why it should be so. 
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We now give hints how to identify Eq. 74 and Eq. 75, and find passim an interesting simplification of the expressions 

of d&M, ^M- 

Split x = xy + x± into its components xu € H e , CT and xj_ G n4" CT . According to Eq. 71, xy der and dgxy belong to 
_n Q 
il^ a . Hence e 2 e xy (icr = xy der e 2 CT , and: 

d& M = sin ^ ({xy dcr - rfgx||}e 2 CT + {e 2 ^ xj_ da — dg xj_ 6 2 ^jj (77) 
Because xy € n ecr , we have 

xy a — gx\\ =0. 

Likewise, because xj_ £ we have 

Xj_ (7 + £X_L = 0. 

We differentiate these relations, keeping x constant, and get: 

xy da — dgx\\ =0 x± da + dg x± = 0. 

Hence 

_n n 
rf^M = sin y {e 2 e x_l rfcr — dgx± e 2 CT }. (78) 

The terms depending on xy have disappeared. The remaining terms can be transformed in many ways, using the 
equalities just derived and the rotation invariances of Eq. 70 and Eq. 71. We get for db^i, 

dbui = -2 sin § cos § d\ t e x± — 2 sin § cos § dA xj_ £ CT , (79) 

and using similar transformations for gJm, 

c?m = —dfl-p t g x± — dflF x± t a . (80) 

Hence, to summarize: g?&m and describe the same infinitesimal disclination attached to L at M, of angle dH,F, of 
axial plane n 4e)t<r . 

Notice that the definition of 5£l — 2 dX sin ^ is not mere chance; it originates in expressions of the type e 2 ( a + dcr ) ; 
which can also be written 



For the sake of completeness, we write some other expressions for db^ 

n 



db M = sin§{-e 2 e + x ± due 2 CT }, (81) 



We shall make use of Eq. 81. 
Also, because x± 6 n^ CT , 

we can write 



where v± £ n^o-i because 



o — - — CT 

: sin -J {— dg e 2 e x 1 + xj_ e 2 a da}. 



Q Si 

7TO FT <T 

e! ii = x± e 2 . 



o^M = sin ^ (— dgv± + ux do"). 



no n n n 

e~2e(e2e x± )e-2< T = x± e - ^ = e 2^x A 



(the first member of this equation is indeed a rotation of v± with axial plane CT ), 

d. Infinitesimal Burgers vectors and disclination lines. 

Starting from Eq. 74 or Eq. 81, we calculate |d&M|i using the relation |cK>m| 2 — d&M^M- The calculation, not 
reproduced here, uses some of the equalities established above. One finds: 

|cJ&m| = 2 sin 2 § \xj_\d\ (82) 
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This equation compares to Eq. 61 (wedge disclination) by the common presence of the \x±\ term — again, the 
relevant distance is the distance to the axial plane that carries the rotation vector of the disclination. The presence 
of the sin 2 ^ term is interpreted as follows. We can split d&M into two Burgers vectors, appearing in Eq. 81, related 
to two disvections, one right, one left, 

Q _Q 

dbM,cr = sin |r x j_ da e 2 CT , d&M, e = — sin ~ e z e dgxj_. (83) 

db-M.g and d&M,o- have the same Burgers vector modulus, namely |d6M,g| ~ |^m.<t| = sin ^ IxjJdA. The extra 
sin 7^ factor in | cZfojvi | means that the two vectors dh^/i g and dbjyr.o- make a constant angle, equal to fi. One can check 

directly that dbM,p ■ = \ {db^ji tB dbm,a + d&M, e di>M j0 .) = cosfi. 

The description of the infinitesimal defects related to the disclination L in terms of disvections is equivalent to the 
description in terms of attached disclinations. Notice that the disvection lines, which are those lines along which 
db-M.g and g?6m,o- are constant, are not attached to L — which would require that x± = 0, since the points of L are 
characterized by the quaternion coordinate u which obeys u± = 0. The line L itself is a particular disvection line, 
of vanishing Burgers vector. Henceforth the infinitesimal disvection lines form two sets of lines surrounding L, each 
comparable to the set of disvections surrounding a wedge disclination line. 



3. Twist disclination along a great circle 

By constant rotation vector we mean that the local rotation axis of the disclination L, namely in quaternion 
notation, with /ifl = 1, is parallel-transported along a line with the natural connection of the sphere S 3 . If such a 
line is a great circle, i.e. a geodesic, the rotation axis makes a constant angle with this circle, fi ■ m = cos if. This 
configuration is obviously equivalent to a straight disclination line in {am/E 3 }, with constant rotation vector. A point 
M on L depends on the variable u or the arc angle #. Let be the great circle tangent to at M(#) € L; it 
belongs to the plane H g . a = {0, 1 — g a, g + a}; g = g(t?) and a — a(i)) are given by Eq. 63. n p . 9 = {0, 1 — p q,p+q\ 
is the two-plane that contains the great circle L; pandg are constant. 

The parallel transport of p, along L, a geodesic, is also a rotation about the axial plane ELjg = {0, 1 + pq,p — q}, 
completely orthogonal to H p , q - Such a rotation can be written: 

fi + dp, = e 2 p pe 2 q , 

which yields 

d[i=\ (PM + pq)d$. (84) 
Using Eq. 84 and Eq. 63, the following expressions for dg and da obtain: 

dg = | (pg — gp) di), da = —^(qa — aq) (85) 

We also have 

dX = \ dd\p g — gp\ = g d'Q \q a — a q\ 
We have n ■ m = cos tp — | (m fi + p,rh), a constant. We can show, again using Eq. 62 and Eq. 63, that 

(P Q - QP) 2 = (q & ~ & l) 2 = 4 sin 2 tp, 

which yields 

d\ = \ sin <p\dd. (86) 

As indicated by the presence of the coefficient | sin. cp| in the expression of dX, the only component of \x which is 
relevant in the twist properties of the line is its component orthogonal to m. Hence we split the rotation vector as 
follows 

M^H/^ + m". (87) 

where p,^ belongs to the plane and is therefore invariant in any rotation about this plane; = m cost/? belongs 
to the plane H p , q - It is a question of simple algebra to show that the /zy component does not contribute to (a q — qa) 
and to (p g — gp) ■ The only component of the rotation vector that contributes to the dislocation densities attached 
to C is (which is a constant all along the disclination), through the |sin^j| factor. But notice that the rotation 
vector is still £1 fi. 

The same discussion as in VII. D.l and VII. D. 2 is valid, with the simplifications of Eq. 85 and Eq. 86. 
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E. Kirchhoff relations 

Any isometry of S 3 has of the form x' — e~P q x e ap , Eq. 51, and can be split into the product of two commutative 
helix turns, one right, one left. Therefore a product of isometries can be written 

^ = l$...l$^... x ... r £H) r «... r g>, (88) 

where l^^r^ e Q = SU(2). n disclinations meeting at a node can be split into two sets of disvections meeting at 
this node, n left disvections, n right disvections. Thus, we first investigate the Kirchhoff relations for (left or right) 
disvections. 



1. Three disvections meeting at a node 

Let ft,Q = e ap , h^Q = e^ q , h*Q = e ir be the elements of symmetry of three right (say) disvections meeting at a 
node. We then have 

hgh$>h®>={±l} (89) 

Remember that the group of all unit quaternions is 2:1 homomorphic to the group of all the rotations of a 2-sphere 
that leave the origin fixed. In that sense, /iq and —Hq represent the same rotation in SO(3), but they do not represent 
the same transvection in SU(2). {—1} transforms x into —x by an helix turn of angle 7r about any axis w: {—1} = e™. 
It however will appear as a fundamental necessity to introduce the quaternion {-1} in Eq.89, as we shall see. 

Denote by a, b, c the angles between the directions defined by the three pure unit quaternions p, q, r, namely 
a = Z(q, r), b = Z(r,p), and c = Z(p,q). Equation 89, which also reads, e.g., h^? = ±^q , yields: 

- three relations coming from the real part of this relation, of the type 

cos a cos (3 — cos c sin a sin (3 = ± cos 7. (90) 
where cos c = p • q = p 1 q 1 + p 2 q 2 + P3 <73 , 

- nine relations coming from the pure quaternion part, of the type 

qi sin (3 cos a + p\ sin a cos /3 + s c ,i sin c sin a sin (3 = =Fri sin 7 (91) 

where the unit vector s c (i.e. the pure unit quaternion s c ), defined by the cross product s c sine = (p x q) has been 
introduced; similarly s a sin a = (q x r) and s& sin b = (r x p). 

Multiplying Eq. 91 by s Cj i and summing over the three Eq. 91 containing s c 1; s c 2 , s c ,3i one gets 

sin a sin sin c Vp,q,r 
sin 2 a sin 2 (3 sin 2 7 sin a sin (3 sin 7 

where V p>qir = r • (p x q) is the scalar triple product; V v . q , r — P • s a sin a = q • S& sin 6 = r • s c sine. 

One recognizes in Eq. 90 and 92 expressions much akin to those met in 2D spherical trigonometry. The geometric 
interpretations of Eq. 89 in terms of spherical triangles are different, whether one has a { — 1} or a {1} in the 
right-hand member. 

a. Ag>fcg>Ag> = {-!}. 

One finds that there are two conjugate spherical triangles, hereunder denoted T p r and T Sa Sb 8< _, whose angles and 
angular arcs yield { — 1} in the right-hand member: 

(i)- the vertices of T~ q r are the extremities P, Q, R, of the 3- vectors p, q, r; the angles of the triangles are a, /3, 7, 
the angular arcs are a,b,c. See Fig. 30. Standard results in spherical trigonometry yield (Wcisstcin, 1999): 

— cos 7 = cos a cos (5 — cos c sin a sin (3 (93) 

(94) 



sin a sin f3 sin 7 Vp,q,r 



sin a sin b sin c sin a sin b sin c 
Equation 93 obtains by applying the 'cosine rules for the angles', specialized to the angle 7 
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FIG. 30 Spherical 2-D representation of Kirchhoff relations for disvections. P, Q, and R are on the sphere centered at O. See 
text 

The construction of a triangle of angles a, (3, 7, of angular arcs a, b, c, requires that be obeyed the following inequal- 
ities 

a + /3 + 7^7r, a + b + c^2n (95) 

(ii) the vertices of T s ~ Sb Sc are the extremities A,B,C, of the 3-vectors s a ,st,,s c ; the angles of the triangles are 

tt— a, n—b, 7T— c, the angular arcs are tt— a, TT—j3, tt—j. One checks that s„ ■ sj = (q x r) ■ (r x p) = C ° S sina S sinfc OSC e * c -' 
i.e. s a • Sf, = — cos 7, etc., by applying the 'cosine rules for the angles', specialized to the angle tt — c. Equation 93 is 
also satisfied (it obtains by applying the 'cosine rules for the edges', specialized to the edge of arc tt — 7), and Eq. 94 
is replaced by 

sin a sin& sine V s s . s ,„„. 

sin a sin (3 sin 7 sin a sin j3 sin 7 

Eq. 94 and 96 are identical and both identical to Eq. 92, as can be shown by employing the standard identities 

Vp,q,r = sin a sin 6 sin 7 = sin b sin c sin a = ••• , 
V Sa ,s b ,s c = sin a sin /3 sin c = • • ■ , 

which yield 

Vp,q,r Vs ai s b ,s c = sin ol sin f) sin 7 sin a sin b sin c 
Vp qr = sin a sinfe sincy Sa . SbiSt; , 

V s a ,s b ,s c = s ' ma s ' m P sinjV p , qt r. 

The construction of a triangle of angular arcs n — a, tt — (3, it — 7, of angles n — a, it — 6, 7r — c, requires that be 
obeyed the same inequalities as for T~ r . 

Remark. The two tetrahedra OABC and OPQR are in conjugate positions on the 2-sphere; the edge s a of OABC 
is perpendicular to the facet (OQR) of OPQR, and the edge p of OPQR is perpendicular to the facet (OBC) of 
OABC, etc. Hence the appearance of the arc angles tt — a, conjugate to the angles a, etc., of the angles tt — a, 
conjugate to the arc angles a, etc. 

Both geometrical representations OABC and OPQR of the Kirchhoff relation for three disvections, are equivalent; 
we retain the first one. They express that the set of disvections terminate on a disvection of angle tt, whose direction 
W is not given a fixed value. In that sense, the terminal node can be interpreted as a singular point. 

Notice also that the same {—1} disvection can be assigned to the spherical triangle (PQR)* outside the smaller 
triangle (PQR). This shows that, in the present representation of transvections, the full 2-sphere S 2 has to be 
assigned the identity disvection {+!}, with angle 2tt. 

b. 4g>Ag>ftg>-{l}. 

Notice that e 7173 = e vq = e 7rr = { — 1}. Hence the {+1} Kirchhoff relation of Eq. 89, for the angles a,/3, 7 and 
the angular arcs a, b, c, can be transformed to the folllowing one 

e (7r-7)r e (Tt-0)q e {TT-a)p = (97) 

which can be discussed as the previous one for the angles tt — a, tt — (3, tt — 7 and the angular arcs tt — a, tt — b, tt — c. 
The inequalities of Eq. 95 are now replaced by 

a + b+c^TT, a + /3 + 7^27r (98) 
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FIG. 31 The two sketches represent conjugate disvections. See text 

1. e. the complementary ones to the case { — 1}. 

For the spherical representation of Eq. 97, see Fig. 30. Now the relevant angles are the external angles. There is 
no { — 1} disvection at the point where the three disvections merge. 

2. Orientation vs. handedness of a disvection 

The handedness of a disvection is a topological concept; its orientation is related to the Volterra process, it does 
not make sense topologically. 

We haven't yet taken care of the orientations (of the disvection lines, of the edges of the triangles T^ r ) and of the 
signs of the angles a, /3, and 7. In the usual dislocation theory, the orientation of a dislocation line is fixed arbitrarily, 
from which orientation the sign of the related Burgers vector is deduced, but still depending on a convention, generally 
the so-called FS/RH convention (Nabarro, 1967). For a given orientation, the change of Burgers vector has a topolog- 
ical meaning, because b and — b are different translation symmetries. We have here a somewhat more subtle situation. 

a. Topological considerations only. 

(i) The change of sign of a transvection, viz. x — * xe ap changed to x — > — xe ap , does not connote the same 
disvection. For in such a case the triple node e ap e l3q e~ ,r — {±1} changes sign: e ap e l3q e~ ,r = {t!} is a different 
disvection. Hence a change of sign does not correspond to an arbitrariness in the orientation. 

(ii) The change of sign of the angle of a transvection, viz. x — > x e ap changed to x — ► x e~~ ap , does not connote the 
same disvection. The triple nodes e ap e^ q e ir = {±1} and e~ ap e~@ q e~ jr ^ {±1} are different disvections, because 
of the non-commutativity of the transvections. A change of sign of the angle does not correspond to an arbitrariness 
in the orientation. 

(iii) Two conjugate disvections, viz. x — > xe ap and x — > e~ ap x are equivalent disvections. They are not equal, 
because they differ by the disclination e~ ap x e~ ap , which breaks a proper rotation, but a proper rotation does not 
modify chirality in the usual sense. This emphasizes the concept of conjugacy (change of handedness). The two 
equations e ap e l3q e ir — {±1} and e~ 7r e~@ q e~ ap — {±1} are conjugate. The two oriented triangles of Fig. 31 
(illustrating the case {—1}) are conjugate triangles. 

b. Volterra process. 

The topological considerations above indicate that the lips of the cut surface of a disvection are displaced at x by 
the quantity b(x) = x e ap — x (expressed as a quaternion) , but there is no hint whether this process adds or removes 
matter. The choice being made, it can be related to an arbitrary orientation and a convention on the sign ±|b|, in 
the manner of dislocations. The situation is then much alike, and is independent of the handedness. 

The situation is different with disclinations, because the opposite rotations e~ ap xe ap and e ap xe~ ap are both 
allowed topologically. 

3. More than three disvections meeting at a node 

Let 

e ap e f3q e^ r e Ss = {1} (99) 
be four disvections meeting at a node, P, Q, R, S the terminations of the directions p, q, r, s on the unit 2- 
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sphere, a = Z(SPQ), f3 = Z(PQR), 7 = Z(QRS), 6 = Z(RSP). They are represented Fig. 32, top, as a spherical 
quadrangle that can be split into two triangles T~ q r and T~ g _. Equation 99 also holds if the triangles are of types 
T pqr and T+ s p , i.e., if the angles are all complementary to the inside angles of the quadrangle (not represented 
Fig. 32). The two other sketches of Fig. 32 correspond to other possibilities mixing inside and outside angles. 

More generally, a n-gon is the graphical representation of the meeting at a node of n disvections; the n-gon can be 
split into n triangles (in many ways), each of them of type T+ p . pk or T~ pjPk . 



4. Kirchhoff relations for disclinations 

It suffices to consider one set of right disvections and one set of left disvections, separately. To each set the results 
of the previous section apply. 

a. Three disclinations meeting at a node. 

Assume first that the right and left transvections in Eq. 88 are complex conjugate, Iq — fg , and that the angles 
cti and the directions Pi are such that they form a {—1} triangle. The same triangle represents as well the left and the 
right transvections, Fig. 31. Hence, because { — l} 2 = {1}, the left and the right transvections all together amount to 
an identity isometry. The representative triangle, which is at the same time 'left' and 'right', has then to be thought 
of as a {1} triangle, i.e. a disclination of angle 2tt. Taking into account the complementary triangle on the 2-sphere, 
the full S 2 has to be assigned the identity isometry {+1} 2 , with angle An. 

Also, the rotation I^q xTq (with Iq — fg ) is of angle 2aj, each disvection taking part by an angle on . Let us 
consider the Euclidean limit in the OABC representation. Keeping in mind that the rotation Iq xtq (i'q = fg'') is 
of angle 2a j, (a, = — half the angle of the related disclination), and noticing that, in this limit, the quantities sin a;, 

etc, are proportional to the edge lengths of a flat triangle, it appears that in this limit the vectors i ai = 2 sin pi, 
etc. obey a Kirchhoff relation, as already observed for disclinations in {am/E 3 }. 

b. Extension to the generic case, when Iq and Vq are not complex conjugate. 

(i) (i) 

It is possible, as noted in VII. A. 3, to split a rotation Iq xrQ into two rotations with complex conjugate transvec- 
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FIG. 33 Node relations for disvections and disclinations, mixed. A disvection is either left (black lines) or right (shaded lines), 
(a)- disclination split into two disvections of opposite hands; (b)- disvection split into two disclinations 

tions. Such a splitting having been performed, we are back to the previous case. 
5. Mixed case 

Figure 33 indicates on two simple examples how disclinations and disvections can merge at nodes, in the particular 
case when the involved disvections all have the same strength, up to chirality (complex conjugacy). Notice that three 
disclinations of the same strength cannot meet at the same node, a seemingly obvious result, easily demonstrated by 
using the splitting of each disclination into two disvections of opposite hands. 

These results are entirely topological, and do not take into account the orientations and signs of the disvections 
that have been previously discussed (VILE. 2). For example, in Fig. 33(b), one can assume that the two bundled right 
disvections are of opposite signs and same 'Burgers vectors'. Therefore the bundle annihilates and the two remaining 
disclinations are oriented the same way, making them together a unique disclination. Hence Fig. 33(b) is not at 
variance with the result of VII. D. 2, established in the extended Volterra process fashion, which states that a twist 
disclination carries continuous disclinations, not disvections. 

F. The {3,3,5} defects 

A short reminder of the {3,5} and {3,3,5} groups is given in Appendix F. 

We distinguish between disclinations, disvections, and defects that combine both types. In this section, we restrict 
to defects in the non-decurved {3,3,5} polytope. 

1. Disclinations 

There are 60 Volterra disclinations, when identifying 'antipodal' elements in Y. This number goes up to 120 in the 
TS classification, for which IT ^{3^5}) ~ Yx Y. Each element of rotation symmetry leaves invariant a great circle 
C w . For instance, for \ui\ — 27r/5, follows a sequence of 10 spherical edges of the {3,3,5} polytope; in the locally 
Euclidean version of {3,3,5}, in which each spherical facet is approximated by a planar facet, each spherical edge by 
a straight segment in E , C w is a polygon with 10 edges. Observe that a wedge disclination is a disclination loop 
along such a G u . The local rotation vector art is along the tangent to G w \ no attached dislocations (disvections) are 
necessary to curve the disclination line. See Nicolis et al. (1986) for a quantitative description of a disclinated {3,3,5} 
crystal in its habit 3-sphere. 



57 



2. The disvection 'Burgers vectors' 

There are 60 left disvections Ziyi {1} = h^x, and 60 right disvections {l}x h=y = x /iy, when identifying 'antipodal' 
elements in Y. These numbers go up to 120 right disvections and 120 left disvections in the TS classification. 
Disvections leave no fixed point in {3,3,5}. 

Consider the right screw x e aq , x = xq + xi i + X2 j + X3 k being some {3,3,5} vertex M on the habit sphere S 3 of 
radius R. The displacement of M is Ax = xe aq — x, whose absolute value | Ax| is \x\ \e aq — 1| = 2i?sin ^. The left 
screw e~ aq x has the same Burgers modulus | A x\ as the right screw; it is a constant independent of x. However, the 
direction of Ax varies with M in {3,3,5}: it is tangent to the Clifford parallels belonging to the Hopf fibration related 
to q, Fig. 36, Appendix D. 

There are two such Burgers vectors at each point, one 'right', one 'left'. Thereby, Ax is not a Burgers vector in 
the ordinary sense, but it would tend to a constant vector if the radius R of {3,3,5} is allowed to increase without 
limit; however this process has no physical meaning, and we shall see later on another mean to transform Ax into a 
genuine Burgers vector. 

Taking 0:5 = J, we find | A x|s = iJr _1 , which is precisely the length of an edge of the polytope {3,3,5}. This 
'Burgers vector' is the smallest possible. It displaces x to one of its twelve nearest neighbors, which are at the vertices 
of an icosahedron. The orbit of x under the repeated action of the right screw is a geodesic circle C^^ight that is 
approximated by a ten-edge polygon. The coordinate xo measures along the 0-axis the distance to the geodesic circle 
C p . The angle 05 made by the orbit of x with C p is ar ^°' ) , (Sommerville, 1967). The same considerations apply to a 
left screw, yielding the same absolute value for the Burgers vector, but directed now along another great circle Cg^e/t 

at x, with angle ~ a ^ x °- > . Hence the spherical angle between the two circles Cs^igfot and Cg^e/t at x is arc jf ; which 
has to be equal to 27r/5 (five ten-edge great circles emanating from any vertex and belonging to a geodesic 2-sphere). 
Observe that the notion of right or left is relative to C p ; a C5 has an intrinsic meaning, there are 72 different Cg's, 
since there are 120 vertices and since 6 Cg's are incident at each vertex (6x120/10 = 72). There are 48 right (48 left) 
screws whose angle is a multiple of ir/5 (72 degrees of arc, 288 degrees of arc, 144 degrees of arc, 216 degrees of arc) 
along six possible directions from each vertex. The Burgers vectors have respective lengths i?r _1 , i?(3 — r) 1 / 2 , Rt 
and R(2 + r) 1 / 2 . They all join x to the vertices of icosahedra of respective edge lengths i?r _1 for the first and R 
for the three others. The Burgers vectors of length R(3 — r) 1 / 2 join x to its twelve 3rd nearest neighbors, which also 
form an icosahedron (Coxeter, 1973). 
We also have: 

I A x\s — R for 0:3 = ^ —this is the distance from the centre x of a dodecahedron (edge length Rt^ 1 ) to its twenty 
vertices, which are the next nearest neighbors of x; these are not the vertices of the {3,3,5} polytope, but the centers 
of the {3,3,5} triangular facets; 

I A x\2 — V2R for o?2 = § —this is the distance from the centre x of an icosidodecahedron (edge length \/2R) to its 
thirty vertices, which are the 4rth nearest neighbors of x; observe that the edge length can be larger than the radius 
of the habit sphere (Coxeter, 1973). These vertices are the midpoints of the {3,3,5} edges. 

3. Screw disvections 

In analogy with the Euclidean case, we define a screw disvection as a disvection line parallel to the Burgers vector, 
which in our case points out a great circle C„, n = 2, 3 or 5. We have emphasized above the n = 5 case. The cut 
surface can be any surface bound by C5, for example half a great 2-sphere. Now, again in analogy with the classical 
Euclidean case, we are interested in an ideal cut surface S s that slips along itself in the Volterra process. An ideal 
cut surface of a screw disvection (i.e., a loop) has to be parallel to the Burgers vectors, which are supported by a set 
of Clifford parallels. Several possibilities exist: 

(i) C p being the axis about which C5 is twisted, drop a segment of geodesic Cj_ (an arc of great circle) perpendicular 
to C p between any point of C5 and C p ; this operation determines a unique point on C p if C5 and C p are not conjugate, 
any point on C p if C5 = C p . 

(ii) the surface generated by the great circles Clifford-parallel to C5 and lying on Cj_. In effect, consider two 
conjugate great circles C5 and Cg, generated by a five fold axis; these conjugate geodesies are indeed ten-edge great 
circles, as shown in Coxeter (1991). Assume that the screw disvection line is along C5. The Burgers vector varies in 
direction along C5 and the cut surface, but has a constant modulus, as shown above. It is along the tangent to C5. 
The cut surface is made of Clifford parallels to C5 and Cg which stretch on the skew surface E s between C5 and Cg. 
The Burgers vector attached to each point of S s is by construction along the local Clifford parallel, and the Volterra 
process consists in a movement of E s in its surface. 
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4. Edge disvections 

An ideal cut surface S e of an edge disvection loop has to be orthogonal to the Burgers vectors, which are supported 
by a set of Clifford parallels. Therefore, according to the previous discussion, E e is a cap of a great sphere, bound by 
a loop of any shape. 

VIII. DISCUSSION 

The foregoing pages present a general view of disclinations and dislocations. The emphasis is put on 

- their interplay with continuous defects (themselves characterized as continuous dislocations, disclinations, or 
dispirations) in the constitution of defect textures in all media with continuous and frustrated symmetries, 

- their various relaxation processes as they can be approached by the consideration of defects, in particular contin- 
uous ones. 

These questions are not entirely new, as can be seen from the citations, but they have been largely ignored in the 
ill-ordered media community. The purpose of the present essay is therefore to stress the changes of perspective which 
occur in the theory of defects in mesomorphic phases (liquid crystals) and other 'ill-ordered' media, when the notions 
of disclination and continuous defects and the correlated use of the Volterra process, are fully used, in concurrence 
with the topological stability theory of defects. Central to that analysis is the concept of extended Volterra process, 
essential to the understanding of disclinations: this considers, in the last stage of the process, a viscous or plastic 
relaxation of the elastic stresses developed. 

The other approaches to the defects in ill-ordered media are of two sorts. Only the second one is directly 
comparable to the Volterra process approach. 

a. In liquid crystals with quantized translation symmetries, like smectics or columnar phases, a large part of the 
observations can be described in pure geometrical terms as isometric singular textures, in a way largely inspired by the 
pioneering work of Georges Fricdcl in the first decades of the twentieth century (Friedel, 1922; Friedel and Grandjean, 
1910); focal conic domains and developable domains, which obey rather restrictive geometric properties, are of the 
sort. But the geometrical approach to defects is limited to situations where the role of continuous defects can be taken 
aside. For a review, see Klcman et al. (2004). It is a rather remarkable fact that some mesomorphic liquid phases 
adopt, on rather huge sizes compared to molecular sizes, configurations that obey very precise geometrical rules 
and that these configurations are directly related to the structural symmetries of the medium. These configurations 
are worth being investigated for themselves, as they continue to be to-day. But these investigations did not inspire 
considerations leading to the present theory of defects. 

b. In contrast to the emphasis on the geometrical point of view, and since the late seventies, when the theory appeared, 
the defect classification basic concepts in mesomorphic phases and in frustrated media are those of the topological 
theory. The topological theory relates the stability of defects to the topological properties of the order parameter; it 
only concerns discrete, quantized, defects. That topological stability leads to energetic metastability is in many cases 
a reasonable assumption, which however cannot distinguish between different topologically equivalent configurations. 

A. The extended Volterra process 

1. Pure Volterra process in the absence of plastic relaxation. Constitutive defect densities 

The pure Volterra process has been thoroughly employed to the study of material deformation (Frank, 1950b; 
Friedel, 1964); it applies to the construction of a Volterra disclination in a solid, but only to a limited extent. One 
distinguishes two cases. 

a. Perfect disclinations. The only case is that one of a straight, wedge, disclination line in a Volterra continuous 
elastic body or in a crystalline solid, with rotation vector parallel to the line, SI being a symmetry operation of the 
crystal. 

b. Imperfect disclinations. This includes 

(i) straight wedge Frank grain boundaries, that are analyzable: in an amorphous medium, in terms of parallel 
infinitesimal dislocations, uniformly distributed, and in a solid crystal, in terms of finite dislocations, at least for small 
rotations, 
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(ii) straight twist lines or lines of a more general shape, which absorb/emit constitutive dislocations attached to 
the line or in its vicinity. In a solid medium with no relaxation, all these disclinations require very large energies to 
be created and also to move: their creation involves large stress concentrations extending over large regions; their 
motion would, whatever their nature, leave behind a stream of dislocations that could only be dispersed by slip and/or 
climb. These properties are related to the facts that, in such a solid without relaxation, the long range clastic energy 
is very large as soon as f2 is finite and the solid is of macroscopic size; the core energy is also very large as soon as 
the position of the lines deviates from the axis of rotation. Except for border cases involving disclinations of very 
small rotations, connected with weakly polygonized grain boundaries in crystals, one has not observed disclinations 
in solids as isolated objects, but only regrouped in close parallel strands with compensating strengths: parallel wedge 
pairs of opposite rotations in single crystals (described e.g. in Fricdcl (1964)) or triplets with rotations following a 
Kirchhoff relation, at the edges of the grains in relaxed polycrystals, as analyzed above in III. In this last case, the 
twist components of the involved disclinations exchange their constitutive dislocations. In the absence of possible 
relaxation by slip or climb, all these are sessile defects. 

Remark The pure Volterra process does apply only if O > — 2ir, since below this limit there is no matter left. But 
any value in the range so defined is allowed, even if the configuration is topologically unstable. Notice also that there 
is no limitation on the choice of the cut surface S with respect to crystalline symmetries, but this yields different 
dislocation configurations with different energies, e.g.: a tilt grain boundary if the rotation vector f2 (Eq. 1) is in the 
plane of S, with a set of parallel edge dislocations; a twist grain boundary if the rotation vector is perpendicular to 
the plane of E, with two crossing sets of parallel screw dislocations. 

2. Extended Volterra process. Relaxation defect densities 

Some plasticity of the medium can allow dislocations to relax. Disclinations are relaxed either by Nye's dislocations, 
or /and by absorbed/emitted attached defects. 

a. Nye dislocations compensate as best as possible the dislocations of the grain boundaries associated to the disclina- 
tion. They widen the singular core of topologically unstable disclinations (as they do for a k — n line in a nematic). 
They can be included in the extended Volterra process, in a last step of the process that takes into account the plastic 
relaxations allowed by the boundary conditions and the symmetries. 

Notice that the extended Volterra process allows f2 < — 2ir. 

b. Emitted/absorbed dislocations play a prevailing role when the disclination moves or changes its shape, in interplay 
with Nye's dislocations. And, as Nye's dislocations, they can be included in the extended Volterra process, also taking 
into account the allowed plastic relaxations. They can be induced: in a solid, either by high temperature diffusion, 
or by plastic deformation at low temperature (or twinning, or crack formation), in a liquid crystal by the anisotropic 
flow of matter. 



3. Mostly liquid crystals 

It has been known for long, well before the emergence of the topological theory, that the geometrical variability, 
the energetic (as opposed to topological) stability, and the elastic relaxation of defects in cholesterics (Fricdel and 
Kleman, 1969) is largely controlled by relaxation defect densities. This occurs through mechanisms, not yet enough 
investigated, that involve the continuous symmetries of the medium. A continuous defect , just like a quantized defect 
in the Volterra classification, is related to a symmetry element of the liquid crystal, a continuous symmetry precisely. 
In terms of the topological theory, it always belongs to the identity clement {1} of the first homotopy group IIi(V); 
in that sense, a continuous defect is never topologically stable. Continuous defects are often attached to quantized 
disclinations, whose flexibility and relaxation they control. Similar results apply to all mesomorphic phases, as they 
all display continuous symmetries. The extended Volterra process was present in Fricdcl and Kleman (1969), but 
has not been developed since. It is interesting to note that even the Georges Friedel's 'rigid' defect geometries get a 
different light when discussed in terms of continuous dislocations. 

Remark Some non-continuous Volterra disclinations also belong to {1}. For instance, in uniaxial nematics N the 
Volterra process provides two types of quantized disclinations, those whose rotation angle is an odd multiple of 
ir, and those whose rotation angle is a multiple of 2tt. Only the defects of the first type are recognized as being 
topologically stable (i.e. differing from {1}), the defects of the second type belong to {1}. This is so because the 
liquid crystal symmetries contain a continuous rotation subgroup (not present in solid crystals) whose presence drives 
the first homotopy group finite, whereas the group of Volterra defects is infinite. In practice quantized defects that 
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are not allowed by the topological theory are also of interest, as it may happen, for instance with suitable boundary 
conditions, that the energetic stability of the defect prevails over its topological stability. In this sense these quantized 
disclinations which carry an angle of rotation multiple of 2n are true disclinations. 

A deformed mesomorphic phase being a liquid, there are no strains at rest; but obviously there are large curvature 
deformations that can be analyzed in terms of Nye's densities (Klcman, 1982b), of continuous dislocations or continuous 
dispirations. Since continuous defects are topologically related to the continuous symmetries of the mesomorphic 
phase; as already stressed, this puts limitations to the possible curvature deformations of the medium in the vicinity 
of quantized disclinations, and to the shape modifications and mobility of those defects. 

Indeed a continuous transition between Hooke strain elasticity for solids and Frank-Oseen curvature elasticity for 
mesomorphic phases is insured by a continuously growing density of dislocations whose Burgers vectors decrease and 
vanish, in such a way that the total Burgers vector keeps constant. 

4. Extended Volterra process vs. topological stability 

a. The topological stability theory only considers defects that cannot be suppressed by plastic relaxation. It provides 
an a posteriori description that results from a mapping of the distorted medium onto the order parameter space 
V = E 3 /H. Quantized dislocation and disclination invariants {a^} belong to the first homotopy group of the order 
parameter space IIi(V), which is usually non-abelian. 

The topological approach has the advantage to be extendible to defects of any dimensionality (point, line and 
wall defects, and configurations). It is a more condensed process but a poorer one than the Volterra process, as it 
neglects the boundary conditions and makes equivalent all the configurations that can be deduced one from the other 
by continuous deformations. For such deformations to occur, it needs a high degree of plasticity which is not met 
in solids at low temperatures, whether crystals or glasses and other amorphous systems, but which is met in liquid 
crystals and ferromagnetism. It does not always allow to predict the stabler configuration in ordered media, still 
less to foresee the finer details of any configuration. By neglecting the topologically unstable situations, TS makes 
the bet that the corresponding singularities do reduce their energy by dissipating in some way. This bet is most 
often justified; however it might fail in the presence of special boundary conditions, as the k — 1 nematic line in a 
capillary tube (Cladis and Kleman, 1972; Meyer, 1973) and more generally in thin samples, where their broad core 
was well-known by G. Friedel (1922). Also, in the presence of strong material constant anisotropics, a k = 1 discli- 
nation core may be singular — this latter situation is met in nematic main-chain polymers (Mazclct and Klcman, 1986). 

b. The extended Volterra process is an a priori description of defects, as it gives a process of creation of line defects 
(only) whose characteristic invariants are classified by the elements of the symmetry group H of the medium) in a 
medium free of stresses. 

It yields the same conclusions as the topological stability theory, but at a finer level, by properly handling all the 
plastic relaxations, including those related to line-attached defects. This approach can be particularly useful when 
investigating dynamical aspects, when the viscosity is large as in solids, in most smectics, and in polymeric liquid 
crystals (Friedel, 1979; Klcman, 1984). It is also useful when dealing with nanostructures. This is at the price of an 
often much more complex analysis, as already stated. 

5. Flashback on a posteriori Volterra description of a defect in an amorphous medium 

A Volterra constructed disclination is qualified by the line direction with respect to the rotation vector: it is a 
wedge, or a twist, or a mixed disclination line. Similarly a dislocation is a screw, edge, or mixed dislocation line. Such 
a specificity does not belong to a line defined by its topological invariant. The question therefore arises whether it is 
possible to qualify any given line of the deformed medium in the same terms. This requires an a posteriori description 
of the defect, i.e. some kind of mapping that brings back the deformed medium to a stress-free, reference, medium, 
from which the line has been supposedly constructed in the Volterra mode. 

This program seems feasible when the reference medium is endowed with a 3D crystal lattice L whose equivalent 
directions and equivalent points can be recognized in the deformed medium. Directions can be followed along circuits 
surrounding disclinations, and a rotation or Frank vector obtained, which is unique in the reference lattice. This is 
how the Kirchhoff relations are established 6 . 



6 Notice that the Burgers vector of a dislocation surrounding a disclination fl is not uniquely determined this way; it is defined but to a 
rotation f2, as first noticed by Slecswyk (1966), see also deWit (1971); Harris and Scrivcn (1970, 1971). This is a small (and interesting) 
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These considerations, namely the possibility of fully characterizing an isolated dislocation or disclination, extend 
to the continuous theory of defects of Bilby, Kroner, etc, which always assumes explicitly the presence of a lattice, 
whose repeat distances are most certainly taken infinitesimally small, but this does not invalidate our reasoning. The 
continuous theory, in a way, contemplates continuous sets of isolated defects, whose Burgers vector is nonetheless 
infinitesimally small. For this purpose it employs the methods of differential geometry on manifolds. It is precisely the 
existence of a lattice that justifies the use of Cartan-Levi-Civita's distant parallelism method in the continuous theory 
of defects. For, if equivalent points and directionsdid not exist in the deformed medium, the comparison of vectors 
at two distant points would make no sense. It then follows that there are quantized disclinations in an amorphous 
medium, but their characteristics depend on the gedanken lattice drawn on it; therefore they do not present physical 
interest. We have considered only continuous disclinations in amorphous media, in this essay. 

B. Volterra processes in various media, compared 

The Volterra process applies to amorphous and frustrated media as well as to liquid, solid crystals, and templates 
of frustrated media. There are however some differences worth stressing. 

a. In an amorphous medium, the translation vector and the rotation vector that define the relative displacement 
of the cut surface are elements of the full Euclidean group E 3 . The related defects are not topologically stable. 
Dislocations and finite disclinations are non-independent defects. A finite disclination carries a field comprising 
constitutive and relaxation dislocations. It also possibly carries infinitesimal disclinations. The constitutive and 
relaxation dislocations relate to the strength of the line and to a part of its curvature (through the kink mechanism) ; 
the infinitesimal disclinations relate to another curvature component. These infinitesimal disclinations correspond 
to infinitesimal rotations belonging to the rotational part of the Euclidean group, whereas the strength of a finite 
disclination is characterized by the invariants associated to the constitutive and relaxation dislocations, i.e. belonging 
to the translational part of the Euclidean group only; these translations generate the Frank vector associated with 
the finite disclination. 

b. In a solid crystal, the characteristic Volterra invariants of the defects are the translation and the rotation elements 
of its Euclidean symmetry group H, whose elements classify dislocations and disclinations (necessarily finite) as, now, 
independent defects. However, the concepts of constitutive and relaxation dislocations still make sense. Because there 
are no infinitesimally small rotations in H, a line can be curved only by the kink mechanism, which makes use of finite 
relaxation dislocations. On the other hand, a straight disclination line is necessarily of wedge character, as there are 
no infinitesimal dislocations to give it a twist character. 

Notice that this description applies better to coarse-grained crystals with a quasi perfectly polygonized, annealed, 
structure. Most grain boundaries in a polycrystal are strongly misorientatcd, like polynanocrystalline materials below. 

c. In a polynanocrystal, the nanograms are separated by generally large misorientation grain boundaries which 
thereby are not analyzable in terms of quantized dislocations (those of the crystal) but rather in terms of continuous 
dislocations. In an ideal picture, three grains meet along a segment of (continuous) disclination, and these segments 
meet at quadruple nodes, forming then a 3D disclination network. The plastic deformation of polynanocrystals is at 
variance with that one of usual polygonized crystals; it is governed by the disclinations, which yields considerable 
stresses. Disclination-governed plastic deformation has long been a favorite subject of the Saint Petersburg russian 
school (see Romanov and Vladimirov (1992) for a review), without attracting much attention elsewhere. 

d. In a liquid crystal, the characteristic invariants of the defects are as above the translation and the rotation elements 
of a Euclidean group H, whose elements classify dislocations and disclinations as, again, independent defects. But, 
according to the case, the corresponding translation and rotation symmetries are finite or/and continuous, including 
infinitesimally small elements. In a nematic (N), translation symmetries are continuous, including infinitesimally 
small elements, and rotation symmetries are of both types. In a SmA, there are quantized and continuous translation 
and rotation symmetries. Line defect curvatures are therefore related to both the kink mechanism through (quantized 
or continuous) dislocations and the presence of infinitesimally small strength disclinations. 

e. In the curved habit spaces of frustrated media, like {am/S 3 }, {3,3,5/S 3 }, and {dtw/S 3 }, dislocations (which 



complication stressed in II.E.l, whose effect on the Kirchhoff relations has not been investigated. 
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we call there disvections) are not commutative, and disclinations can still be denned, with the same differences as 
above between amorphous media and crystalline media. Generically, disvections do not attach to disclinations, but 
disclinations do. Three-D disclination networks are therefore important features of the habit curved spaces, if one 
draws one attention to the plastic deformation of such spaces at constant Gaussian curvature. 

Disclinations also form 3D networks in actual, Euclidean, amorphous systems, liquid crystals such as cholesteric 
blue phases, Frank and Kasper phases, and probably in undercooled liquids or quasicrystals. In such 3D networks, the 
disclinations have to be somewhat flexible, which is possible, whether these disclinations are quantized or not, only if 
other defects, dislocations (defined in the Euclidean decurved medium) or disclinations, continuous or not, attach to 
them. 
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APPENDIX A: Continuous dislocations in solids and Nye's dislocation densities 

This is a very simplified presentation of a topic that has been much inspiring for the theory of continuous dislocations. 

We consider a deformed body described in terms of deformations /3y or of strains ey = |(/3y + we assume 
that these quantities are so small that it is licit to ignore second order terms. As it is well known, it suffices to know 
the strains ey in order to derive the stresses and the elastic energy, so that the deformations /3y are generally ignored. 
In the absence of elastic singularities (dislocations, disclinations), the /3ys derive from a displacement: /3y = Uj : i (we 
du ' 

use the notation Ujj = — — -). The conditions (due to Beltrami, 1888) to which the strains have to obey in order to 
oxj 

insure the existence of a displacement function are 

(inCG^jij = Sijk^-lmn^jm.kn — (Al) 

Here, we have introduced the incompatibility tensor (ince)y = £ijk^imn^jm,kn i it does not vanish in the presence of 
elastic singularities. 

The condition for the existence of a displacement function takes a much simpler form in terms of the /3y-s 

otij = ~e ik iPij,k = (A2) 

The theory of continuous dislocation densities emphasizes the role of the tensor ay ; when there is no displacement 
function, it does not vanish, and it can be interpreted in terms of dislocation densities; we refer for details to Kroner 
(1981), see also Kroner (1958) and Nabarro (1967), chapt. 1 and 8. 

The integral Aui — — J da Pijdxi = — JJ a otkjdSk on a loop da that bounds the surface element a vanishes if the 
strain derives from a displacement function (/3y = Ujj); if this is not the case, the 1-form fiijdxi is not integrable 
and the integral measures a displacement vector Auj. Au is interpreted as the Burgers vector of the dislocation 
densities that pierce the surface bounded by da; ay is the density of dislocations along the a^-axis; it measures the 
total Burgers vector component along the Xj-axis of the set of dislocations through the unit area perpendicular to the 
Xj-axis. For a dislocation L of finite Burgers vector b along a direction t, one gets ay = — tibjS(L). 

Now, again after Kroner, we write ct^ as follows: 

ay = -Eikieij.k + u>ij - LUk.kbij (A3) 

by splitting /3y — e.ij+ujij into a symmetric part ey and an antisymmetric part wy We have introduced uJi = \eiki ojm, 
(reciprocally wy — EijkUJk), which yields UJi — \sikiPu- Hence the relation A3. If is 'compatible', i.e. ay = 0, 
we have = \etki u^k- In the generic incompatible case, the /3ys and the w^s have to be given new interpretations, 
which we detail below. 

For the time being, notice that u>i measures the rotation of an element of volume with respect to the lattice directions 
(ct>3 is a rotation about the a^-axis with respect to the Xi-axis and the a^-axis); according to Eq. A3 it contains two 
contributions to the rotations, one coming from the strains and the other from the dislocation densities ay. This 
splitting into two contributions makes sense, because it is possible to conceive an ordered medium with no strains 
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FIG. 34 Genealogy of Nye's dislocations. These figures represent the gradual passage from a situation where (a) a random 
density of finite dislocations imposes curvature and strain to a situation where (c) the Burgers vectors becoming infinitesimal, 
the strain vanishes but the curvature subsists. In-between (b), the dislocations of (a) have organized into subgrain boundaries 



but a density of dislocations, the Nye's dislocation densities (Nye, 1953), at the origin of the rotation distortions. An 
example is given in Fig. 34, directly inspired by Nye's analysis. In a strainless medium, these densities read: 

&ij = - uj kyk Sij (A4) 

which can also be written 

Wjj = ay - \a kk 5ij (A5) 
In the general case, the gradient tensor of the rotations is 

= £ikietj,k + a ij - 2 a kkh,j (A6) 

where the two contributions (from the strains, from the dislocation densities) are separately non integrable. 

Of course, this model is valid as long as the rotations w^s are small. In the smectic case discussed in the text, the 
WjS measure the rotation of the orthonormal Darboux-Ribaucour trihedron from its equilibrium position along the 
principal axes of the shell (Darboux, 1894). 

There is no equivalence between a description in terms of Nye's dislocation densities and a description in terms of 
strain dislocation densities; in other words the solutions in (3ij of the equation onj = — e,-fcj (3ij, k — uiij — to kik Sij are 
such that flij + Pji = . Reciprocally the solutions in /Jy of the equation = —Ei k i f3ij. k = —sn-i eij, k are such that 
flij — f3ji = 0. Descriptions in terms of Nye's dislocations and of strain dislocations are exclusive one from the other. 

The gradual passage depicted in Fig. 34 can also be thought of as a gradual transformation of a solid (with quantized 
dislocations carrying strains) into a SmA phase (with continuous Nye's dislocations only). In this transformation, the 
distorted solid crystal elasticity gradually becomes the elasticity of the distorted smectic phase. More generally, the 
prevalence in mesomorphic phases of Nye's dislocations over quantized dislocations is the sign of their liquid character, 
and the origin of the Frank— Oseen elasticity (Frank, 1958) replacing the Hookean elasticity. 

Coming back to the physical meaning of the /Jys and of the u>iS in the generic incompatible case: 

/?2i describes a plastic distortion that displaces an element of volume with a finite section AxiAa^ along the xi-axis, 
by a quantity /^lAa^. This distortion can be managed by the glide along the same axis of a set of edge dislocations 
parallel to the X3-axis, with Burgers vectors along the xi-axis. There is no change of volume density. 

P22 describes a plastic distortion that elongates (or shortens, according to the sign) an element of volume with a 
finite section Ax±Ax2 along the £2-axis, by a quantity fiii^xi- This distortion can be managed by the climb, along 
the same axis, of a set of edge dislocations parallel to the X3-axis, with Burgers vectors along the X2-axis. There is a 
relative change of volume density equal to P22 ■ 
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^3, discussed above, can also be thought of as the result of two glide operations, foi an d —0X2- 
In the present theory, no account is taken of disclination densities, since the rotation dwi = u>i t jdxj is of course 
integrable: the density of disclinations vanishes. 

In fact this theory computes the minimum density of dislocations necessary to produce a given distortion of the 
medium. It assumes that this results from dislocation creation and annihilation, as well as from motion by glide and 
climb. Such hypotheses arc well fulfilled in magnets or liquids or (partly) in liquid crystals. They are not generally 
fulfilled in solids, except partly at very high temperatures. As in TS, it assumes in fact perfect plastic relaxations but 
for given boundary conditions. 

APPENDIX B: Topological stability and Volterra process compared. Conjugacy classes and homology 

To start with, notice that L i— > (0, a) and L" i— ► (2, a) are in the same conjugacy class of the first homotopy 
group; in fact, all the elements (2p,a), p <E Z, belong to the same conjugacy class and represent the same core-type 
of disclination; the other type of core corresponds to the conjugacy class (2q + l,o), q € Z. Each core-type of 
disclination, i.e. each conjugacy class can be identified in the Volterra classification with two types of disclinations, 
k = ij, differing by the sign: the topological classification does not distinguish k — \ and k = —\. There is an 
infinite number of other conjugacy classes, corresponding to the dislocations, each defined by a positive integer r and 
containing two elements {(r, e), (— r, e)}, i.e. two Burgers vectors equal in modulus, opposite in signs. All together 
one has oo + 2 conjugacy classes: 



Each conjugacy class therefore corresponds to an clement of the Volterra classification, although it does not specify 
whether the line is twist, wedge,... This assertion can be made more accurate as follows. The classes C 2r - they 
represent dislocations whose Burgers vector is even -, are commutators of Hi(Vs m A)', it is indeed not difficult to 
prove that any element of the form uvu v , where u,v = (n,a), belong to a C 2r . The commutators generate a 
normal subgroup of ^i{Vs m A) ~ the so-called derived group D — , such that IIi( Vs m A) can be partitioned into cosets Ci 



with Co = D; ci = (2p + l,a)D; c 2 = (2q, a)D; C3 = (2r + l,e)D. It is easy to check that the content of each coset 
is independent of the values of the integers p, q, and r. Furthermore, all the elements of C\ belong to ci, those of 
C 2 belong to c 2 , and those of all the C 2r +i belong to C3. Finally, Hi(VsmA)/^> = {co, c±, c 2 , C3} is an abelian group 
whose identity element is c . It is the dihedral group D 2 ; its multiplication rules C1C3 = c 2 , C2C3 = C\ reproduce 
the effects described above of the emission/ absorption of a dislocation with an odd Burgers vector at the core of a 
disclination. Two disclinations of different core types give an odd dislocation, when merging: cic 2 = C3. All the odd 
dislocations appear as one and the same element C3; the even dislocations appear in the identity element cq. More 
generally, IIi(V)/D is an abelian group, called the first homology group, noted Hi(V, Z). 

Let f^i be some rotation symmetry axis in a SmA liquid crystal; an axis f2 2 obtained by any rotation of f^i about 
the normal to the layers is also a rotation symmetry axis. However the two wedge disclinations Li and L 2 are not 
distinguished by the topological classification, which assigns indifferently the element (0, a) to both. The latter result is 
another weakness of the topological classification, which originates in the presence of a continuous rotation symmetry 
about the layer normal in a SmA. The same difficulty does not arise with a SmC phase, whose first homology group is 
also Z? 2 , whereas the conjugacy classes are twice as numerous as in a SmA phase; IIi(V smc) = ZOZ4, see Bouligand 
(1974); Bouligand and Klcman (1979). Observe that the symmetry point group of a SmC is HsmC = ZDZ 2 (reflections 
excluded), and that D 2 is precisely the quotient of the point group by its derived group. ZOZ2 is the Volterra group 
of defects (it is easy to check that it contains the two kinds of rotation axes, those in the layer planes, and those in the 
mid-layer planes, all perpendicular to the mirror plane), and D 2 can be understood as the Volterra group restricted 
to its abelian properties. The SmC point group is a quantized group, and the relation just indicated between the 
Volterra group and the homology group can be extended to any medium with a quantized point group (reflexions 
excluded). The same relation does not hold when the point group is continuous, as it is for SmA (Ks m A = ZODqc), 
sec Klcman (19821)), chapter 10. 
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FIG. 35 Distribution of the rotation vectors tt = 7rt on an ellipse in a focal domain. The layers intersect the plane of the 
ellipse along circles centered in one of the foci 

But in both cases the even dislocations are lumped in the identity element of Hi(V, Z). In both cases, eventually 
in all cases of ordered media, the Volterra classification has a stronger kinship with the set of conjugacy classes of 
ili(V) than with Hi(V, Z). The lumping of even dislocations in the identity element of the first homology group 
is of course a weakness of the homology classification; a tentative explanation is given in Kleman (1982b), chapter 
10, where can be found a more general discussion of the conjugacy classes and of the homology classification in any 
ordered medium(see also Kleman (1977)). 

APPENDIX C: The ellipse in a FCD as a disclination 

In a generic FCD (Fig. 26), the hyperbola H is completely embedded in the domain and is a wedge disclination of 
strength k = 1; indeed any direction along it is a 2tt symmetry axis. On the other hand the ellipse E, which is usually 
on the boundary of the domain, is a disclination of strength k — | for the full layer geometry, which contains the FCD 
and the two sets of outside layers that meet on the plane of E. But E is not in general of wedge character, because 
the tangent to the ellipse is not a folding axis for the local smectic layers. An obvious solution is given Fig. 35, where 
the local rotation vector fit, with = w, is tangent to the layers, which are folded symmetrically with respect to the 
plane of the ellipse (Kleman et ai, 2006). Other solutions are possible, with i~2 off this plane, (but then necessarily 
non-symmetric with respect to the ellipse plane), which would correspond to a different set of attached dislocations, 
most probably of higher energy. 

We apply Eq. 1, rewritten 



to the ellipse, Fig. 35; db tr is perpendicular to the plane of the ellipse, and one easily gets db tr = 2dr; therefore the 
total Burgers vector attached to one side of the ellipse (0 < <fi < ir) is (Kleman and Lavrentovich, 2000) 



Taking dr = do, — an approximation which makes sense (up to second order), since do is so small compared to 
the size a of the ellipse — it is visible that the points M, of polar coordinates {r,<fi}, and N, of polar coordinates 
{r + dr, <fi + d4>}, are on two parallel smectic layers at a distance <i n , and the total Burgers vector attached to a layer is 
equal to 2g?o- There are no dislocations attached to the singular circle of a toric FCD, as the eccentricity e vanishes, 
and Sr = 0. The Burgers vector attached to a layer that is transverse to the perimeter of the ellipse is a constant, 
2e?o- Consider one such dislocation of Burgers vector 2e?o; it spreads outside the ellipse in the shape of a quantized 
dislocation, on both sides of the ellipse, and goes across it in the shape of curved layers, in a manner akin to the 
k = \ case (Fig. 21). The full balance of Burgers vectors has to take into account the Nye's dislocations, including 
those related to the infinitesimal disclinations, 



dh tr = 2t x MN, 



(CI) 




(C2) 



dt = dt 



{cos </>, sin (/)} dc/), 



(C3) 



dt is a vector along the normal of the layer, as required by symmetry (the normal to the layer is an axis of continuous 
rotational symmetry). But there is no a priori reason that the total Burgers vector carried by the Nye's dislocations 
balances topologically the quantized attached dislocations; it is the attachment that provides a balance equation. 
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If the dislocation segments outside the ellipse belong to a plane, they build a grain boundary limited by the ellipse. 
It is a simple matter to check that the misorientation angle is lj = 2 sin -1 e. This grain boundary is a pure tilt 
boundary if it occupies the plane of the ellipse. Notice that u> (the tilt vector, along the direction of the ellipse minor 
axis) and f2 1 (variable along the ellipse) define both the same tilt grain boundary; we are in the case described in 
II.E.l, Eq. 9. 

We know from topological stability theory that the defect geometry of the ellipse as a disclination does not break the 
constancy of its conjugacy class in Hi(V smA) along E. The variation of the representative homotopy class (inside the 
same class of conjugacy) for a circuit embracing the line depends on the quantized dislocations (not the infinitesimal 
dislocations which, as already stated, belong to the identity class of the homotopy group) attached to it, which it also 
embraces. 



APPENDIX D: A few geometrical characteristics of the 3-sphere 

Details on the topics that follow can be found in Sommerville (1967), Montesinos (1987), Coxeter (1998). 

1. Clifford parallels and Hopf fibration 

The trajectory of a point M S S 3 under the action of a right screw, say, x' = xe aq , < a < 2n, is a great circle 
Cright of S 3 . Two right (resp. left) great circles that rotate helically about the same q are equidistant all along their 
length and mutually twisted with a pitch w = 2tt (resp. — 2ir): they form a congruence of so-called right (resp. left) 
Clifford parallels that fill S 3 uniformly. All the great circles of this congruence are equivalent. Because any of those 
great circles can be considered as the axial line of cylindrical double-twist geometry for the other parallel great circles, 
the geometry of Clifford parallels has inspired a non-frustrated model of the Blue Phase (Sethna, 1984). The great 
sphere S 2 defined by the intersection of the hyperplane q\X\ + 92^2 + 93^3 = and S 3 intersects the great circles of 
the congruence orthogonally; this geometry defines S 3 as a fiber bundle of great circles S 1 over a great sphere S 2 , 
the Hopf fibration, see Fig. 36. Notice that great circles are the geodesic lines of S 3 and great spheres the geodesic 
surfaces of S 3 . 



2. Spherical torus 

The axis p and the circle of unit radius in its equator plane (Fig. 29) just play identical roles in the double rotation: 
p is the stereographic projection of a great circle of S 3 , here denoted C p , which is the intersection of the sphere of unit 

X\ X2 X$ ~ 

quaternions x% + x\ + x\ + x\ — 1 with the 2-plane — = — = — ; the equator circle, here denoted C„, is a great 

Pi P2 P3 

circle intersection of the unit quaternion sphere with the 2-plane xq = 0, piX\ + P2X2 + P3X3 = . These two great 
circles are conjugate; they have this remarkable property that the arc of any great circle C g joining any point on C p 
to any point on C p is perpendicular to both and has arc length tt/2. Select any C g and transform it under the action 
of the right screw x' = xe ap , < a < 2tt, say; this yields a set < C g > of great circles, all stretching between C p and 
C p . The surface so obtained is called the spherical torus; it is generated by two families of great circles: < C g > just 
defined, and < C p > made of the Clifford parallels that are parallel to C p (and to C p ). These two families form a 
rectangular network on the spherical torus. 

3. Clifford surfaces 

Consider now the action of a left screw y' — e^ v x' about the same axis p, i.e. y' — e^ p x e a p on Q r i g ht when M(a) 
traverses the geodesic line Crighu and (3 varies in the range < (3 < 2ix: each point on C r i g ht develops into a full 
geodesic trajectory Gieft (a great circle) whose entire set forms a Clifford surface, which is a closed surface with two 
conjugate axes of revolution; these axes are precisely the conjugate geodesic circles of the double rotation p. The 
geometric properties of the Clifford surface depend on the angular distance 1? of G right to C p . This is illustrated 
Fig. 36. 

The particular Clifford surface $ = 7r/4 is a spherical torus, whose axes of revolution are C p and C p . 
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FIG. 36 (a)- C p and C p are conjugate great circles (represented as straight lines in reason of their geodesic character); C r i g ht 
and Ci e ft are equidistant to C p and C p (OA= 1?, ACl = n/2 — All right and left Clifford parallels at the same distance •& 
generate a ruled Clifford surface with two axes of revolution C p and C p ; (b)- representation on the basis S 2 of a Hopf fibration 
of S 3 , with C p and C p as particular fibers. Each fiber is represented by a point on the basis; Clifford surfaces (the small circles 
on S 2 ) # = constant; spherical torus ■& — 7r/4 

APPENDIX E: Geometrical elements related to a great circle in S 3 

1. Great circle defined by two points 

Consider two points M and M' on the 3-sphere S 3 , in quaternion notation u and u', in vector notation u and u'. 
We have 



u • u' = R 2 cos 1? uu' + u'u = 2R 2 cos 1?. 



(El) 



d is the angle between the directions u and u'. 

We wish to define the geometric elements related to the great circle C (centered at the origin {0}) going through u 
and u' . 

We introduce the two quaternions u + u' and u — u'\ they correspond to orthogonal vectors. We now show that 
there are two pure unit quaternions g and a such that 



u + u' oc 1 — ga, 
u' — u oc g + a, 



(E2) 



with the same coefficient of proportionality, a real number. In effect, multiply the second equation by g (on the left) 
in order to get an expression for g, and by a (on the right) in order to get an expression for a, substitute then in the 
first equation, it comes: 



g oc (it' — u)(u' + u) = u'u — uu' 



a oc (u + u)(u' — u) = uu — u'u (E3) 
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where we have used uu, = v!v! (= Ft 2 ). It is easy to check that g and <r are pure quaternions. It remains to renormalize, 
in order to get pure unit quaternions. Using 

\uu' - u'u\ 2 = (uu' - u'u)(u'u - uu!) = 4i? 4 sin 2 ■& 

and the same expression for \u'u — uu'\ , g and o eventually read 

u'u — uu' uu' — u'u /t?a\ 

9= 2i? 2 |sintf| ° = 2i? 2 |sintf| ( ' 

The 2-plane in which lies the great circle C, defined by u, u' , and the origin of the coordinates (which is the center 
of C), contains the directions 1 — go and g + o. It is denoted II e!Cr = {0, 1 — go, g + o}. Any rotation of angle $ that 
leaves invariant this plane leaves invariant the great circle C that contains u and u' ; it transforms any point y £ E 
into y' £ E 4 according to the relation 

y = e 2 e ye2 . 

H g a is the axial plane of the rotation. The plane n^ CT {0, 1 + go, g — o} is completely orthogonal to n e <T . 



2. Great circle defined by the tangent at a point 

Let us look for the great circle going through two very close points u and u + du, when du vanishes continuously; 
one gets: 

duu — udu 1 du _ dil . 
5 = 2i? 2 dtf = 2i? ( ^"" Ms) 
udu — duu 1 du du 
2R 2 d-d 2R y ds ds ' 

where ds = Rdd is the arc element on C. 

Let \x be the unit tangent at u to C, = 1, we have 

8 = 2~R^ ~ ° = 2~R^^ ~ ^ 
Because of the relation of orthogonality /2 • u — 0, we can also write: 

1111 , , 

= ~R^ iU = ° = R Ufl = ~R^ U ' ^ * 

One easily checks that g and o are two pure unit quaternions; g 2 — o 2 = — 1. 



APPENDIX F: The {3,5} and {3,3,5} symmetry groups 



The symmetry group of the {3,3,5} curved crystal is related to the symmetry group Y of the icosahedron, since any 
vertex is the centre of an icosahedron. 



1. The group of the icosahedron {3,5} 

Y is a finite group with 60 elements, represented by rotations about the six 5-fold axes, the ten 3-fold axes, and 
the fifteen 2-fold axes of the icosahedron. The order parameter space of the icosahedron {3,5} is V^s} = SO(3)/Y, 
whose first homotopy group IIi(V{ 3 5 j) ~ Y is the lift of Y in the covering group of SO(3), i.e. in SU(2). Y is also 
known as the binary group < 5, 3, 2 > (Coxeter, 1973); it has twice as many elements as Y. 

The topological theory defect classes of the icosahedron have been investigated by Nelson (1983b). Y is perfect (i.e. 
the commutator subgroup D[Y] = Y), so that in principle all defects can mutually annihilate (Trcbin, 1984). 

The icosahedron vertices themselves can be expressed in terms of quaternions, so that we have all the analytical 
tools to represent the symmetry actions on {3,5}. A caveat: since all the symmetry elements hy £ Y are rotations 
about axes having the centre of the icosahedron as a fixed point, going through the two-sphere S 2 on which {3,5} 
lives, any action on a point x of the icosahedron has the form x h-Y and connotes a disclination, in the Volterra 
process terms. Therefore there are no disvections generated by the point group Y in {3,5}. 
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2. The group of the {3,3,5} polytope 

The 120 elements of Y, in quaternion representation, occupy on S 3 precisely the locations of the 120 vertices of 
a {3,3,5}. Thus the binary icosahedral group constitutes a group isomorphic to the group built by these vertices 
(Coxeter, 1973). The {3,3,5} symmetry point group, (indirect isometries excluded), is 

H = Y x Y (Fl) 

whose lift in G(S 3 ) = S 3 xS 3 is 

H = Y x Y. (F2) 
H has 7200 elements. The displacement of a point x of {3,3,5} under the action of the symmetry group is 

x' = hy xh Y h Y eY, hy G Y (F3) 

(or, taking the conjugate, a point y' = h^y hy), where and hy are not necessarily conjugate, i.e. Eq. F3 composes 
any right action and any left action. 
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